AD-A071  663 


UNCLASSIFIED 


WASHINGTON  UN IV  ST  LOUIS  MO  OEPT  OF  ELECTRICAL  ENGI—ETC  F/8  *0/12 
AN  AC  LARGE  SIGNAL  MOOEL  FOR  THE  GAAS  MESFET.(U) 

AUG  79  A MAOUAR*  F J ROSENBAUM  N0001A-78-C-0256 


/o 3 &0(o 


ii 


TABLE  OF  CONTENTS 

No . Page 

1.  FET  Analysis  Methods 1 

1.1  Introduction . 1 

1 . 2 Shockley  Model 3 

1.2.1  Device  Structure  and  Basic  • 

Properties 3 j 

1.2.2  The  Gradual  Case 7 

1.2.3  The  Expop  Region 13 

1.2.4  Discussion 17 

1.3  Modern  Analysis  Methods 18 

1.3.1  High  Frequency  FET  Structure 18 

1.3.2  General  Approaches 20 

1.3.3  Turner  and  Wilson  Model 23 

1.3.4  Kennedy  and  O'Brien  Simulations 25 

1.3.5  Barnes  Simulations 28 

1.3.6  Yamaguchi  and  Kodera  Simulations 30 

1.3.7  Grubin  and  McHugh  Simulations 33 

1.3.8  Semi-Empirical  Analysis 34 

1.4  Scope  of  This  Report 35 

2.  The  Basic  Approximate  Analytic  Solution 37 

2.1  Introduction 37 

2.2  Problem  Definition 37 

2.3  The  Approximate  Charge  Carrier 

Density  Function 43 

2.4  The  Solution  to  Poisson's  Equation 47 

2.5  Displacement  Current  Expressions 54 

3.  The  Novel  FET  Model  for  The  Active  Region 58 

3.1  Introduction..... 58 

3.2  Case  Classification  Induced  by  Velocity- 

Electric  Field  Curve 59 


iii 


TABLE  OF  CONTENTS 
(continued) 


No . Page 

3.3  Exact  Conduction  Current  Expressions 60 

3.3.1  Source  Conduction  Current 60 

3.3.2  Drain  Conduction  Current 66 

3.3.3  Numerical  Solution  for  V1 68 

3.4  Analytic  Expressions  for  Conduction  Current 

and  70 

3.4.1  Source  Conduction  Current 71 

3.4.2  Drain  Conduction  Current... 74 

3.4.3  Analytic  Solution  for  78 

3.5  Expressions  for  Derivatives 81 

3.6  Calculation  of  g and  g. 84 

3m  a 

3.7  Miscellaneous  Topics 85 

3d.  (0) 

3.7.1  Correction  of  -577 85 

3VSG 

3.7.2  Charge  Accumulation 87 

3.7.3  Smooth  Velocity-Field  Curve 91 

3.7.4  Criteria  for  Operational  Mode 

and  Case 93 

4.  The  Complete  FET  Computer  Model  Structure 100 

4.1  Introduction 100 

4.2  The  Model 100 

4.3  The  Charging  Resistor 106 

4.3.1  Introduction 106 

4.3.2  Charging  Resistance  of  a Distributed 

Network 107 

4.3.3  Mode  A Ill 

4.3.4  Modes  B and  C 113 

115 


4.4  Time  Delay 

4.5  Computer  Package  Structure 


119 


iv 


TABLE  OF  CONTENTS 
(continued) 

» No.  Page 

5.  Simulation  Results 127 

5.1  Introduction 127 

5.2  Number  of  Iterations 129 

5.3  Charge  Correction  Coefficients 132 

5.4  General  Properties 135 

5.5  Scattering  Parameters 14  3 

5.6  Large  Signal  Analysis 151 

5.7  Oscillator  Analysis 155 

6.  Summary  and  Conclusions 176 

7 . Acknowledgments 180 

8.  Appendices 181 

Appendix  8.1  Electric  Field  Expressions 182 

Appendix  8.2  Displacement  Current  Expressions...  187 

Appendix  8.3  Ratio  Function  for  Mode  C 192 

Appendix  8.4  Integral  Function  for  Mode  C 197 

9.  Bibliography 204 


l>  i m J • • - 

- - 


V 


LIST  OF  TABLES 


No.  Page 

4.1  Values  of  INVSAT  for  Each  Mode  and  Case 122 

5.1  Typical  Microwave  FET  Parameters  Used  in 

the  Computation 128 

5.2  Oscillator  Performance 175 


Appendices : 

8.3- 1  The  Ratio  Function 194 

8.3- 2  Coefficient  for  the  Polynomial 

Approximation 196 

8.4- 1  Numerical  Evaluation  of  the  Integral 

for  Jl2=0. 15 199 

8.4- 2  Numerical  Evaluation  of  the  Integral 

for  £2=0.4 199 

8.4- 3  Evaluation  of  the  function  in  8.4-9 

for  i2=0*15 202 

8.4- 4  Evaluation  of  the  function  in  8.4-9 

for  i2=0*4 202 


8.4-4  Evaluation  of  the  function  in  8.4-9 

for  i2=0*4 202 


vi 

LIST  OF  FIGURES 


No . Page 

1.1  FET  Structure 4 

1.2  Current-Voltage  Characteristics  for  FET 9 

1.3  Shockley's  Small  Signal  Equivalent  Network 

for  FET 11 

1.4  Typical  High-Frequency  FET  Structure 19 

1.5  A Three  Dimensional  Graph  of  Electron 

Distribution 29 

2.1  Gate  Active  Region  Considered  in  the  Analysis....  39 

2.2  Velocity-Electric  Field  Curves 41 

2.3  Subdivision  of  FET  Active  Region 44 

2.4  Operational  Modes  of  FET 4 6 

3.1  Partitioning  of  the  Transition  Region  at  the 

Source  Side 65 

3.2  False  Negative  Slope  in  ID~VDS  Curves 89 

3.3  A "Sharp-Corner"  in  I^-V^g  Curves 92 

3.4  "Smooth"  Iq_vds  Curves 94 

4.1  Assembled  FET  Chip  Model 101 

4.2  (a)  FET  Channel  as  a Distributed  Network 

(b)  Discretized  Distributed  Network  Considered 
in  Text  (c)  Equivalent  Network 108 

4.3  FET  Subdivision  for  Charging  Resistor 

Calculation 112 

4.4  Flow  Chart  for  the  FET  Computer  Model 124 

5.1  I-V  Curves  of  the  FET 130 

5.2  g„  Plots 133 

m 

5.3  "Conversation"  between  Operator  and  Computer 

PDP-11/45 


1 


136 


vii 


LIST  OF  FIGURES 


(continued) 


No . Page 

5.4  I-V  Curves  of  the  FET 138 

5.5  FET  Properties 139 

5.6  I-V  Curves  for  the  FET  using  Shockley  Model 140 

5.7  Measured  I-V  Curves  of  the  FET 142 

5.8  Scattering  Parameters  of  the  FET 149 

5.9  Scattering  Parameters  for  Reduced  Bias 152 

5.10  Saturation  Characteristic  of  the  FET 154 

5.11  Oscillator  Network 156 

5.12  Oscillator  Equivalent  Network 158 

5.13  Oscillator  Build-Up 165 

5.14  Oscillator  Output  Voltage  in  Steady  State 167 

5.15  Drain  Current  Waveform 169 

5.16  Drain-to-Source  Voltage  Waveform 170 

5.17  Oscillator  Waveforms  (Second  Design) 172 


Appendices : 


8.3- 1  The  Ratio  Function 195 

8.4- 1  Numerical  Evaluation  of  the  Integral 

for  &2=0.15 200 


viii 

LIST  OF  SYMBOLS 

a 

Epitaxial  layer  thickness 

D 

Diffusion  coefficient 

°L 

Low  field  diffusion  coefficient 

dh 

High  field  diffusion  coefficient 

d 

Transition  region  height 

/ 

dL(x) 

Neutral  channel  height 

E 

Electric  field 

Ec 

Critical  electric  field  for  velocity 
saturation 

^m 

Transconductance 

^d 

Drain  conductance 

I 

con 

Conduction  current 

I 

Terminal  current 

IDIS 

Displacement  current 

XDIF 

Diffusion  current 

Ig 

Gate  current 

XD 

Drain  current 

xs 

Source  current 

J 

Conduction  current  density 

Jtot 

Total  current  density 

4 

Displacement  current  density 

1 

ddif 

Diffusion  current  density 

k 

Boltzmann's  constant 

ix 


e =£oer 


LIST  OF  SYMBOLS  - 
(Continued) 

Gate  length 

Donor  impurity  concentration 

Charge  carrier  concentration 
Electron  Charge 
Resistance 

Absolute  temperature 
time 

charge  carrier  velocity 

charge  carrier  saturated  velocity 

Drain  to  source  voltage 

Source  to  gate  voltage 

Pinchoff  voltage 

MESFET  width 
Dielectric  constant 

Vacuum  dielectric  permittivity 

Dielectric  permittivity 

Debye  length 

Low-field  mobility 

Conductivity 

.4 

Built-in  potential  of  the  Schottky  bairier 
Electric  potential 


1.  FET  ANALYSIS  METHODS 

1.1  INTRODUCTION 

As  suggested  by  the  title,  this  report  presents  a 
novel  practical  AC  large-signal  model  for  the  GaAs  MESFET . 
The  new  model  is  a "true"  AC  large  signal  one  in  the  sense 
that  the  total  terminal  currents  are  given  as  explicit 
functions  of  the  total  terminal  voltages  VgG,  VDg , and 
their  time  derivatives.  Previous  models  and  simulations 
were  either  DC  alone,  or  a DC  analysis  with  a small-signal 
incremental  AC  model  derived  from  it.  Few  authors 
attempted  "true"  AC  simulation  by  numerical  time-dependent 
solution  to  the  device  differential  equations.  This 
approach  is  ordinarily  very  expensive  and  extremely  time- 
consuming  on  a digital  computer. 

The  present  model  is  circuit  design  oriented  and  can 
be  used  to  analyze  and  design  large-signal  components  (os- 
cillators, power  amplifiers,  frequency  multipliers,  etc.). 
As  such  it  is  very  efficient  and  fast  on  a digital  com- 
puter. To  achieve  this  an  approximate  analytic  solution 
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to  the  device  differential  equations  was  derived.  Since 
the  model  is  derived  from  basic  principles,  it  is  directly 
dependent  on  device  geometry  and  the  semiconductor  pro- 
perties. 

Our  analysis  is  based  on  the  general  method  suggested 

* l 

by  Yamaguchi  and  Kodera  II]  . Their  t>asic  approach  was 
extended  and  improved  in  two  ways:  i]  by  extending  it  into 
a "true"  AC  large  signal  model  2)  by  making  it  almost  com- 
pletely analytic  by  solving  for  an  unknown  parameter  in 
their  equations  in  an  analytic  fashion  instead  of  their 
numerical  iterative  method.  This  improves  dramatically  the 
computation  time,  which  is  0.05-0.1  seconds  per  point  in 
our  model  compared  to  1 second  for  the  iterative  method. 

The  model  has  only  modest  memory  requirements  ( 50K  words 
on  IBM  3601 . 

This  chapter  includes  a review  of  the  main  FET  analy- 
sis methods  presented  in  the  literature.  The  review  starts 
with  the  classical  Shockley  approach  and  ends  with  the  more 
recent  modern  approaches.  The  various  advantages,  disad- 
vantages and  limitation  of  each  method  are  discussed. 
Special  attention  is  given  to  the  validity  of  the  analysis 
methods  for  modern  high-frequency  GaAs  devices.  These 
devices  are  characterized  by  very  short  gates  (micron  or 

* The  numbers  in  parenthesis  in  the  text  indicate 
references  in  the  Bibliography. 
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sub-micron) , thus  yielding  relatively  small  length  to 
height  ratios  for  which  the  Shockley  gradual  approximation 
may  be  invalid.  Also,  the  electric  field  in  the  device  is 
normally  very  high,  thus  making  the  device  properties 
dominated  by  saturation  of  the  velocity  of  the  charge  car- 
riers. 

1.2  SHOCKLEY  MODEL 

The  so-called  field-effect  transistor  was  first  intro- 
duced by  Shockley  [2]  in  1952.  Shockley  suggested  the  con- 
struction of  the  new  device  and  developed  a detailed  analy- 
sis of  it.  This  was  basically  a DC  analysis,  and  as  a 
result  the  I-V  curves  of  the  new  device  were  deduced. 
Shockley  also  derived  a simplified  incremental  small  signal 
AC  model.  The  basic  limitations  of  this  analysis  are: 

1)  charge  carriers  are  assumed  to  have  constant  mobility  - 
no  consideration  of  velocity  saturation.  2)  one  dimen- 
sional approach  - only  one  component  of  the  electric  field 
is  considered  (the  quadrature  component  is  assumed  zero) . 

3)  abrupt  transition  assumed  between  the  conducting  chan- 
nel and  the  depletion  region.  To  avoid  confusion  in  the 
review  of  Shockley's  model  below,  his  notation  is  used.. 
1.2.1  Device  Structure  And  Basic  Properties 

The  device  introduced  by  Shockley  is  illustrated  in 
Figure  1-1.  It  is  a three-terminal  device  consisting  of  a 
p-type  layer  sandwiched  between  two  layers  of  n+  type  mat- 
erial. The  electric  current  is  composed  of  holes  flowing 
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in  the  x-direction  in  the  p layer.  The  terminal  contacts 
are  formed  by  p+  material.  In  normal  operation  reverse 
bias  is  applied  to  the  two  junctions,  which  form  the  space- 
charge  regions  in  which  the  carrier  concentration  is  prac- 
tically zero.  So  the  carriers  flow  in  a channel  of  p-type 
material  bounded  by  two  space-charge  regions. 

As  a first  step  assume  that  the  applied  voltages 
across  the  two  junctions  are  equal.  This  means  that  no 
current  is  in  the  channel  and  so  the  channel  is  uniform. 

The  charge,  electric  field,  and  potential  distributions 
in  the  space  charge  region  are  derived  from  solution  to  the 
one  dimensional  Poisson  Equation: 


dy 

which  yields: 


dEy  = qNa 


(Ey=0  at  y=b) 


qN. 


Cy-b) 


Ey  = - 

c 

qW  2 2 

v = TTl  (y"b)  " (a~b)  J 


(1-1) 

(1-2) 

(1-3) 


The  constant  of  integration  was  chosen  to  make  V=0  at 
y=a  (gate  as  the  reference  voltage) . 

The  potential  in  the  channel  is: 

qN  2 

V (y=b)  = - (a-b)  (1-4) 

To  avoid  using  negative  numbers  in  the  derivation,  Shockley 
defined: 

2 

W = - V (y=b)  = II- (b/a) ] Wq  (1-5) 


y 
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qNa  2 

Wo  - IT  3 U_6) 

WQ  is  the  voltage  required  to  make  the  space  charge  pene- 
trate the  entire  p-region.  He  termed  this  as  the  “pinch- 
off"  voltage. 

Now  assume  there  is  a small  difference  between  the 
voltages  applied  across  the  two  junctions.  In  that  case 
the  channel  can  be  assumed  almost  uniform  as  before.  The 
channel  conductivity  is  qN.^,  where  ^ is  the  hole  mobi- 
lity (assumed  constant) . Thus,  the  current  per  unit  width 
in  the  channel  for  an  applied  Ex  is: 

I = 2bqNayQEx  = g(W)  Ex  (1-7) 

gCW)  = 2qNayob  (1-8) 

g(.Wl  is  the  conductance  of  a unit  square  of  the  layer  2b 
thick.  So,  the  total  conductance  per  unit  width  is 

G = g(W)/L  (1-9) 

g(.W)  is  a function  of  the  reverse  bias  W.  Using  Equations 
1-9,  1-5  one  gets: 

1/2 

g(.W)  = II-  ($-)  ] gQ  (1-10) 

O 

where  g m 2qN  u a 
o aHo 

Equation  1-10  clearly  demonstrates  the  basic  principle 
of  FET  operation,  namely,  a voltage  controlled  conductance. 
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Shockley  suggested  the  following  names  to  the  elec- 
trodes : 

source  for  electrode  1 through  which  the  carriers  flow  into 
the  channel. 

drain  for  electrode  2 through  which  the  carriers  flow  out 
of  the  channel. 

gate  for  the  control  electrode  3. 

These  names  were  accepted  and  are  commonly  used. 

1.2.2  The  Gradual  Case 

When  the  applied  source  and  drain  voltages  are  not 
equal  the  channel  is  not  uniform,  and  Equation  1-5  is  in- 
valid, since  it  was  derived  for  the  one-dimensional  case. 
32V 

However,  if  -r— =■  in  the  depletion  region  is  very  small  com- 
3X2 

pared  to  qNa/e,  then  the  one  dimensional  approximation  is 
valid  and  so  is  Equation  1-5.  For  the  above  to  be  true 
the  conditions  along  the  channel  should  vary  gradually, 
thus  this  assumption  is  called  by  Shockley  the  "gradual 
approximation".  The  gradual  approximation  holds  if  g^<<l. 
In  the  gradual  case  the  electric  field  in  the  space  charge 
region  is  assumed  to  be  in  the  y direction,  while  in  the 
channel  it  is  assumed  to  be  in  the  x direction. 

The  basic  equation  for  the  gradual  case  is: 

1 - g(w)  (1-11) 

which  is  equivalent  to  Equation  1-7. 
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Rearrange  and  integrate  Equation  1-11  and  get: 


g (W)  dW 


(1-12) 


■t  r 

VIJ 


g (W)  dW 


(1-13) 


Define  the  function: 


J(W) 


1/2 

= / g (W)  dW  = g WI1-  §(§-)  ] 

-A  o 


(1-14) 


Then: 


I = IJC Wd)  - J (Wg ) ] /L 


= U(Vg-Vd)  - J(Vg-Vs)]/L 


(1-15) 


A plot  of  I vs.  |Vd|  with  Vg  as  a parameter  and 

V = 0 is  given  in  Figure  1-2.  Equation  1-15  is  not  valid 
s 

for  > Wq  because  pinchoff  occurs  and  the  gradual  approxi- 
mation fails.  The  continuation  of  the  curves  in 
Figure  1-2  is  based  on  the  analysis  developed  in  the  next 
section.  The  critical  condition  for  pinchoff  (Wd  = WQ)  is 
indicated  by  the  dashed  line.  The  equation  for  this  curve 
is 


1 

y 


I = [J(Wo)  - J(Wo  -|Vd[)]/L 


(1-16) 
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In  the  case  of  a small-signal  AC  superimposed  on  the  DC 
it  is  possible  to  calculate  the  small  signal  drain  current 
as  a function  of  the  small  signal  voltage  by  taking  the 
time  derivative  of  Equation  1-15.  The  result  is: 

Ad  = -i  = Gdg[(VVs’  + (vd-vg)/ydg]  (1"17) 

where : 

Gdg  = 9<Ws)/L  (1-18) 

Rd  = L/g(Wd)  (1-19) 

^dg  " Gdg  Rd  (1-20) 

Equation  1-17  suggests  a small  signal  AC  equivalent 
circuit  for  the  FET  as  shown  in  Figure  1-3.  This  is  an 
over-simplified  equivalent  circuit  and  does  not  include  any 
parasitic  elements. 

Shockley  also  derives  the  channel  shape  (i.e.  b vs.  X) . 


From  Equation  1-5 

-(a-b)qN 

dw  = db  (1-21) 

and  insertion  into  Equation  1-11  leads  to 

IdX  = -IQa(l-u)u  du  (1-22) 

where : 

u = b/a  (1-23) 

Xo  - 90  Eq  (1-24) 

Eo  " 2Wo/a  d-25) 
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Integrating  Equation  1-22  yields 
al  2 3 

X = - (-^)  tj  - j ] U-26) 

Equation  1-26  gives  the  shape  of  the  channel.  Once  it 
is  known  the  voltage  along  the  channel  is  readily  calcu- 
lated using  Equation  1-5. 

Using  Equation  1-26  it  is  possible  to  derive  the  con- 
dition for  the  validity  of  the  gradual  approximation:  over 
a distance  a the  fractional  change  in  channel  width  is 
small: 


a db 
b dX 


(l-u)u  I 


< 1 


(1-27) 


From  the  above  equation  it  is  obvious  that  the  gradual 
approximation  fails  for  u=0  and  u=l.  For  u=l  (full  chan- 
nel) the  failure  is  of  little  importance  since  the  electric 
field  is  small  and  the  conductance  large.  Also  u=l  does  not 
happen  in  practice,  since  there  is  always  some  residual 
depletion  region.  For  u=0  the  electric  field  is  high,  and 
this  case  is  always  approached  when  > WQ.  For  u 0 
Equation  1-27  becomes: 


I 

u2I 


< 1 


(1-28) 


Also  Equation  1-26  becomes: 

2 

I au^ 
o 


IT 


1 


-x 


(1-29) 


-13- 


constriction  at  y=0.  Since  the  charge  due  to  these  holes 
is  very  small  compared  to  the  space  charge  in  the  depletion 
region,  that  part  of  the  electrostatic  potential  which  is 
due  to  charges  is  well  approximated  by 

2 

VQ(x,y)  = - WQ[1  -(£)  ] (1-31) 

which  satisfied  Poisson's  equation  and  goes  to  zero 
at  y = ±a  and  to  -Wq  at  y=0.  This  potential  is  independent 
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of  X and  produces  no  field  which  causes  the  holes  to  flow. 
The  field  is  produced  by  the  drain  voltage.  The  poten- 
tial distribution  due  to  drain  voltage  is  the  solution  of 
Laplace's  equation  with  the  boundary  condition  that  it 
vanishes  at  y = ±a.  The  general  solution  is: 

cr> 

T(x,y)  = ^ A exp[ir  (2n+l)X/2a]  cos  [ir  (2n+l)y/2a]  (1-32) 

n=0  n 

The  zero  order  term  (n=0)  is  usually  an  adequate  ap- 
proximation to  the  potential  function,  due  to  the  exponen- 
tial decay  (the  decay  per  each  interval  of  a is  e-1T//2  =0.21 

3ir 

“ o 

for  n=0,  and  e = 0.009  for  n=l) . So,  in  the  vicinity  of 
the  expop  the  potential  is  approximated  by 

V (x, y)  = VQ(x,y)  - ATQ(x,y)  (1-33) 

where  VQ(x,y)  is  given  by  Equation  1-31  and: 

T0(x,y)  = exp  (irX/2a)  cos(7ry/2a)  (1-34) 

The  only  problem  left  is  to  determine  the  coefficient 
A.  Shockley  determines  it  by  a graphical  method.  He  draws 
a graph  of  the  potential  versus  X.  For  x < - ^ he  uses 
Equations  1-26  and  1-5.  For  X > a,  Equation  1-33  is  used 
with  several  values  of  A.  By  observing  the  graphs 
(Figure  12  in  [2]),  Shockley  chose  that  value  for  A,  which 
gave  a "smooth"  transition  from  the  gradual  to  the  exponen- 
tial solution.  This  value  of  A is 


i 
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1/2 

A = Wo(^-)  (1-35) 

o 

Using  Equation  1-33  for  the  drain  V (Xd  ) = - it  is 
possible  to  solve  for  X^,  the  distance  between  the  drain 
and  expop : 

xd  - ir  ln  (<W  (2Vi)1/2/wo>  (1'36) 

The  logarithmic  dependence  of  Xd  on  drain  voltage  and  cur- 
rent suggests  that  Xd  is  insensitive  to  these  values.  As 
an  example,  for  IQ/I  = 18,  X^  changes  from  1.15a  at 
Wd  = 2Wq  to  2a  at  = 5Wq.  In  view  of  the  above  it  is 
possible  to  explain  the  continuation  of  the  curves  in 
Figure  1-2  for  | | > Wq.  Equation  1-11  holds  only  for  the 
gradual  portion  of  the  channel,  and  the  integration  must  be 
taken  between  and  Wq.  Accordingly,  X varies  from  zero 
to  a distance  Lg  up  to  the  expop,  where: 

Ls  = L - Xd  (1-37) 

The  result  of  the  integration  is 

I = [J(Wq)  - J(Vg-Vs)]/Ls  (1-38) 

Therefore,  for  Wd  > W the  current  is  almost  indepen- 
dent of  Vd,  thus  explaining  the  horizontal  line  continua- 
tion in  Figure  1-2.  There  is  a weak  dependence  on  Vd  via 
Lg.  In  fact,  as  |vd|  increases,  Xd  increases,  Lg  decreases, 
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and  I increases.  Thus,  the  curves  are  not  horizontal  but 
have  a slight  positive  slope. 

The  corresponding  small  signal  AC  equation  can  be 
derived  in  a way  similar  to  that  in  the  previous  section. 
Starting  with  equation: 


I (L  - Xd)  = 


' P 

W J 


g (W)  dW 


(1-39) 


take  appropriate  derivatives  and  get: 


i = I 

where : 


-li  - 2*  + v (-i i_ 

Exs  Exd  g Exd  Exs 


Exs  - J/g(Ws> 


E , = -4-  (W.  - W ) 
xd  2a  d o 


/(L-v  f> 


(1-40) 

(1-41) 

(1-42) 


Comparing  Equation  1-40  to  Equation  1-17  one  can  conclude: 
E. 


dg 


3dg 


xd 


xs 


- 1 


(L-Xd+  f) (1  + pdg) 
1 


Gdg  ^dg 


(1-43) 

(1-44) 

(1-45) 


So,  the  equivalent  circuit  is  the  same  as  in 


Figure  1-3  but  the  component  values  are  calculated  from 
Equations  1-43  to  1-45. 
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1.2.4  Discussion 

The  Shockley  equations  as  derived  above  are  generally 
valid  and  useful  for  low  frequency  devices.  These  are 
usually  made  of  silicon  or  germanium,  their  gate  length  is 
typically  several  microns,  and  their  thickness  is  large 
compared  to  Debye  length: 


Under  these  conditions  the  electric  field  in  the  de- 
vice is  low  enough  so  that  no  velocity  saturation  exists. 
Also,  the  transition  region  between  the  channel  and  the 
depletion  region  (which  is  in  the  order  of  a Debye  length) 
is  negligible,  so  that  an  abrupt  transition  assumption  is 
reasonable.  Also,  since  the  device  length  to  height  ratio 
is  large,  the  gradual  approximation  holds  and  a one 
dimensional  approach  is  acceptable. 

Modern  high  frequency  devices  are  made  of  GaAs  or 
InP  materials  with  low  threshold  field  for  velocity  satur- 
ation. The  gate  length  is  on  the  order  of  1 ym  and  even 
less,  and  the  height  is  usually  less  than  0.5  ym.  For 
these  devices  the  electric  field  is  usually  much  higher 
than  the  threshold  for  velocity  saturation,  thus  the  device 
properties  are  strongly  affected  by  velocity  saturation 
effects.  Also,  the  length  to  height  ratio  is  not  large,  so 
the  gradual  approximation  is  generally  invalid,  and  a two- 
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dimensional  approach  is  necessary.  Also,  since  the  device 
height  is  not  very  large  compared  to  the  Debye  length,  it 
is  necessary  to  take  into  account  a non-zero  transition 
region  between  the  channel  and  the  depletion  region. 

In  addition  to  the  above,  for  the  high  frequency 
devices  it  is  necessary  to  calculate  the  displacement  cur- 
rents and  the  capacitive  effects  associated  with  them,  as 
well  as  other  high  frequency  effects  such  as  transit  time 
through  the  device  and  a depletion  region  charging  time 
constant.  Shockley  did  not  consider  these  effects. 

In  the  last  decade,  when  the  advantages  of  the  FET 
for  high  frequencies  were  realized,  many  researchers 
faced  the  problem  of  a lack  of  a valid  model  for  the  high 
frequency  FET.  Therefore,  much  work  has  been  done  during 
the  last  few  years  to  develop  a useful  model.  The  next 
paragraphs  of  this  chapter  contain  descriptions  of  these 
efforts,  as  well  as  their  advantages  and  disadvantages. 

1.3  MODERN  ANALYSIS  METHODS 
1.3.1  High  Frequency  FET  Structure 

The  modern  high  frequency  FETs  are  constructed  in  a 
different  way  than  the  low  frequency  devices  proposed  by 
Shockley  (Figure  1-1) . They  are  usually  fabricated  in  a 
planar  structure  as  shown  in  Figure  1-4.  The  device  is 
built  on  a semi-insulating  substrate  (GaAs  or  InP) . On 
top  of  this  substrate  an  epitaxial  layer  of  n-type  material 
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is  grown.  A typical  doping  level  is  1017  cm-3.  The  elec- 
trode metalization  is  deposited  on  top  of  the  epitaxial 
layer,  thus  achieving  a planar  device  which  is  easier  to 
process  compared  to  the  Shockley  structure.  Note,  that  the 
device  in  Figure  1-4  is  not  symmetrical  as  in  the  Shockley 
device  (only  one  gate  electrode) . Also  it  has  a Shottky- 
barrier  junction  rather  than  the  p-n  junction  in  the 
Shockley  device.  The  high  frequency  devices  are  always 
n-channel  to  utilize  the  much  higher  mobility  of  elec- 
trons compared  to  holes.  Although  most  of  the  modern 
analysis  methods  described  below  refer  to  the  high 
frequency  structure  (Figure  1-4) , they  are  equally 
applicable  to  the  Shockley  device. 

1.3.2  General  Approaches 

Some  of  the  important  papers  dealing  with  high  fre- 
quency FF;T  modeling  are  listed  in  the  reference  list 
( 1 1 1 r [3]  — [21] ) . Before  describing  the  special  features 
of  some  of  these  models  it  is  useful  to  present  the  general 
approaches  used  by  the  authors.  The  various  simulations 
belong  in  one  of  two  major  categories:  a)  extension  and 
correction  of  Shockley  model,  b)  numerical  solution  to 
the  device  differential  equations.  References  [9],  [17]  do 
not  belong  to  either  of  the  above  approaches.  They  deal 
with  special  methods  to  be  described  below. 
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References  [3]  , [4]  , [6]  , [8]  , [10]  , [12]  , [13]  , [18]  use 
approach  (a),  that  is,  they  try  to  improve  Shockley's  model 
and  adapt  it  to  high-frequency  devices.  The  general  char- 
acteristics of  these  methods  are: 

1)  basically  a one  dimensional  approach  (with  some  two 
dimensional  analysis  for  the  high  field  region  like 
Shockley) . 

2)  assumption  of  an  abrupt  transition  from  channel  to 
depletion  region. 

3)  incorporation  of  velocity  saturation  effects  on  the 
I-V  curve  (basically  follow  Turner  and  Wilson  [3]  approach) . 

4)  extension  of  Shockley's  analysis  for  the  high  field 
drain  region. 

From  this  group  of  papers  the  basic  paper  of  Turner 
and  Wilson  [3]  is  discussed  in  some  detail  later  in  this 
chapter. 

References  [1]  , [5]  , [7] , [11] , [14] - [16] , [19]  - [21]  use 
approach  (b) , namely,  a numerical  solution  to  the  device 
differential  equations.  The  characteristics  of  this  ap- 
proach are: 

1)  Start  from  basic  principles  and  write  the  exact 
partial  differential  equations  describing  the  device 
behavior  (two-dimensional) . 

2)  State  the  boundary  conditions  applicable  to  the 
geometric  structure. 
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3)  Solve  the  equations  under  the  stated  boundary  con- 
ditions. 

Since  the  above  equations  are  usually  quite  compli- 
cated, no  analytic  solution  has  been  published,  and  it  seems 
unlikely  that  such  a solution  will  be  found.  However,  with 
existing  powerful  digital  computers  and  the  advanced  state 
of  numerical  analysis  it  is  quite  feasible  to  solve  the 
equations  numerically.  This  has  been  done  quite  exten- 
sively in  the  last  few  years,  thus  yielding  accurate 
characterization  of  FETs.  Since  no  major  assumption  or 
approximation  is  done  in  formulating  the  equations,  the 
simulation  is  accurate  (provided,  of  course,  a good  and 
stable  numerical  method  is  employed) . The  greatest  disad- 
vantage of  this  approach  is  the  large  computation  time  (and 
large  expense  associated  with  it) . The  new  model  to  be 
described  in  this  report  removes  this  disadvantage. 

The  FET  differential  equations  are  (see,  for  example, 
Yamaguchi  and  Kodera  [1],  Kennedy  and  O'Brien  [5]  or 
Reiser  [7] ) : 


V24;  = 
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Equation  1-47  is  Poisson's  equation.  Equation  1-48  is  the 
current  continuity  equation.  Equation  1-49  defines  the  con- 
duction current  (drift  and  diffusion) , and  Equation  1-50 
defines  the  total  current  (conduction  and  displacement) . 

The  numerical  solution  to  the  above  equations  is  per- 
formed using  the  following  steps: 

1)  define  a grid  of  points  in  the  two  dimensional 
region. 

2)  approximate  derivatives  by  differences. 

3)  use  (2)  in  the  above  equations  for  each  point  of 
the  grid. 

4)  get  a set  of  linear  algebraic  equations 

(the  number  of  equations  is  equal  to  the  number 
of  grid  points) . 

5)  solve  the  above  set  using  standard  computer  sub- 
routines. 

Some  of  the  important  results  reported  in  the  litera- 
ture using  this  approach  are  described  in  the  following 
sections. 

1.3.3  Turner  And  Wilson  Model 

A good  example  of  an  approach  attempting  to  correct 
the  classical  Shockley  model  is  the  Turner  and  Wilson  model 
[3] . Many  of  the  recent  publications  are  basically  an 
improvement  of  this  model.  The  main  contribution  reported 
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in  [3]  is  the  incorporation  of  velocity  saturation  effects 
on  device  properties.  Although  Turner  and  Wilson  consider 
the  non- symmetrical  device  (Figure  1-4) , the  analysis  is 
applicable  to  the  Shockley  device  too. 

The  basic  idea  is  to  use  the  Shockley  equations  for 
low  VQS,  where  the  electric  field  is  below  the  critical 
field  for  velocity  saturation.  At  the  onset  of  velocity 
saturation  the  current  through  the  device  can  be  expressed 
both  using  the  Shockley  expression  and  also  as: 


**  “ Xo  (1-u) 

Xo  = ^NDV3Wa 
where 


(V, 
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" VGS)/VP 


(1-51) 

(1-52) 


Equating  the  two  expressions  gives  an  equation  for 
the  onset  of  velocity  saturation: 

3Ec  3u2  - 2u3 




u at 


(1-53) 


Thus  the  I-V  curves  of  the  device  are  calculated  by 
using  the  Shockley  equation  from  VD£,  = 0 to  the  point 
determined  by  the  critical  u,  at  which  point  the  current 
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remains  constant  (independent  of  VDg) . This  is  demon- 
strated in  Figure  3 of  [3] . As  argued  in  the  Shockley 
model,  here  too  the  current  does  not  actually  remain  con- 
stant, but  has  a slight  positive  slope.  Beside  the 
correction  to  the  I-V  curve.  Turner  and  Wilson  also  derive 
corrections  to  the  transconductance  and  input  capacitance 
as  function  of  the  parameter  u.  For  V^s  = 0 the  equations 
are: 


g 


m 


(3-2u) u 
6-9u+4u2 


(1-54) 
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(1-55) 


gmo  is  the  transconductance  calculated  by  Shockley's 
equation. 

This  model  has  the  same  basic  limitations  as  the 
Shockley  model,  except  that  the  effect  of  velocity  satura- 
tion on  some  of  the  device  properties  is  considered. 

1.3.4  Kennedy  And  O'Brien  Simulations 

One  of  the  first  attempts  to  solve  the  FET  differen- 
tial equations  numerically,  and  thus  predict  the  device 
behavior,  was  made  by  Kennedy  and  O'Brien  [5].  They 
present  an  extensive  amount  of  data  generated  by  the 
computer  simulation.  Some  of  their  important  results  are 
described  below.  All  their  simulations  are  done  for 
silicon  devices,  thus  avoiding  the  negative  differential 
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mobility  of  GaAs  with  the  charge  accumulation  acompanying  : 

it.  The  device  analyzed  is  similar  to  Shockley's  { 

structure  (Figure  1-1)  . 

The  data  are  presented  in  three  types  of  Figures:  map 

of  charge  carriers  distribution,  map  of  electric  potential 

distribution,  and  I_  vs.  V_c  curve  with  constant.  The 

D DS  sg 

) 

authors  compare  their  data  to  data  derived  from  Shockley's 
equations  and  also  check  some  of  Shockley's  assumptions. 

The  important  observations  and  conclusions  are: 

1)  The  transition  from  the  neutral  channel  to  the 

depletion  region  is  gradual  and  not  abrupt  as  Shockley  f 

assumed.  There  exists  a very  clear  transition  region,  ] 

which  may  be  important  for  thin  devices. 

2)  In  the  case  of  field-dependent  mobility,  namely, 
velocity  saturation,  there  is  some  charge  accumulation  in 
the  device  in  the  high  field  region  near  the  drain, 
accompanied  with  a depletion  region  further  toward  the 
drain  to  preserve  charge  neutrality.  The  explanation 
suggested  by  the  authors  to  this  phenomenon  is  simply 
current  continuity  in  the  device.  Toward  the  drain  the 
channel  narrows.  In  the  case  of  constant  mobility  the 
current  is  kept  constant  by  the  increase  in  electric  field. 

However,  when  the  velocity  is  saturated,  the  only  way  to 
keep  the  current  constant  in  spite  of  the  narrowing  channel 
is  by  charge  accumulation.  It  is  worthwhile  mentioning  here 
that  there  is  another  mechanism  to  keep  the  current 

I 
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constant.  Since  the  problem  is  two-dimensional,  the 

velocity  vector  can  rotate.  Thus  the  x component  of  the  , 

velocity  can  increase  inspite  of  the  fact  that  the  magnitude 
of  the  velocity  vector  is  constant.  This  mechanism  was 
suggested  by  Yamaguchi  and  Kodera  [1]  and  is  discussed 

/ 

further  in  this  chapter. 

j 

3)  In  the  case  of  velocity  saturation  the  drain 
current  vs.  VDg  curve  saturates  for  voltages  below  the 
Shockley  "pinchoff  voltage".  Thus,  the  main  reason  for 
current  saturation  is  the  charge  carrier  velocity 
saturation  and  not  pinchoff  effects.  In  fact,  by  observing 
Kennedy  and  O'Brien's  charge  distribution  maps  it  is  seen 
that  for  many  devices  the  channel  is  never  pinched  off 
even  for  very  large  V'Dg. 

4)  By  observing  the  electric  potential  maps  one  con- 
cludes that  the  one  dimensional  approach  of  Shockley  is  not 
valid  for  gate  length  to  epitaxial  layer  thickness  ratio 
smaller  than  5.  Even  when  the  ratio  is  5 the  constant 
potential  contours  are  not  parallel  to  the  device  axes, 
thus  suggesting  that  the  electric  field  is  two-dimensional. 

Thus,  the  extensive  computer  simulations  of  Kennedy 
and  O'Brien  yield  useful  data  which  help  in  the  under- 
standing of  some  of  the  physical  properties  of  the  FET. 

These  simulations  are  also  very  useful  for  checking  the 
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approximations  made  by  many  authors.  Since  1970  much  more  j 

work  has  beendone  on  numerical  simulations,  including  those  j 

for  GaAs  devices.  Some  of  this  work  is  described  in  the 

next  three  sections.  ^ 

i 

1.3.5  Barnes'  Simulations 

/ 

One  of  the  recent  extensive  numerical  simulations  of 
a GaAs  MESFET  was  performed  by  Barnes  [14] . Like  most  of 
the  previous  simulations  this  one  is  based  on  the  numeri- 
cal solution  of  the  partial  differential  equations  of  the 
device.  The  report  contains  a large  number  of  graphs 
displaying  the  calculated  characteristics  of  several 
devices.  The  author  uses  three  dimensional  graphics  to 
display  functions  like  charge  carrier  concentration, 
voltage,  electric  field  etc.  vs.  both  x and  y 
(as  an  example,  see  Figure  1-5) . 

Most  of  the  data  generated  is  for  DC  conditions,  thus 
allowing  the  calculation  of  the  static  I-V  curves.  The 
author  uses  his  simulation  program  to  investigate  the 
physical  properties  of  the  device  and  their  dependence  on 
various  parameters  such  as  doping  level,  geometry,  non- 
uniform  doping,  substrate  effects,  etc.  It  is  worthwhile 
noting  that  the  simulations  are  performed  for  doping  levels 
in  the  range  10'1'5-10'*'6  cm  ^ . No  simulation  is  presented 
for  more  practical  devices  in  the  range  lO^7  cm”^.  The 
reason  for  this  is  probably  the  much  larger  number  of  grid 
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points  necessary  in  the  high  doping  case  to  obtain  a stable 
numerical  solution.  Obviously,  when  the  number  of  grid 
points  is  too  large  the  simulation  is  not  practical  due  to 
both  excessive  memory  requirements  and  extremely  large 
computation  time. 

Barnes'  simulation  is  a good  example  of  the  usefulness 
of  numerical  methods  for  the  purpose  of  FET  investigation. 
Although  it  is  expensive  and  time  consuming,  the  data  gen- 
erated is  accurate,  and  useful  in  the  understanding  of  the 
device  properties. 

1.3.6  Yamaguchi  And  Kodera  Simulations 

Another  interesting  simulation  is  the  one  developed  by 
Yamaguchi  and  Kodera  [1]  and  its  extension  with  S.  Asai 
[15].  Paper  [1]  serves  also  as  the  basis  for  the  practical 
model  described  in  this  thesis.  This  aspect  is  presented 
in  detail  in  chapter  2.  Here  only  some  results  of  their 
numerical  simulation  are  presented.  The  simulation  in  [1] 
neglects  the  negative  differential  mobility  of  GaAs , and 
the  assumption  is  made  that: 

yQE  / for  E < E critical 
v (E)  = ) (1-56) 

vg  I for  E > E critical 


i 


Some  results  of  the  numerical  simulation  are  given  in  the 
form  of  charge-carrier  density  plots,  I-V  curves,  etc.  An 
interesting  result  emerging  from  their  simulation  is  an 
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explanation  of  current  mechanism  in  the  device  for  the  case 

of  carrier  velocity  saturation  on  the  drain  side,  but  no 

l-  saturation  on  the  source  side.  Kennedy  and  O'Brien  [5] 

helve  already  shown  that  to  maintain  current  continuity 

there  is  a possibility  of  charge  accumulation  on  the  drain 

side.  Y^maguchi  and  Kodera  have  pointed  out  another 

mechanism:  rotation  of  the  velocity  rector.  Even  though 

the  velocity  is  constant  at  v„  on  the-,  drain  side,  as  V„ 

s DS 

increases  the  velocity  vector  rotates  such  as  to  increase 
the  longitudinal  component  of  the  velocity.  This  increases 
the  drain  current.  In  practice  both  mechanism  exist.  For 
thick  devices  (compared  to  about  6 times  Debye  length)  the 
charge  accumulation  effect  is  dominant,  while  for  thin 
devices  the  velocity  vector  rotation  is  dominant. 

The  aspect  of  instability  in  FETs  due  to  Gunn  domain 
formation  is  analyzed  in  [15]  . As  in  [1]  here  too  a 
numerical  simulation  is  performed,  however,  here  the 
negative  differential  mobility  is  considered.  Generally, 
the  simulation  shows  that  instability  is  possible  in 
relatively  thick  devices.  The  simulation  actually  shows  the 
formation  of  a domain  traveling  from  gate  to  drain.  As  the 
device  is  made  thinner  the  traveling  domain  becomes  a 
static  domain  (appreciable  charge  accumulation  near  the 
drain) . For  even  thinner  devices  the  domain  disappears, 
and  the  device  properties  are  practically  like  those  of  a 
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silicon  device  (the  authors  call  this  a "pentode-like" 
operation) . Based  on  their  numerical  analysis,  the  authors 
developed  an  approximate  analytic  criterion  for  determining 
the  stability  of  a given  device.  Define: 

. I_2s_  )1/2 

( 9HD  (*B+VSG','Eth  Wf  f1'571 

4>b  - built  in  potential  of  the  Schottky 

barrier  junction. 

- critical  field  for  negative 
differential  mobility. 

VSQ  - source  to  gate  voltage. 

The  device  is  absolutely  stable  (no  Gunn  domains)  if  the 
device  thickness  is  less  than  W^.  For  device  thicknesses 
between  and  2Wfc^  a static  domain  exists  (this  shows  up 
in  the  I-V  curves  as  a negative  resistance  region) . For 
device  thicknesses  over  2Wfc^,  travelling  domain  exists. 

The  above  criterion  is  very  useful  and  enables  the 
design  of  stable  devices.  It  also  sets  limits  on  the 
validity  of  simulations  which  neglect  the  negative 
differential  mobility  of  GaAs. 

Shur  and  Eastman  [22]  have  attempted  to  improve  the 
original  Yamaguchi  and  Kodera  model  [1]  by  considering  a 
static  domain.  They  have  subdivided  the  device  into  two 
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regions:  1)  under  the  gate  - simulated  by  the  Yamaguchi 

and  Kodera  model;  and  2)  gate  to  drain  region  - simulated  by 
a static  domain.  There  is  some  inaccuracy . involved  in  this 
assumption  since  usually  the  domain  penetrates  substantially 
into  the  gate  region. 

1.3.7  Grubin  And  McHugh  Simulations 

An  interesting  numerical  analysis  was  presented  by 

Grubin  [16]  and  Grubin  and  McHugh  ([20] , [21] ) . The 

emphasis  in  their  simulation  is  the  investigation  of  the 

instabilities  due  to  Gunn  domain  formation.  The  analysis 

is  a "true"  time  dependent  AC  large  signal  analysis.  The 

problem  is  three  dimensional  (two  spatial  coordinates  and 

time) . As  mentioned  above  numerical  simulation  of  this 

type  is  time  consuming  and  expensive.  To  reduce  time  and 
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cost  the  analysis  was  made  for  10  cm  doping  level 
(instead  of  the  realistic  value  around  10  cm  ) , thus 
requiring  many  fewer  grid  points  for  the  numerical  analysis. 

The  analysis  performed  is  of  the  transient  response 
type.  The  gate  to  source  and  drain  to  source  voltages  are 
changed  rapidly  from  one  DC  value  to  another,  and  the 
transient  response  in  current,  voltage  and  charge  is 
observed.  For  "thin"  devices  the  transient  response  is 
normal,  in  the  sense  that  the  charge  in  the  device  was 
rearranged  to  accomodate  the  new  DC  conditions,  and  this 
was  done  gradually  according  to  the  appropriate  time  con- 
stant. However,  for  "thick"  devices,  the  DC  bias  change 
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triggered  the  formation  for  Gunn  domains,  and  oscillations 
were  observed.  A travelling  domain  exists  between  gate 
and  drain.  The  analysis  shows  that  for  each  potentially 
unstable  device  there  is  a range  of  DC  conditions  in  which 
the  device  is  unstable,  but  if  the  drain  to  source  voltage 
is  increased  or  decreased  to  bring  the  device  outside  of 
this  region,  the  oscillation  stops. 

1.3.8  Semi-Empirical  Analysis 

As  mentioned  above,  the  numerical  simulations  are 
accurate  and  useful  but  unpractical  for  circuit  design. 

A semi-empirical  approach  to  achieve  a circuit  design 
oriented  simulation  was  developed  by  Willing  and 
Rauscher  [17]  . Their  approach  is  pursued  as  follows: 

1.  Stipulate  an  equivalent  circuit  for  the  FET 
containing  several  fixed  elements  and  some  bias-sensitive 
(non-linear)  elements.  Willing  and  Rauscher  chose  a 
topology  similar  to  the  commonly  used  small  signal 
equivalent  circuit  (see  [17]). 

2.  Measure  the  small-signal  s-parameters  of  the 
device  for  the  full  bias  range  required . 

3.  Calculate  the  fixed  and  bias  sensitive  elements  of 
the  equivalent  circuit  by  curve  fitting  the  measured  and 
calculated  small-signal  s-parameters,  thus  getting  the  full 
quantitative  representation  of  the  device. 
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Large  signal  analysis  is  performed  by  embedding  the 
; above  equivalent  circuit  in  the  complete  network  to  be 

L 

analyzed  and  solving  the  network  time-domain  differential 
equations . 

The  advantage  of  the  above  approach  is  its  efficiency 
and  low  cost  (once  the  equivalent  circuit  is  known) . The 
disadvantages  are:  1)  a long  and  costly  process  of  measure- 
ments to  establish  the  equivalent  circuit  for  each  device; 
and  2)  no  direct  dependence  on  device  geometry  and  other 
properties,  since  this  is  a phenomenological  approach,  and 
is  not  based  on  first  principles. 

1.4  SCOPE  OF  THIS  REPORT 

As  mentioned  in  the  beginning  of  this  chapter  it  is 
devoted  to  a review  of  the  main  FET  analysis  methods 
published  in  the  literature.  The  rest  of  the  report 
includes  the  following  topics: 

Chapter  2 develops  the  basic  analysis  approach 
employed  in  the  new  model  as  suggested  by  Yamaguchi  and 
Kodera  [1] . Chapter  3 includes  the  actual  derivations  for 
the  gate-region  analytic  model.  Emphasized  are  the 
improvements  made  to  make  the  model  "true  AC".  The 
complete  FET  model  (including  internal  and  external 
parasitics)  is  presented  in  chapter  4.  The  application  of 
the  model  is  illustrated  in  chapter  5.  The  data  generated 


A 


by  the  model  are  presented  in  graphical  form,  and  include 
static  properties,  s-parameters,  waveforms,  harmonic  i 

distortion,  etc.  Chapter  6 contains  conclusions  and 
suggestions  for  future  work. 

/ 
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2.  THE  BASIC  APPROXIMATE  ANALYTIC  SOLUTION 

2 . 1 INTRODUCTION 

The  fundamental  ideas  and  derivations  upon  which  the 
present  FET  model  is  based  are  presented  in  this  chapter. 

First  outlined  is  the  basic  approach  of  Yamaguchi  and 

Kodera  [1],  The  idea  of  postulating  a "reasonable"  / 

functional  behaviour  for  the  charge  carrier  concentration  is 

| 

explained,  and  the  derivation  of  the  resulting  parametric 
analytic  solution  is  presented.  The  analytic  determination 
of  the  parameter,  which  is  one  of  the  improvements  reported 
here,  is  derived  in  chapter  3. 

Once  the  parameter  is  known,  the  problem  is  solved 
and  the  electric  field  and  charge  carrier  density  can  be 
calculated.  Thus,  by  calculating  the  current  density  and 
integrating  over  the  proper  electrodes,  the  drain,  source 
and  gate  currents  can  be  calculated.  These  are  the  total 
currents  (conduction  and  displacement) , and  thus  depend  on 
the  terminal  voltages  and  their  time  derivatives.  The 
expressions  for  the  currents  are  derived  in  this  chapter. 

Also,  from  the  displacement  current  expressions,  the  small 
signal  input  and  output  capacitances  are  computed  as 
functions  of  the  DC  bias  voltages. 

2.2  PROBLEM  DEFINITION 

The  device  analysed  in  this  report  has  a structure  as 
presented  in  Figure  1-4,  namely  it  is  a planar  structure, 
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with  the  three  electrodes  deposited  on  top  of  the  epitaxial 
layer.  The  source  and  drain  electrodes  form  an  ohmic 
contact,  while  the  gate  electrode  forms  a Schottky  barrier 
junction.  The  semi-insulating  substrate  is  considered  a 
perfect  insulator.  Of  course,  it  is  possible  to  account 
for  substrate  losses  by  an  external  effective  resistor 
between  the  source  and  drain  contacts. 

For  simplicity,  only  the  active  region  directly  under 
the  gate  is  considered  in  the  analysis.  The  effect  of  the 
two  regions  outside  the  gate  electrode  (source-gate  and 
gate-drain  spacings)  is  considered  in  chapter  4.  Therefore, 
the  geometrical  structure  to  be  analysed  is  as  given  in 
Figure  2-1.  The  source  reference  point  is  (0,0).  At  this 
point  the  potential  is  taken  as  zero.  The  drain  point  is 
at  (lg,  0) , and  the  potential  there  is  assumed  VDg  (a  given 
boundary  condition) . The  gate  extends  from  x=0  to  x=lg  at 
y=a.  The  potential  at  the  gate  is  "(<I>B+VSG),  where  <f>B  is 
the  built-in  potential  of  the  Schottky  barrier  junction  and 
VgG  is  a given  boundary  condition.  The  differential 
equations  to  be  solved  are  (1-47)  to  (1-50) , which  are 
reproduced  here  for  convenience: 


2'l>  = 

- f <ND-n> 

(2-1) 

1"  — 
3t 

V • J 

(2-2) 

= 

-qnV  + qDVn 

(2-3) 

tot 

j + e P 

0 t 

(2-4) 
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The  above  equations  represent  two  coupled  partial 
differential  equations  in  the  variables  \ [>,  n.  The  boundary 
conditions  are: 


-e- 

o 

o 

= 0 

(2-5) 

V (x,  a) 

" “(*B  + VSG) 

(2-6) 

dll)  (x,  0) 

3 y 

3n(x,0)  _ „ 

3 y 

(2-7) 

<|»  (lg,  0) 

ii 

< 

a 

cn 

(2-8) 

(source  conduction 

current)  = (drain 

conduction  current) 

(2-9) 


Condition  2-9  is  strictly  correct  only  for  DC  since 
then  there  is  no  conduction  current  via  the  gate.  For  AC 
drive  the  above  is  incorrect  for  instantaneous  currents 
since,  for  example,  in  transient  response  some  of  the 
source  conduction  current  is  used  to  charge  the  depletion 
region.  In  the  present  analysis  the  assumption  is  made 
that  (2-9)  holds  for  AC  too.  In  chapter  4 a method  is 
presented  to  compensate  for  the  error  introduced  by  using 
Equation  2-9. 

Another  condition  to  be  specified  before  the  equations 
can  be  solved  is  the  function  of  charge  carrier  velocity 
versus  electric  field.  This  function  is  a property  of  the 
material.  Figure  2-2 (a)  presents  the  velocity  field 
relation  for  GaAs.  Typical  values  are:  Ec  = 3.2  KV/cm, 
v = 107  cm/sec,  v = 2.107  cm/sec.  For  simplicity,  and 
in  view  of  the  discussion  in  section  1.3.6  the  function 
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Figure  2.2  Velocity-electric  field  curves 
(a)  GaAs  (b)  functions  assumed 
in  the  present  analysis 
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used  in  the  present  analysis  is  given  in  Figure  2-2 (b) 
(curve  1) . This  is  a piecewise  linear  approximation,  which 
neglects  the  negative  differential  mobility  region.  This 
approximation  is  appropriate  for.  "pentode  like"  operation 
devices  (section  1.3.6),  namely,  for  a < Wfch  (eq.  1-57) 
(for  these  thin  devices  Gunn  domains  do  not  form) . For 
thicker  devices  Gunn  domains  form  near  the  drain  electrode. 
An  approach  to  handle  this  situation  is  outlined  later  in 
this  report.  The  function  represented  by  curve  2 in 
Figure  2-2 (b)  is  a smooth  approximation  with  a continuous 
derivative  at  2Ec.  This  curve  is  used  in  chapter  3. 

The  various  curves  in  Fig.  2-2  represent  the  relation- 
ship between  the  magnitudes  of  velocity  and  electric  field. 
It  is  further  assumed  that  V and  E are  colinear,  thus 
yielding  the  relationship  between  the  cartesian  components: 


E 


VE) 

- irv(E) 

(2-10) 

vy(E) 

Ey 

= / v(E) 

(2-11) 

where. 

E 

-/  E 2 + E 2 

(2-12) 

V x y 

Yamaguchi  and  Kodera  II]  solved  the  above  equations 
numerically  as  mentioned  in  Section  1.3.6.  From  these 
results  they  suggested  a parametric  function  for  the  charge 
carrier  concentration  n(x,y).  Using  this  function  in 
Equation  2-1  yields  an  analytical  solution  to  the  potential 
function  i|»(x,y).  The  details  are  given  in  the  next  section. 
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2.3  THE  APPROXIMATE  CHARGE  CARRIER  DENSITY  FUNCTION 


As  described  in  chapter  1 much  work  has  been  done  in 
the  last  few  years  on  the  numerical  simulation  of  GaAs  FET 
characteristics.  Examination  of  charge  carrier  density 
plots  generated  by  many  workers  (i.e.  Yamaguchi  and 
Kodera  II],  Kennedy  and  O'Brien  [5],  Barnes  [14])  shows 
that  the  Shockley  assumption  of  subdividing  the  device  into 
two  regions  [channel  and  depletion)  is  invalid  for  the  high 
frequency  devices  considered  here.  In  fact,  there  is  an 
obvious  transition  region  in  which  the  charge  carrier 
density  changes  gradually  from  NQ  in  the  channel  to  zero  in 
the  depletion  region.  Based  on  the  above,  Yamaguchi  and 
Kodera  suggested  subdividing  the  device  into  the  three 
regions  shown  in  Figure  2-3: 


1. 

channel  (n=N^) : 

0 

1 A 

< 

dx(x) 

2. 

transition  : 

d1  (x) 

< 

y < d1(x)  + d 

3. 

depletion  (n=0): 

dx  (x) 

+ 

d < y < a 

The  charge  carrier  density  is  approximated  in  the 
transition  region  by  the  function: 

1+cos  5-(y-dn  (x)  ) 

n(x,y)  = ND[l-a(x-y)]  § ± (2-13) 

where. 


d = 6 XD  (see  eq.  1-46) , a 


= 0 for  x < y 
/ 0 for  x > y 
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The  function  d^ (x)  and  the  parameters  a,y  are  determined 
by  using  the  boundary  conditions,  as  described  in  the 
next  section.  Now  by  substituting  the  assumed  n(x,y)  in 
Equation  2-1  one  obtains  a relatively  simple  differential 
equation,  which  can  be  solved  analytically.  There  are,  of 
course,  three  different  equations,  one  for  each  of  the 
regions. 

It  is  important  to  note  that  the  suggested  n(x,y) 
implies  three  operational  modes  for  the  device  as  shown  in 
Figure  2.4: 

Mode  A:  all  three  regions  exist  across  the  entire 
device.  This  mode  occurs  for  low  enough  V„_,  such  that  the 
channel  is  not  pinched  off.  tY>lg,d^(lg)  > 0,  d^(0)  > 0) 

Mode  B:  The  neutral  channel  exists  on  the  source  side, 
but  not  on  the  drain  side  (0<y<lg) . This  occurs  for  large 
enough  VgG  to  pinchoff  the  drain  side,  but  not  enough  to 
pinchoff  the  source  side.  (d1(lg)=0,  d^(0)>  0) 

Mode  C:  The  neutral  channel  does  not  exist  anywhere 
across  the  device  (y<0) . This  occurs  for  large  enough 
Vgr  to  pinchoff  the  channel.  The  device  is  not  cutoff, 
since  there  is  some  charge  transport  via  the  transition 
region  (d^O)  = d1(lg)  = 0). 

There  is,  of  course,  also  operational  mode  D in  which 
the  device  is  cutoff.  The  exact  driving  conditions 
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(V„„,  V__ ) for  each  mode  are  derived  in  the  next  sections. 
Do 

The  properties  of  the  device  in  each  one  of  the  modes  are 
different,  and  must  be  derived  separately. 

It  is  important  to  note  that  the  assumed  n(x,y)  does 
not  allow  for  charge  accumulation  in  the  channel,  as  was 
proven  to  exist  by  several  authors  (see  sections  1.3.4, 
1.3.6).  The  accumulation  considered  here  is  not  related 
to  Gunn  domains  but  to  current  continuity.  Neglecting  this 
accumulation  may  cause  some  problems  in  simulating  thick 
channel  (d^ (lg)  > d)  devices.  A method  to  compensate  for 
this  error  is  suggested  in  chapter  3. 

2.4  THE  SOLUTION  TO  POISSON'S  EQUATION 

The  solution  of  Poisson's  Equation  2-1  with  n(x,y) 
from  the  previous  section  is  discussed  here.  As  a first 
step  assume  that  the  solution  can  be  written  in  the  form 

(2-14) 

(2-15) 


V = + ’I'i 

where, 

v2»0  = 0 

V2„1  - - f |nd  - n(x.y)} 


(2-16) 
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the  boundary  conditions: 


t|/Q(0,0)  = ^(0,0)  = 0 


(x' a)  = 0 


^1  (x,  a)  = ~ (0>B  + VSG) 


3(|/0(x,0)  S^^O) 


3y 


3y 


= 0 


(2-17) 

(2-18) 

(2-19) 

(2-20) 


Also  define  VQ  as: 


vx  = 'J'ji  dg> 0)  VQ  = i^0(lg,0) 


V1  + V0  “ VDS 


(2-21) 

(2-22) 


<|/g  is  the  solution  to  Laplace's  equation  and  can  be 
written  generally  as: 

*0=2  *omsinh(IJ  x)  cos  (2-23) 

m 

For  simplicity  in  this  analysis,  only  the  first  term  of  the 
series  is  used  (see  [1]): 

*0  = .t„!%ia,  3inhl7?lco3(%1  (2-24) 

sinh  2ii 

Since  the  solution  to  depends  on  n(x,y),  it  has  to 
be  treated  separately  for  each  region: 
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1.  neutral  channel  (.n=N^ , excess  charge  = zero) 
(0<y<d1(x)).  Due  to  charge  neutrality  and  boundary 
condition  (2-20)  is  a linear  function  of  x: 

ip1  = f (x)  = CjX  + C2  (2-24) 

To  simplify  the  following  derivations,  an  assumption 
is  made  that  d1  (x)  decreases  gradually  in  the  x direction, 
namely: 

92d,  (x) 

% < < 1 (2-25) 

9 x^ 

Also  along  the  interface  between  the  channel  and  transition 
region  (y=d^(x))  and  should  be  continuous. 

2.  Transition  region  (d^x)  < y < d^x)  + d) 

Using  (2-25)  and  the  continuity  requirements  the  solution: 


qNr 


qN 


D ,d. 


*1  = * TT(*-di(x))  + 2TW  1-cos 


irfy-d-^x)  J 


qN. 


+ -^a(x-Y) 


(l<! 


2 irfy-d^x))  (y-dx (x) ) 


2 


) cos 


| + f (X) 


(2-26) 
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3.  Depletion  region  (d^ (x) +d<y<a) 

Similarly,  the  solution  for  this  region  is:  i 

| 

gND  2 gN 

= — 2^—  {y- (dj^CxJ+d)  } - d {y- (^  (x) +d)  } {l+a(x-y)}+  ^ 

gN_.  o i 9 

+ ~T  d' {L2  - V~  ~T  f {1+  ^2)atx-Y)  + f(x)  (2-27)  ‘ 

IT  IT  f 

The  validity  of  2-26,  2-27  can  be  easily  verified  by 
substitution  in  2-1  and  using  2-25.  Condition  2-25  turns 
out  to  cause  some  problems  in  the  solution.  This  condition 
is  generally  satisfied  everywhere  except  at  x = y 

3d.  (x) 

(see  Figure  2-4  (b)).  At  this  point  — jumps  from  a 
non-zero  value  to  zero;  thus  the  second  derivative  is 
infinite  there.  It  was  found  that  this  discontinuity 
causes  discontinuities  in  the  capacitance,  and  other 

functions  which  depend  on  derivatives  of  d^(x)  as  the 
device  switches  from  mode  B to  mode  C.  Thus,  some  compen- 
sation is  required  to  avoid  this  discontinuity.  It  is 
described  in  chapter  3 . 

Equations  2-24,  2-26,  2-27  include  the  following 
parameters:  C^,  C2,  a,  y and  the  function  d^ (x) . These 
must  be  determined  before  the  solution  is  considered 
complete.  Note  also  that  in  mode  A and  the  region 
0 < x < y of  mode  B,  the  third  term  of  2-26  and  the  fourth 
term  of  2-27  are  zero  since  a=0.  Also,  in  the  region 

. 
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Y < x < lg  of  mode  B and  mode  C,  d^(.x)  = 0, 

The  quantities  and  C2  are  determined  by  using  the 
boundary  condition  at  x=0  and  lg  (2-17,  2-21),  which  yield: 


C2  = 0 (2-29) 

Now  use  condition  2-19  in  eq.  2-27  (with  ct=0)  and  get 
the  function  d^(x): 


(2-30) 

The  pinchoff  voltage,  V , is  defined  as  the  value  of 
V1  which  sets  d^lg)  = 0.  From  Equation  (2-30)  get: 

VP  - |(«-d)2+d(a-d)  + (l  - 2j)d2j  - <VVS(.)  <2-311 

Thus,  for  vi  < vp  “ operational  mode  A 

for  > Vp  - operational  mode  B,  C or  D 

\ 

Also,  in  the  region  y<x<lg  of  mode  B,  ^ at  y=a  is 
equal  to  -(<J>B+VgG),  from  condition  (2-19).  Inserting  this 
into  2-27  and  using  2-31  yields: 
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qNn  A2  A2 

fix)  = -^(-2-  ~ f + -%)a(x-Y)  + Vp 

2tt 


(2-32) 


However,  f(x)  was  already  evaluated  as  X(Equations  2-28, 
2-29) , so  now  a and  y can  be  evaluated  by  equating  the  two 
expressions. 

Define: 

qNr 


o _ D,a  d d . . 


(2-33) 


Then: 


a = 


Y = 


Big 


(2-34) 


(2-35) 


Thus  far,  all  the  parameters  in  the  solution  of  i|^  are 
expressed  in  terms  of  a single  parameter  V^.  The  solution 
to  tf>0  is  expressed  in  terms  of  the  parameter  VQ  = VDg  - 
So,  the  solution  to  |i  = comPlete  once  is 

known . 

The  parameter  is  determined  by  using  condition  2-9, 
namely,  conduction  current  continuity.  To  apply  this 
condition  one  must  derive  the  electric  field  expressions 
at  the  source  and  drain  side.  Then  use  Equation  2-3  to 
derive  the  current  density  expressions.  This  is 
complicated  by  the  nonlinearity  of  the  v-E  relationship 
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(velocity  saturation  - Figure  2-2) . Then  integrate  the 
current  density  expressions  along  the  source  and  drain 
electrodes  and  obtain  the  source  and  drain  current 
expressions,  expressed  in  terms  of  the  parameter  V^. 

Now,  find  such  as  the  two  currents  are  equal. 

The  above  process  is  quite  lengthy,  and  is  generally 
done  in  a numerical  fashion,  namely  the  equation  Ic  = In 
is  solved  by  using  numerical  methods  such  as  Newton-Raphson. 
This  is  the  approach  employed  by  Yamaguchi  and  Kodera  [1] . 
Thus,  even  though  an  approximate  analytic  solution  has  been 
achieved  for  the  FET,  the  determination  of  the  parameter  is 
done  numerically. 

This  situation  has  been  rectified,  and  in  chapter  3 
a method  of  achieving  an  almost  analytic  solution  for 
is  presented.  The  method  is  considered  almost  analytic, 
since  it  has  been  found  that  the  results  are  improved  by 
employing  one  additional  numerical  iteration  after  the 
approximate  analytic  data  is  calculated.  This  new  method 
was  found  to  be  about  one  order  to  magnitude  faster  than 
Yamaguchi  and  Kodera' s numerical  approach.  A typical 
computation  time  for  a given  pair  of  (V„^,V„0)  is  about 

ovj  Do 

0.1  seconds  on  an  IBM  360.  This  fact  makes  the  new  model 
practical  for  circuit  design,  in  which  the  model  has  to  be 
called  typically  tens  of  thousands  of  times. 
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For  the  rest  of  this  chapter  it  is  assumed  that  the 
parameter  has  been  evaluated,  thus  yielding  a complete 
solution  for  the  potential,  and  charge  carrier  concen- 
tration, n.  In  the  next  section  expressions  are  derived  for 
the  displacement,  currents  at  the  device  terminals.  The 
expressions  for  conduction  current  are  derived  in  chapter  3, 
since  they  are  closely  related  to  the  procedure  of 
determining  the  parameter  v^. 

2.5  DISPLACEMENT  CURRENT  EXPRESSIONS 

The  expressions  for  the  electric  potential  tp(x,y), 
which  were  derived  in  the  previous  section,  can  be  used  to 
derive  expressions  for  the  electric  field  E = - 7i(i.  The 
procedure  is  straight  forward  and  appendix  8.1  contains  a 
full  listing  of  the  electric  field  expressions  at  the  drain, 
source  and  gate  for  each  one  of  the  three  operational  modes 
and  in  each  region  (channel,  transition,  depletion) . 

The  above  expressions  can  be  used  to  calculate  the 
displacement  current  at  each  of  the  electrodes  for  all  the 
operational  modes.  The  displacement  current  density  is 
defined: 


(2-36) 


By  taking  the  time  derivative  of  the  expressions  in 
appendix  8.1  ^ can  be  calculated.  Integrating  jQ  over 

I 
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the  proper  electrode  yields  the  displacement  current  via 
this  electrode.  This  procedure  is  demonstrated  here  for 
the  gate  displacement  current  in  operational  mode  A.  The 
expressions  for  the  other  electrodes  and  operational  modes 
are  derived  in  a similar  fashion. 

The  displacement  current  density  at  the  gate  for 
mode  A is  calculated  by  taking  the  time  derivative  of 
Equation  8.1-12: 


3E 


3D~  e-3t  e 


ir/2a 
,tt1 


sinh(^) 


sinh  (2a)  3t  + 


qNr 


2e  ^Ig  +VSG+<^B^ 


(*-  Hi  + !!§£, 

vlg  3t  3 t 


(2-37) 


In  2-37  the  expression  for  d1(x)  (Equation  2-30)  was 

used. 

The  total  displacement  current  at  the  gate  is 
calculated  by  integrating: 


DIS 


■/’ 

x=0 


jD  dx 


(2-38) 


Using  also  Vq  = VDS-V^  (W  = device  width)  get: 


dV  dV  dV. 

IDIS  = ew  [KGS  -ST  + KDS  -gL  + KG 


dt 


dt“J 


(2-39) 
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where : 


KGS 


A^e  (7V  1+VSG+<(,B  “7VSG+4>B) 


KDS  = [cosh(5|2)  -lJ/sinhCjl2) 


(2-40) 


(2-41) 


KG  = 


2_l2  (rp- 
2e  V,  WV1 


+VSG++B  - 3^<V1+VSG+V 


3/2 


3/2 


‘(VSG+*B)  W-KDS 


(2-42) 


Equation  2-39  is  not  the  final  expression  since  it 
includes  the  time  derivative  of  the  parameter  V, . As 
mentioned  above,  the  procedure  for  determining  V,  is 
described  in  chapter  3.  Since  V,  is  an  implicit  function 

°f  VSG  and  VDS: 


dVl 

dt 


3V1  dVS(,  3V,  dVDS 

3VSG  dt  3VDS  dt 


(2-43) 


define; 


3V, 


3V, 


K1  = 


3V 


K2  = 


SG 


3V. 


(2-44) 


DS 


Then,  Equation  2-39  becomes: 


dV  dV 

XDIS  = EW[  (KGS+KG*K1)  -g|£  + (KDS+KG-K2) 


(2-45) 


v 
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The  variables  Kl,  K2  are  calculated  by  taking  the 
time  derivative  of  the  defining  equation  for  V^.  This 
procedure  is  described  in  chapter  3.  Equation  2-45  is  the 

I— 

final  expression  for  the  gate  displacement  current  in 

operational  mode  A.  It  expresses  the  above  current  as  a 

function  of  the  terminal  voltages  (,V__,  V-„)  and  their  time 

oG  DS 

derivatives. 

The  displacement  current  expressions  for  the  other 
electrodes  and  other  operational  modes  are  listed  in 
appendix  8.2. 
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3.  THE  NOVEL  FET  MODEL  FOR  THE  ACTIVE  REGION 
3.1  INTRODUCTION  . 

This  chapter  presents  the  method  used  in  the  new  FET 
model  to  determine  the  parameter  V^,  which  appears  in  the 
solution  for  and  n,  and  therefore,  in  all  relevant 
expressions.  As  mentioned  in  chapter  2,  V1  is  calculated 
by  equating  the  source  and  drain  conduction  current. 
Therefore,  the  exact  expressions  for  these  currents  as  a 
function  of  are  derived.  Due  to  the  complexity  of 
these  equations,  it  is  impossible  to  solve  for 
analytically.  Next,  some  approximations  are  made,  which 
allow  simplification  of  the  conduction  current  expressions, 
and  eventually  enable  an  approximate  analytic  solution 
for  to  be  obtained. 

An  improvement  to  the  approximate  solution  is  achieved 
by  adding  an  additional  numerical  corrective  iteration. 

For  typical  microwave  devices  one  corrective  iteration  is 
sufficient.  As  seen  in  section  2.5  expressions  are  needed 
for  the  derivatives  of  V.  with  respect  to  Van  and  V_„. 

These  expressions  are  derived  using  the  exact  conduction 
current  expressions,  because  this  yields  more  accurate 
results  (derivatives  are  much  more  sensitive  to  errors 
than  the  function  itself) . Also,  since  is  already 
known  it  is  reasonable  to  use  the  exact  expressions. 
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Also,  considered  in  this  chapter  is  a modification  in 
the  equations  caused  by  the  use  of  a smooth  velocity-electric 
field  curve,  and  another  modification  to  allow  for  a small 

L. 

charge  accumulation  at  the  drain  side  of  the  channel. 

3.2  CASE  CLASSIFICATION  INDUCED  BY  VELOCITY-ELECTRIC 
FIELD  CURVE 

The  assumed  charge  carrier  velocity  dependence  on  the 
electric  field  is  as  shown  in  Figure  2-2 (b)  (curve  1). 

This  is  a piecewise  linear  approximation,  which  considers 
a constant  mobility,  yQ,  for  E < Ec  and  a constant 
velocity  for  E > Ec . Since  the  conduction  current 
expressions  are  different  for  constant  mobility  and 
constant  velocity,  there  are  three  different  cases  to  be 
considered: 

case  1:  E < Ec  no  velocity  saturation  at  either 
source  and  drain. 

case  2:  E < E„  at  the  source,  but  E > E at  the  drain, 
c c 

i case  3:  E > E at  both  source  and  drain. 

c 

! The  case  E > E at  the  source,  but  E < E at  the  drain  is 

j c c 

identical  to  case  2,  except  for  the  interchange  of  drain 
and  source.  Therefore,  it  is  not  considered  a separate 
case.  It  is  important  to  note,  that  the  "test-point"  for 
the  source  electric  field  is  at  (x=0,y=0),  and  for  the  drain 
electric  field  is  at  (x=lg,y=0) . As  seen  in  the  equations 
in  appendix  8.1  the  electric  field  generally  increases  with 
y,  such  that  even  if  at  the  "test  point"  E < Ec , the 
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electric  field  may  exceed  Ec  for  a large  enough  y.  This 
is  considered  further  in  the  next  section. 

Since,  there  are  three  different  operational  modes  for 
the  device  (as  described  in  chapter  2) , and  for  each  mode 
there  are  three  different  cases  (as  described  in  the 
present  section) , there  is  a total  of  9 different  cases  in 
each  of  which  the  FET  is  characterized  by  a different  set 
of  equations. 

3.3  EXACT  CONDUCTION  CURRENT  EXPRESSIONS 

To  demonstrate  the  calculation  of  the  conduction 


current,  outlined  below  is  the  derivation  for  operational 
mode  C-case  2,  namely,  the  case  of  velocity  saturation  at 
the  drain,  but  not  at  the  source.  From  Figure  2-4 (c)  it 
is  obvious  that  the  contribution  to  the  conduction  current 
in  mode  C comes  only  from  the  transition  region. 


3.3.1  Source  Conduction  Current 

The  electric  field  on  the  source  side  (transition 
region)  is  given  by  Equations  8.1-19  and  8.1-20;  which  are 
duplicated  here  for  convenience: 

V1  2 

E = - [y— (1+Ky1  {2cos6-0z})+  K V cos(d.0)]  (3-1) 

x X x y2  o x 


qNDd 

Ey  = 2ir£  - [0  - sin©  + Kc  (0  + sin©)] 


(3-2) 
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where , 


9 = 

IZ 

d 

(3-3) 

;yi  = 

qND  .d  2 

Teg  7T 

(3-4) 

tr/2a 

(3-5) 

y2  = 

sinh(^2.) 

Kc  " 

-Vp/6 

(3-6) 

di  - 

d/2a 

(3-7) 

In  case  2,  'which  is  considered  here  E(o)  < Ec,  namely, 

V. 

E(o)  = |Ex(o)  |=  I±.(l+2Kyl)  + Ky2-VQ  < Ec  (3-8) 

However,  as  y increases  Ey  increases,  and  for  large 
enough  y,  E(o,y)  may  be  larger  than  Ec . The  critical 
9=9  at  which  E(o,y)  = Ec  is  defined  by  the  equation: 

2^2  2 

Ex<9)  + Ey  (9)  = Ec  (3-9) 

Equation  3-9  can  be  solved  for  9 numerically  (Newton- 
Raphson)  if  is  known.  When  employing  a numerical  solu- 
tion for  in  an  iterative  fashion  a value  for  is 
assumed  to  start  with  and  this  value  is  corrected  each 
iteration.  Thus  in  equation  (3-8)  the  current  value  of 
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V1  is  used  to  solve  for  For  the  analytic  solution  F is 
estimated,  as  described  later  in  this  chapter. 

In  the  range  0 < 0 < (T  the  velocity  is  not  saturated 
and  the  contribution  to  the  drift  conduction  current  is 
calculated  as  follows.  The  drift  current  density  is  given 
by  (Equation  2-3) . 

jx  " -9nVx  “ (3"10) 

The  electron  density  n at  the  transition  regi  i s given  by 
Equation  2-13.  Using  the  appropriate  expressi . ns  for  a,  y 
and  d1(x)=0  the  expression  for  n at  the  source  is: 


n 


source 


ND[1“Kc] 


1+COS0 

2 


(3-11) 


Thus,  the  contribution  to  the  drift  current  is 


6d 


"/  V*  - “ / 


jxae 


(3-12) 


Inserting  3-1  and  3-11  in  3-12  yield  the  solution: 

IL  = CDlCP  (1-Kc)  (?+sin(T)  Vx  + 

+CDOC  (1-Kc)  [sin(d1'IT)+d3sin(d2l’)+d5sin(d49')  ] »V 


where. 


CDlCP 


qNdWdyc 
2 irlg 


(3-13) 


(3-14) 
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CDOC 

qNDM0w 

2sinh(g2.) 

(3-15) 

d2 

= 1 + d1 

(3-16) 

d3 

= d1/(2d2) 

(3-17) 

d4 

" 1 “ d! 

(3-18) 

d5 

= dj/(2d4) 

(3-19) 

In  the  region  0 < 0 < tt  the  electron  velocity 

is 

saturated; 

thus  (see  Equation  2-10) : 

‘ -^nVx  - -C‘nva  r,  X , 

vE  +E  Z 
x y 

(3-20) 

Inserting  n from  3-11  into  3-20  and  integrating  over 
the  range  0 < 0 < ir  yields  the  drift  current  contribution 
from  the  upper  region: 


IT 

I„  = CD2CP • (1-K  )•  j (l+cos0)d0_  (3-21) 

U C / 2 

0 \ll+  (Ey/Ex)  z 

where 

qN  v Wd 

CD2CP  = — (3-22) 

In  Equation  3-21  the  expressions  for  Ex,  Ey  are  taken 
from  equations  3-1  and  3-2.  The  integration  in 
equation  3-21  can  be  performed  numerically  when  is 
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known.  In  the  analytic  solution  this  integral  is  approxi- 
mated by  an  analytic  function,  as  described  later  in  this 
chapter.  Figure  3-1  illustrates  the  above  mentioned 
partitioning  of  the  transition  region  at  the  source  side. 

The  diffusion  current  density  is  given  by  (see 
Equation  2-3) 


= qDVn 


(3-23) 


DIFX  = qD(Vn)x 


(3-24) 


The  diffusion  coefficient,  D,  is  in  general  a function 

of  the  electric  field.  To  simplify  calculations,  and 

since  the  electric  field  at  the  source  for  the  case 

considered  is  quite  low  (less  than  3.2  KV/cm) , an  average 

2 

low-field  constant  value  is  assumed;  DT  = 250cm  /sec.  The 

Li 

expression  for  (Vn)x  = is  derived  from  Equation  2-13 
with  the  appropriate  expressions  for  a,  y and  d1(x)  = 0: 


IS  i - - w(1+cose) 


(3-25) 


source 


Using  3-25  in  3-24  and  integrating  yields: 


/*d  * 

= W J Wx  = F/  W/e  = CDIFF.Vl  (3-2  6) 
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where. 


CDIFF 


qDLNDwd 

201g 


(3-27) 


Thus,  the  total  source  conduction  current  is  given  by 
the  sum  of  equations  3-13,  3-21,  3-26 


IT  + I„  + 

L U DIF 


(3-28) 


3.3.2  Drain  Conduction  Current 

The  electric  field  on  the  drain  side  (transition 
region)  is  given  by  Equations  8.1-14  and  8.1-15,  which  are 
duplicated  here  for  convenience: 

V1  2 

Ex  = -[j^(l+Ky1{2cos9-e'‘})+  Kh*VQ  cos(d19)]  (3-29) 


qNDd 


V1'VP 


Ey  = 2i ~ t0_sin0+  — g (0+sin0)]+  sin  (d^ ) (3-30) 


where  0,  Ky1  are  defined  by  Equations  3-3  and  3-4  and 


Kh  = coth  (^) 


(3-31) 


In  the  case  considered  here  (E.  . > E ) the  electron 

drain  c 

velocity  is  saturated,  thus  the  drift  current  density  is: 


-qnVx  = -qnVg 


/e  2+E  2 
' x y 


(3-32) 
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The  expression  for  the  electron  density  is  derived 


from  Equation  2-13  with  the  proper  a,  y and  d1(x)=0: 

V -V 

_ „ fn  1 P,  1+COS0 
ndrain  " ND[1 T“]  


(3-33) 


Inserting  3-33  into  3-32  and  integrating  over  the  drain 
(y=0  to  d)  gives: 

IT 


DRIFT 


= CD2CP*  (1-  V--g-P)  J (1+CQS9) 


de 


o Vl+ (Ey/Ex) “ 

CP2CP  is  given  by  equation  (3-22)  and  the  expressions  for 
Ex,  Ey  to  be  used  in  3-34  are  given  in  Equations  3-29,3-30. 

The  diffusion  current  density  is  given  by 
(see  equation  2-3) 


(3-34) 


; _ 3n 

JDIFx  ~ qD  3x 


(3-35) 


For  simplicity,  the  diffusion  coefficient  on  the  drain 
side  is  assumed  a constant  equal  to  a high-field  average 
Djj  = 50cm2/sec.  The  expression  for  is  derived  by  the 
same  method  as  in  the  source  side.  The  result  was  found 
to  be  the  same  as  for  the  source  (3-25) . 


Integrating  3-35  yields  the  solutions 


IDIF  = CDIFFH.Vl 


where , 


CDIFFH  = 


qDHWNDd 

23lg 


(3-36) 


(3-37) 
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So,  the  total  drain  conduction  current  is  given  by  the 
sum  of  Equations  3-34,  3-36: 


Icond  IDRIFT  + IDIF 


(3-38) 


3.3.3  Numerical  Solution  For  V. 


To  solve  for  the  parameter  define  the  function 


f(Vl>  " Icon_  - 1 


con . 


(3-39) 


-s  d 

The  right  value  for  sets  f(V1)=0,  namely,  it 
equalizes  the  source  and  drain  conduction  currents.  From 
the  expression  derived  above  it  is  obvious  that  ffV^  is  a 
quite  complex  function  of  V^,  so  that  an  analytic  solution 
to  f(V^)=0  is  unlikely  to  be  found.  However,  a numerical 
solution  is  straight  forward  using  the  Newton-Raphson 
approach:  choose  an  initial  guess  for  and  proceed 
iteratively  using  the  iteration  equation: 

f (vx) 

df/dV. 


V 


1i+l 


= V 


Xi 


(3-40) 


This  approach  was  adopted  by  Yamaguchi  and  Kodera  [1] . 
There  is  another  approach  which  can  be  employed,  especially 
if  the  initial  guess  is  quite  close  to  the  final  solution. 
The  approach  is  to  substitute  for  in  the  equation 
f(V^)=0  the  initial  guess  only  in  the  expressions,  which 
are  complicated  functions  of  (i.e.  the  integrals) . 
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In  the  other  expressions  of  the  above  equation  keep  the 
unknown  V^.  The  equation  thus  derived  is  a simple  equation 
for  V1  (either  linear  or  quadratic),  and  is  readily  solved. 
This  process  can  be  repeated  iteratively  by  each  time 
substituting  for  in  the  exact  equation  f(V^)=0  the  last 
value  computed  and  solving  the  resulting  simple  equation. 
This  method  is  not  guaranteed  to  converge  generally. 
However,  if  it  converges  it  can  save  computation  time, 
since  it  does  not  employ  derivatives.  In  the  computer 
model  presented  in  this  work,  the  above  method  was  employed 
successfully  for  the  corrective  iteration,  as  discussed 
later  in  this  chapter. 

The  iteration  equation  in  the  above  approach  for  the 
operational  mode  C-case  2,  which  was  employed  in  the  model, 
is: 


+ 


(1-Kc) [CD2CP  (DINT-SINT)-CDOC*KDO'Vds] 

(1-Kc)  [CD1CP  ('e+sine’)-CDOC*KDO]+CDIFF-CDIFFH+^^DlNT- 

(3-41) 

DINT  is  the  drain  integral  (Equation  3-34) , SINT  is  the 
source  integral  (Equation  3-21) . Both  of  them  are 
calculated  numerically  using  . 

Kc,  CD2CP , CDOC , CD1CP , CDIFF,  CDIFFH  are  defined  in 
Equations  3-6,3-22,3-15,3-14,3-27,3-37,  respectively,  and 


i 
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KDO  = sin  (d^Q') +d3sin  (d2'S') +d^sin  (d^)  (3-42) 

Equation  3-41  employs  computation  of  the  source  and 
drain  integrals  by  using  the  previous  iteration  value  for 
V^.  Also  an  averaging  process  is  included,  namely,  taking 
the  average  of  two  consecutive  iterations.  This  was  found 
to  improve  convergence  and  avoid  possible  oscillations. 

The  integrations  were  performed  numerically  by 
partitioning  the  integration  region  into  20  parts  and 
using  the  simple  rule: 

0 20 

I f (0)  d0  a f(0i)A0i  (3-43) 

o i=l 

Twenty  partitions  were  found  to  give  adequate 
accuracy.  For  testing  purposes  100  partitions  were 
employed,  and  the  result  thus  obtained  was  within  1-2%  of 
the  result  obtained  with  20  partitions. 

3.4  ANALYTIC  EXPRESSIONS  FOR  CONDUCTION  CURRENT  AND  V1 

In  this  section  some  approximations  are  made  which 
enable  the  source  and  drain  conduction  currents  to  be 
expressed  as  analytic  functions  of  V^,  and  thus  permit 
Vx  to  be  solved  analytically.  The  approach  is 
demonstrated  for  mode  C - case  2,  as  in  the  previous 
section. 
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3.4.1  Source  Conduction  Current 

This  subsection  follows  the  same  steps  as  3.3.1,  except 
that  approximations  are  introduced  when  necessary  to  achieve 
an  analytic  solution.  First  an  estimate  is  derived  for  the 
critical  angle  0.  The  exact  defining  equation  is  3-9.  For 
Ex=0  3-9  becomes  Ey=Ec  or: 

e^sine^^  + Kc(91+sin01)  = S1 

where:  (3-44) 

2ttE  e 

S = c 

1 qNDd 

Equation  3-44  is  solved  for  0^  numerically  by 
Newton-Raphson.  This  equation  is  very  suitable  for  this 
method,  and 
becomes  Ey  = 

©2  - sin02  + Kc(02+sin02)  = s,/^2~  (3-45) 

This  equation  is  solved  for  02  as  above.  The  estimate 
for  the  actual  critical  angle  is: 

? = (01+92)/2  (3-46) 

Since  Equation  3-13  for  the  lower  region  contribution 
is  analytic  in  V^,  it  is  kept  unchanged.  The  contribution 
of  the  upper  part  is  given  by  Equation  3-21,  which  is  a 
complicated  function  of  Therefore  it  must  be 


the  convergence  is  very  fast.  For  Ex=E^  3-9 


Ec  /Jl~ or : 
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estimated  by  an  analytic  expression.  To  achieve  this  *oal, 
consider  the  ratio  of  the  total  drift  current  to  the  lower 
part  contribution  (Equation  3-13)  with  H from  equation 
(3-46).  The  following  assumptions  are  made: 

1.  In  the  E expression  (3-1)  the  V term  is 

X o 

negligible.  For  a typical  device  this  introduces 
negligible  error  since  Ky2  is  very  small  (division  by  a 
sinh  function) . 

2.  Consider  the  above  ratio  for  E =E 

x c 


vi 

(or  ^ * Ec»  since  VQ  term  is  neglected) . In  this  case  the 
ratio  is  a general  function  almost  independent  of  device 
geometry.  The  ratio  thus  achieved  is  used  as  an  estimate 
in  general,  and  not  just  for  E^=Ec. 

Under  the  above  assumption,  the  ratio  function  becomes 

TT 


R 


/■ 


(l+cos0 ) d0 


o Jl+ (Ey/Ex) 
o’  + sin(T 


(3-47) 


where: 


E 2 0-sin0+K  (0+sin0) 

(»*)  = [A  = ] (3-48) 

x 1+Kyl (2cos0-e“) 

A * 1/S1  (3-49) 

In  general  R depends  on  the  parameters  Kc,  A,  Ky1> 
re,  it  was  calculated  numerically  versus  the  above 
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parameters.  The  results  are  plotted  in  appendix  8.3.  The 

parameter  Kyl  is  normally  very  small  for  GaAs  (0.02-0.03). 

and  the  ratio  function  was  found  to  be  practically 

independent  on  it  in  the  following  range  of  the  other  two 

parameters:  0 < Kc  < 1 in  the  mode  C (see  Equation  3-11, 

for  K =0  n =Nn,  for  K =1  n=0) , A depends  on  the  doping 
c max  u c 

17  —3 

level  (for  example  A=5.76  for  NQ=10  cm  , A=1.8  for 
ND=101^cm  ^) . Thus,  in  appendix  8.3  the  ratio  function  is 
plotted  versus  Kc  with  NQ  (or  A)  as  a parameter. 

The  above  functions  were  approximated  by  a fifth  order 
polynomial.  This  was  achieved  by  using  curve  fitting, 
yielding : 

R (Kc ) +RRO+RR1 • Kc+RR2 • K^+RR3 • K^+RR4 • K*+RR5 • (3-50) 

The  six  coefficients  RRO,  RRl...are  functions  of  Np.  This 
dependence  is  given  in  appendix  8.3.  Thus,  the  estimate 
for  the  total  source  drift  current  is 


IDRIFT 

SOURCE 


R(Kc) -I 


L 


(3-51) 


where  IL  is  given  in  equation  3-13. 

The  diffusion  current  is  given  in  Equation  3-26,  and 
since  it  is  a simple  expression,  it  may  be  kept  unchanged. 
However,  for  consistency  it  was  estimated  in  the  actual 
calculation.  In  the  course  of  similar  approximations 
performed  for  modes  A and  B it  was  necessary  to  estimate 
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the  diffusion  current  in  order  to  achieve  an  analytic 
solution  for  V^.  It  was  found  that  a reasonable  estimate 
was  7%  of  the  term  of  the  drift  current.  In  mode  C this 
assumption  is  not  necessary,  but  for  consistency  and  to 
avoid  discontinuity  it  was  employed  there  too.  Since  in 
FETs  the  diffusion  current  is  generally  much  smaller  than 
the  drift  current,  the  above  assumption  is  not  expected  to 
introduce  any  appreciable  error.  Also,  as  mentioned  above, 
an  additional  corrective  numerical  iteration  is  employed 
using  the  exact  current  expression.  Thus,  finally,  the 
estimated  source  conduction  current  is  given  by: 

Icon  ~ t1*07  CD1CP*  (1-Kc)  (ef+sin?)  *VX  + 
s 

+ CDOC* (1-KC) *KD0*Vo] *R(Kc)  (3-52) 

3.4.2  Drain  Conduction  Current 

The  exact  drift  drain  current  is  given  by 
Equation  3-34.  To  achieve  an  analytic  expression  the 
integral  in  the  above  expression  must  be  estimated  subject 
to  simplifying  assumptions : 

1.  In  the  Ex  expression  (Equation  3-29)  the  term 
is  neglected.  This  is  a reasonable  assumption  in  the  case 
considered  (no  velocity  saturation  at  source,  and  velocity 
saturation  at  drain)  since  is  typically  a small  number 
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(»0.5V)  so  v0=vDs"vi  usually  larger  than  V^.  Also  the 
VQ  coefficient  is  usually  much  larger  than  the  V1 
coefficient,  since  ^ >>j^. 

2.  coth(-^i2)  « 1.  This  holds  for  practical  micro- 
wave  devices. 

Under  the  above  assumptions: 


( 


1Z) 

Ex' 


0-sin0+f3 (6+sin0)  2 

l— +tan(l28)]  (3-53) 


where: 


it2eVo 

gNDda 


coth(ii2) 


(3-54) 


(3-55) 


l3  = 


V1-VP 


Consider  the  integral: 


(3-56) 


I 


(1+cosO ) d0 
J 1+ (Ey/Ex) 2 


(3-57) 


The  expression  from  3-53  is  used  in  3-57.  This 
integral  is  a function  of  the  three  parameters  l 2, 

Note  that  for  ^-*-0,  I-*-0.  Also  for  (|^)-«- tan  (f20 ) 

and  the  integral  can  be  solved  analytically.  The  result  is: 


U1(i2)=  I = j-sin(^2ir)+  2 (1  *-£  y sin[(l+i2)ir] 

/ 2 2 
*-l+oo 


+ 2(l-£2)  sint  d--t2)1T^ 


(3-58) 


The  integral  function  is  investigated  in  appendix  8.4.  It 
is  calculated  numerically  (100  partitions)  in  the  range: 


0 < ^i  < 5 (with  limiting  value  for  °»  from  3-58) 


£2  = 0.15,  0.4 


0 < £3  < 1 


An  attempt  was  made  to  approximate  the  function  by  an 
exponential,  as  described  in  appendix  8.4,  and  the  best 
fit  to  the  data  was  achieved  by  the  function: 


0.3835 


-2.414  £. 


I « p1[i-u2  e 


(3-59) 


is  given  in  Equation  3-58  and: 
U2  = 1 + 4.798£3 


(3-60) 


The  above  approximation  was  found  to  be  reasonably 
accurate  for  £3  > 1.  As  ^-*-0  the  approximation  fails. 
This  is  not  a serious  problem  since  as  mentioned  above, 
for  the  case  consider  here  Vq  is  large  enough,  so  that 
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I 

I 

practically  £ ^ is  larger  than  1.  If  a situation  occurs  in 
which  £^<1,  the  solution  achieved  using  the  following 
equations  is  inaccurate,  and  this  may  require  more  than 
one  corrective  iteration.  In  the  course  of  using  the  model 
to  analyze  typical  high-frequency  devices  the  largest 
number  of  corrective  iterations  which  had  to  be  used  was  2. 
Thus,  using  the  approximation  of  3-59  in  3-34  get: 

0.3835 

-2.414£. 

Icon  = CD2CP*u1(l-£3) [l-(l+4.798£3)e  1 ] 

d 

(3-61) 

where  u1  is  defined  in  3-58  and  £1,  £2,  £ 3 are  defined  in 
3-54,55,56. 

As  done  in  the  approximation  for  the  source  current, 
here  too  the  diffusion  current  is  estimated.  Since 
Dh=0.2Dl,  the  diffusion  current  on  the  drain  side  is  only 
20%  of  that  on  the  source  side.  On  the  source  side  the 
diffusion  current  was  estimated  at  7%  of  the  drift  current. 
Thus,  the  diffusion  current  on  the  drain  side  is  slightly 
over  1%  of  the  drift  current.  So,  for  the  purpose  of  the 
analytic  approximation  the  diffusion  current  at  the  drain 
side  is  negligible.  Therefore,  the  expression  in 
Equation  3-61  is  an  approximation  of  the  total  drain 
conduction  current. 
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3.4.3  Analytic  Solution  For  V1 

The  analytic  expressions  for  the  source  and  drain 
conduction  current  which  were  derived  above  can  be  used  to 
get  an  analytic  solution  for  the  parameter  . 

Define  the  equation: 


(3-62) 


where  I is  expressed  in  Equation  3-52  and  I is 

con„  r n con. 

• s d 

expressed  in  Equation  3-61.  The  I expression  is  a 

s 

simple  linear  function  of  (set  v0=vDs-vi^  * However,  the 

Icon  expression  contains  an  exponential  dependence  on 
d 

(through  Z-^)  and  a quadratic  dependence  on  (through  Z^\  . 
Neglecting  temporarily  the  exponential  dependence,  equation 
(3-62)  is  basically  a simple  quadratic  equation  in  V^, 
whose  solution  is  straight  forward.  To  solve  the 
exponential  dependence  problem,  in  Equation  3-62  substitute 
VDS  instead  of  VQ  in  the  expression  for  This  is 

justified,  since  a)  VQ  for  the  case  considered  here  is 
usually  larger  than  (as  discussed  in  the  previous 
section)  so  VQ  « VD£  ^Vo=VDS~Vl^'  an<^  k)  the  ^rain  integral 
(see  appendix  8.4)  is  a slowly  varying  function  of  for 


large  Z ^ ( Z ^ > 2). 
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Under  the  above  assumptions  Equation  3-62  can  be 
written  in  the  form: 


av£-bv1  + c = o 


where: 

A 


4 . 798f 
6 


B = l-f+  4'g98£-(6+2Vp)  + P|-— 


C = 


f = 


(6+Vp)  <l-f+  l^Hf-vp)-  ^ vDS 


0.3835 


-2.414f. 
e 1 


(3-63) 

(3-64) 

(3-65) 

(3-66) 

(3-67) 


u2eV, 


ll 


qNr 


qUDWvsdUi 

2ug 


(3-68) 

(3-69) 


1 . 07qN  Wdu 

Dl=  j TT (l-KJRtKJ  (9+sin0) 

g 

qNDli0w 


Do= 


2sinh  (j—) 


(l-Kc).R(Kc).KD° 


The  solution  to  Equation  3-63  is: 


vi- 


/I 

i-  B -4»A«( 


2A 


(3-70) 


(3-71) 


(3-72) 


-T 


■? 
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In  the  computer  model  Equation  3-72  is  used  to 
calculate  a reasonable  estimate  for  the  parameter  V^,  which 
is  usually  very  close  to  the  actual  solution.  However,  due 
to  the  approximations  involved  in  the  process  of  achieving 
the  analytic  solution  some  error  is  expected.  Therefore, 
as  a second  step  in  the  model  a better  estimate  for  is 
computed  by  using  the  iterative  approach  with  the  exact 
equations,  as  described  in  section  3.3.3.  This  part  of  the 
model  is  numerical,  however,  due  to  the  good  analytic 
estimate  given  by  3-72  (zero-order-approximation) , the 
numerical  process  is  very  short  - typically  one  iteration. 
Thus,  this  new  model  is  very  fast  (0.05-0.1  sec  for  a given 
pair  of  voltages  VgQ,  VQg  on  IBM  360)  and  practical  for 
circuit  design. 

Once  the  solution  for  is  obtained  the  modeling 
problem  of  the  active  region  is  practically  solved.  The 
conduction  current  is  computed  using  the  exact  source 
conduction  current  expression  (Equation  3-28) . The 
displacement  current  is  calculated  using  the  equations  in 
appendix  8.2  with  expressions  for  derivatives  of  as 
derived  in  the  next  section.  The  expressions  for  transit 
time,  charging  resistor,  transconductance  and  drain 
conductance  are  derived  in  the  following  sections. 
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3.5  EXPRESSIONS  FOR  DERIVATIVES 


As  discussed  in  section  2.5,  the  calculation  of  the 
displacement  current  involves  the  derivatives 


, K2= 


. These  derivatives  can  be  evaluated 


by  using  the  defining  equation  for  V^,  namely 

f(V.)=I  -I  „ =0,  Equation  3-39.  The  I expression 

1 con  con.  con  c 

S a S 

is  given  in  Equation  3-28,  and  I is  given  in 

con^ 

Equation  3-38.  From  an  examination  of  the  =*bove  expressions 
one  concludes  that  the  function  f(V^)  is  in  fact  a function 
of  the  following  variables:  V^  VSQ,  VDg,  ?.  Now,  take 
the  time  derivative  of  the  equation  f(V^)=0  and  get: 


dVl  dvcr  dV„ 

ACVl*^i  + ACSG*  + ACDS  + ACTET  = 0 


(3-73) 


where : 


ACV1=(1-K  ) [CD1CP*  ( Q+sinQ^) -CDOC  • KDO+CD2CP  • 3 (SINT)  ] + 

C o 

V — V 

+ CD  IFF  - CDIFFH  + CD2CP  *DINT/8-CD2CP  • ( 1- 

P oV^ 

(3-74) 

ACSG= [CD2CP (1-KC)  — -f"- - - CD2CP-SINT  - CDOC*KDO*V  - 
C 3KC  ° 


- CD1CP ( 0+sin0 ) V^]  / 8 + CD2CP* DINT/6  - CD2CP • 


V -V 

M 1 3 (DINT) 

» 3VSG 


(3-75) 
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ACDS=(1-K_)  [CD0C*KD0+CD2CP*  3 

c 3vds 

. CP2CP  (1- 

e 3vds 


(3-76) 


ACTET=  (1-K_)  [CDlCP(l+cos?)V.+CDOC*Vrt-  11M21  + 

x ° de 

+ CD 2 CP  3 (SINT)  ] 


(3-77) 


SINT  is  the  source  integral  (Equation  3-21)  and  DINT 
is  the  drain  integral  (Equation  3-34) . These  integrals 
and  their  derivatives  are  calculated  numerically 
(20  partitions) . In  the  computer  model  SINT  is  computed  in 
routine  UPADER  and  DINT  is  computed  in  routine  DRACUR. 

Some  of  the  derivatives  in  the  above  expressions  can  be 
computed  analytically: 


3 (SINT)  1+cos? 


/l+[Ey(0) /Ex(9)  )2 


(3-78) 


The  expressions  for  Ex,  Ey  to  be  used  in  3-78  are  given  in 
3-1,  3-2. 


d (KDQ)  = d1  [cos  (d.  o’)  +0. 5*cos  (d,”0)  +0. 5*cos  (d  .0”)  ] 

d? 

(3-79) 

To  get  rid  of  the  dependence  on  0"  in  Equation  3-73  it 
is  necessary  to  derive  another  equation.  This  is  achieved 


i 


Equations  3-6,3-22,3-15,3-14,3-27,3-37,  respectively,  and 


by  using  the  defining  equation  for  (5*  (Equation  3-9)  . Take 
the  time  derivative  of  equation  3-9  and  get: 


dv  -i  dv  dv  da- 

AIVl • ■—  + AISG*  + AIDS  + AITET  • = 0 


(3-80) 


where : ^ ~ 

1+Ky. (2cos?-(T  ) 

AIV1  = [ - Z • cos  (d^l))  ] • 


V,  {1+Ky,  (2cos0'-0/)  } 

[— p— + Z*VQ*cos(d10)  ] 


(3-81) 


E 2 ~ *** 

AISG  = (^)  ^ + Qin9  [?-sin0+Kc* (?+sin?) ] (3-82) 

p 

V1{l+Ky1(2cos?-?2)  } 

AIDS  = Z -cos  (d^*©)  [— Yg + Z *VQ  *cos  (dj?)  ] 

(3-83) 

E 2 

AITET  = (^)  • {1-cos 0'+KC*  (l+cose'n*  {?-sinlJ+KC  • (IT+sineO  }“ 
S1 


2Ky  -V, 

-{ — (?'+sin?')+Z.d:L.V1.sin(d1e)  } • 


V,  (1+Ky,  (2cos0'-'e2)  } 

{ — p- + Z.Vo.cos(d10)  } 


(3-84) 


2a»sinh(^p) 


(3-85) 


Combining  Equations  3-73,  3-80  eliminates  -rr  and 

dVl  dVSG  dVDS 

thus  expresses  in  terms  of  and  : 
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i- 


dV^ 

dt 


= K1 


♦ K2  "» 


dt 


dt 


(3-86) 


where : 

. AISG*ACTET-ACSG*AITET 
ACV 1 • AI TET - AI V 1 • ACTET 


(3-87) 


_ _ AIDS • ACTET-ACDS  *AITET 
ACVl'AITET-AIVl*  ACTET 

The  above  two  equations  are  the  desired  expressions 

3V.  3V. 

for  the  derivatives  Kl  = — , K2  = — . 

3VSG  3vds 

3.6  CALCULATION  OF  g AND  g. 

^m  3d 


For  the  purpose  of  small  signal  characterization  of 
the  FET  (such  as  s-parameter  calculation)  it  is  necessary 
to  compute  the  transconductance,  gm,  and  drain  conductance. 


Jd" 


^m 


3 


3V 


SG 


(3-89) 


*d  = 


3(1  ) 

con 


3V, 


DS 


(3-90) 


The  computation  is  done  directly  from  the  definition, 
namely,  taking  the  proper  derivatives  of  the  conduction 
current.  In  the  computer  model  the  drain  conduction 
current  expression  is  used  (Equation  3-38)  . The  derivation 
is  straightforward  and  the  expressions  are: 
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DINT  is  the  drain  integral  (Equation  3-34) • The 
integral  and  its  derivatives  are  calculated  numerically  as 
described  in  the  previous  section.  K1  and  K2  are  computed 
from  equations  3-87,  3-88. 

3.7  MISCELLANEOUS  TOPICS 

This  section  describes  some  of  the  problems 
encountered  in  the  development  of  the  new  model  and  their 
solutions.  Most  of  these  problems  are  due  to  some  imperfec- 
tions of  the  analytic  solution  and  the  basic  assumptions. 
Therefore,  generally  a complete  solution  was  not  possible, 
and  only  first  order  corrections  were  employed. 

3d  (0) 

3.7.1  Correction  of  

3VSG 

As  discussed  in  section  2.4,  the  derivative  of  the 
function  d1 (x)  , which  defines  the  boundary  of  the  neutral 
channel  and  is  given  by  equation  (2-30) , is  continuous 
except  at  the  point  x=y,  where  it  jumps  from  some 


I 
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negative  value  to  zero  (see  Figure  2.4(b)).  Thus  the  second 
derivative  is  infinite  at  this  point,  which  violates 
assumption  2-25,  and  in  fact,  is  equivalent  to  a surface 
charge  located  at  x=y . 

The  practical  effect  of  the  above  is  to  induce  a jump 

in  the  capacitances  (or  equivalently  in  the  displacement 

current  coefficients) . The  jump  occurs  as  the  device  goes 

from  operational  mode  B to  operational  mode  C (as  V_„  is 

increased) . The  jump  is  on  the  order  of  a few  percent. 

Since  the  problem  arises  due  to  a basic  assumption, 

it  is  impossible  to  cure  it  completely  without  changing  the 

entire  solution.  However,  since  practically  the  jump  is 

not  too  large,  it  is  possible  to  improve  the  situation  by 

a first  order  correction  as  follows: 

The  main  effect  of  the  function  d^^  (x)  on  the 

3d. (o) 

capacitance  is  via  the  derivative  — ^ . The  expression 

aVSG 

for  this  function  is  derived  from  Equation  2-30: 


3d.  (0) 


IV 


SG 


qNT 


(a_l 


-dL(0)) 


Theory 


(3-93) 


An  attempt  was  made  to  correct  the  problem  by  adding  a 
constant  to  the  above  expression 


V 

I 
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d-L  (0) 


SG 


practice 


qND(a-|  - dL  (0)) 


- CORATl  (3-94) 


CORAT1  is  a constant,  which  is  set  such  as  to  equate  the 
input  capacitances  (common  source)  at  the  boundary  between 
modes  B and  C calculated  from  the  mode  B equations  and 
mode  C equations.  This  approach  was  proven  successful  for 
practical  devices,  and  it  prevents  the  above  mentioned  jump. 
3.7.2  Charge  Accumulation 

As  discussed  in  section  1.3.4,  numerical  simulations 
of  the  FET,  such  as  Kennedy  and  O'Brien's  [5],  show  that 
there  is  some  charge  accumulation  near  the  drain.  This 
happens  even  for  silicon  devices  and,  thus,  is  not  due  to 
Gunn  domain  formation.  The  accumulation  presented  here  is 
based  on  the  assumption  of  a neutral  channel,  and  no 
allowance  is  made  for  any  charge  accumulation.  For  devices 
which  are  not  too  thick,  the  neutral  channel  is  small  and 
current  continuity  is  achieved  mainly  by  rotation  of  the 
velocity  vector  in  the  transition  region.  However,  for 
thick  devices  the  channel  is  thick,  and  most  of  the 
current  is  carried  through  it.  In  this  case  current 
continuity  cannot  be  achieved  only  by  means  of  velocity 
vector  rotation,  since  in  the  channel  the  velocity  vector 
is  almost  entirely  in  the  X direction  even  for  low  fields. 
Thus,  some  charge  accumulation  is  expected. 


f 
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The  effect  of  neglecting  this  charge  accumulation  in 
the  model  is  to  introduce  a false  negative  slope  region  in 
the  curve  beyond  the  "knee"  for  VgG  around  zero, 

as  demonstrated  in  Figure  3.2.  As  V__  is  increased,  the 

bb 

channel  decreases,  and  current  continuity  can  be  achieved 
by  rotation  of  the  velocity  vector,  so  the  negative  slope 
disappears. 

The  problem  exists,  of  course,  only  for  operational 
mode  A,  since  in  the  other  modes  there  is  no  neutral  channel 
at  the  drain  side.  To  solve  the  problem,  a first  order 
correction  is  proposed  by  assuming  a small  charge  accumu- 
lation at  the  drain,  which  depends  on  the  electric  field 
as  follows: 

-(CHADEX) (Edrgin  - i) 

n = Nd  [1+(CHAADR) {1-e  tc  }]  (3-95) 

CHAADR  is  a constant  representing  the  maximum  relative 

charge  accumulation  (for  Earaj^-*10)  » and  CHADEX  is  an 

electric  field  constant.  Equation  3-95  is  used  instead 

of  Nd  in  the  expression  for  the  drain  conduction  current 

in  the  neutral  channel  in  mode  A.  E.  is  the  electric 

drain 

field  at  (x=lg  y=0) . Thus,  this  correction  affects  the 
solution  for  V^,  its  derivatives  and  the  expressions  for 

gm'  gd* 
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CASE  NO.  3 


Figure  3.2  False  negative  slope  in  ID~VDS  curves 


-90- 


The  constants  CHAADR,  CHADEX  are  set  so  as  to  eliminate 
the  negative  slope  region.  Since  the  assumption  of 
equation  3-95  is  only  a first  order  correction  there  is  no 
unique  way  to  set  the  constants.  A convenient  and  effective 
method  to  do  it  is  as  follows: 

(a)  Choose  a value  of  VDg  3-4  times  the  "knee"  voltage. 

(b)  Set:  CHADEX  « 2. 

(c)  For  the  above  V__  plot  g„  versus  (zero  to 

ds  * m SG 

cuttoff) . 

(d)  Repeat  (c)  with  several  values  of  CHAADR 
(typically  0.05-0.2)  and  choose  that  value  which  makes 
gm  almost  constant  in  modes  A,  B (slightly  decreasing 
with  VgG) , and  then  sharply  decreasing  for  mode  c. 

(e)  With  the  value  of  CHAADR  set  in  (d)  plot  g versus 

m 

VDS  ^VSG=0*  ^°r  VDS=0  to  bey°nd  tbe  "knee". 

(f)  Repeat  (e)  with  several  values  of  CHADEX 

(typically  0. 5-1.0)  and  choose  that  value  which  makes 

g„  increase  with  V_e  up  to  the  "knee"  and  then  remains 
m ds 

almost  constant.  If  the  plot  exhibits  a maximum,  choose 
a value  that  minimizes  that  maximum. 

An  example  of  applying  the  above  method  is  given  in 
chapter  5 . The  correction  described  here  was  found  quite 
effective  for  typical  microwave  devices.  For  very  thick 
devices  the  correction  is  ineffective. 
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3.7.3  Smooth  Velocity  - Field  Curve 

As  stated  in  chapter  2,  in  the  entire  theory  presented 
above,  the  assumption  was  that  the  functional  dependence 
of  the  change-carrier  velocity  on  the  electric  field  is  a 
piecewise  linear  function  as  shown  in  Figure  2-2 (b) -curve  1. 
This  assumption  was  found  to  cause  some  "cosmetic"  problems 
in  operational  mode  A when  the  channel  is  thick  enough 
(VgG«0)  . As  discussed  in  the  previous  section  the  rotation 
of  the  velocity  vector  is  quite  small  in  the  neutral  channel, 
since  the  electric  field  is  almost  entirely  in  the 
x-direction.  Thus  the  problem  is  effectively  one 
dimensional  in  the  channel.  Therefore,  if  most  of  the 
conduction  current  comes  via  the  channel  (only  a small 
contribution  from  the  transition  region),  the  ID“VDS  curve 
is  practically  a replica  of  the  velocity-electric  field 
curve.  So,  the  ID“VDS  curve  of  the  FET  has  a "sharp  corner" 
for  VgG k 0 due  to  the  "sharp  corner"  in  the  assumed  v-E 
curve  (see  Figure  3.3) . 

From  experimentally  measured  I -Vnc  curves  of  FETs  it 
is  known  that  the  transition  between  the  two  sections  at  the 
"knee"  is  smooth.  To  get  such  a smooth  transition  for  the 
model  it  is  necessary  to  assume  a smooth  v-E  curve.  In 
modes  B and  C the  transition  is  smooth  even  for  the 
piecewise  linear  v-E  curve  since  there  is  no  neutral  channel 
at  the  drain,  and  the  effect  of  gradual  rotation  of  the 
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velocity  vector  creates  a smooth  transition.  Therefore, 
it  is  necessary  to  use  a smooth  v-E  curve  only  for  mode  A. 
The  smooth  curve  assumed  here  is  the  one  shown  in 
Figure  2-2 (b)  - curve  2.  This  is  a parabolic  function  with 
the  same  initial  slope  as  the  piecewise  linear  curve. 
Physically,  this  means  the  same  low-field  mobility. 

Thus,  the  equations  for  mode  A were  re-derived  with 
the  above  smooth  v-E  curve.  For  simplicity,  the  smooth 
v-E  curve  was  applied  only  in  the  neutral  channel,  while 
in  the  transition  region  the  piecewise  linear  curve  was 
kept.  Also,  since  the  problem  is  practically  one 
dimensional  in  the  neutral  channel,  the  smooth  curve  was 
assumed  to  hold  between  v -E  . The  I_-Vnc  curves  predicted 
by  the  model  with  the  smooth  v-E  curve  are  shown  in 
Figure  3.4. 

3.7.4  Criteria  For  Operational  Mode  And  Case 

As  mentioned  above,  the  FET  can  operate  in  either  one 
of  the  three  operational  modes  (or  cutoff  - referred  to  as 
mode  D) , and  in  each  one  of  these  modes  it  can  operate  in 
either  one  of  three  cases  (see  section  3.2).  Given  a pair 
of  voltages  VgG,  VQS  the  first  step  is  to  determine  the 
operational  mode  and  case,  and  then  apply  the  proper  set  of 
equations.  In  this  section  the  method  to  determine  mode  and 
case  is  outlined. 


actuuocuiv) 
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Figure  3.4  "Smooth"  ID"VDS 
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The  transition  between  mode  C and  mode  D (cutoff)  is 
defined  by  zero  conduction  current  in  both  source  and 
drain.  From  Equations  3-28  and  3-38  it  is  concluded  that 
both  currents  are  zero  if  a)  V =0,  and  b)  Vp=-8.  However,  Vp 
can  be  written  as  (Equation  2-31) : 

00 

VP  - VP  “ VSG  <3-96) 

where : 

00  2 12  2 
vp  = 2~[  (a-dr+d(a-d)  + (|  - ^j)d z]  - (3-97) 

Thus,  the  transition  between  modes,  c,  d occurs  for: 

00 

VSGI  = VP  + 8 (3-98) 

c-*-d 

Equation  3-98  defines  a simple  boundary  which  depends 
on  VgG  only.  The  transition  between  modes  B and  C is 
defined  by  d^(0)=0.  Using  the  equation  for  d-^O) 

(see  Equation  2-30)  get: 

a1(0).  a-|-  J|h^(vsg+*b)  + <^  - i)d2  - 0 (3-99) 

or  after  rearrangement  and  use  of  Equation  3-97: 

00 

VSG|  = vp  (3-100) 

b-k: 

The  transition  between  modes  A and  B is  defined  by 
d1(lg)=0.  Using  Equations  2-30  and  3-97  get: 
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SG 


+ Vll 

A-*-B 


00 


= VT 


(3-101) 


Note  that  3-101  is  different  from  the  other  two 
transition  conditions  3-98  and  3-100  in  that  it  depends  on 
both  VgG  and  VDg  (via  V^) . Since  is  unknown  a priori,  if 
the  pair  VSG ' VDS  is  such  that  the  operational  mode  is  found 
to  be  either  A or  B an  estimate  must  be  used  for  to 
decide  which  equations  to  use.  This  is  obtained  by  using 
the  approximate  expression  for  V1  in  mode  A.  Applying  the 
above  yields  the  approximate  expression 


V1  - 


2?A1£  - D0VPS 

di  - D0 


(3-102) 


where: 


3.256  - 1.063  ~ 


(3-103) 


= i - e""2  * 414<2 


0.3835 


(3-104) 


ir2eV 


iryo(d1  (0)+DOT] 

2av  sinh  (5^2.) 

S Za 

I.07y()[d1(0)+D1T] 


-3, 


.16 


, N_ (cm  J) -10J 

DOT  = | (1.16  + 0.12  -5 — ) RDO 

* 9 1016 


(3-105) 


(3-106) 


(3-107) 


(3-108) 


U' 
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N (cm  3)-1016 

D1T  = S— (1 .16+0 . 12  — ^ ) (6  + sin6)  (3-109) 

2ir  9-1016 


S’  = £ + e2) 


(3-110) 


Q1  is  the  solution  of  ©1  - sin01  = S.^  (3-111) 


02  is  the  solution  of  02  “ sin02  = 


(3-112) 


from  equation  (3-102)  is  used  in  equation  (3-101)  to 
determine  the  operational  mode.  However,  since  this  is  an 
approximation,  when  the  calculations  are  finished  the  exact 
Vx  is  inserted  in  (3-101)  and  if  the  initial  decision  is 
found  wrong,  the  calculations  are  repeated  with  the  correct 
mode  equation.  This  can  happen  only  for  driving  conditions 
close  to  the  transition  between  modes  A,  B. 

To  sum  up,  the  criteria  for  determining  the 
operational  mode  are: 


Vgc  > Vp  + 0 -*■  mode  D 


V + B > VSG  * VP 


mode  C 


VSG<  VP 


Vgc  + > Vp  -*■  mode  B 


VgG  + < Vp  -*■  mode  A 
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Once  the  operational  mode  is  determined  it  is  still 

necessary  to  determine  the  appropriate  case 

(see  section  3.2).  For  a typical  microwave  device  the  most 

common  one  is  case  2,  namely,  E.  . > E and 

drain  c 

Esource<  V Case  1 (Esource,Edrain  * Ec>  oocurs  only  for 
very  low  VDS(0.5-1V),  and  case  3 (Esource_Edrain  > *„> 

occurs  quite  rarely  (typically  for  ratio  small) . 

d 

From  the  equations  in  appendix  8.1  it  is  obvious  that 
the  electric  field  at  the  drain  (x=lg,  y=0)  is  approximately; 


"drain 


V_1  . 
lg  2a 


(3-113) 


For  very  low  V_0  : E,  . < E (case  1) . As  V__ 

2 DS  dram  c DS 

increases,  both  and  Vq  increase  and  eventually 

c 

Edrain  * Ec  ^case  2)  • An  estimate  to  this  critical  VDg  is 
achieved  by  neglecting  the  VQ  term  (V^  « VDS) • 

VDS  ' Vlg  <3-114' 

c 

VDS  is  in  fact  an  upper  limit  to  the  critical  VDg, 

since  actually  Vq  is  not  zero  and  it  contributes  to  the 

c 

electric  field,  so  that  it  reaches  Ec  for  VDg  < VDg . In 

view  of  the  above,  the  method  to  determine  the  proper  case 

% 

is  as  follows: 
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a) 

If 

VDS 

c 

i VDS 

choose  case 

1 

b) 

If 

VDS 

c 

> VDS 

choose  case 

2 

c) 

Proceed 

the  calculation  of 

V..  with  the  chosen  case 

equations.  When  V1  is  calculated  compute  the  exact 
electric  field  at  the  drain. 


d)  If  Edrain  is  compatible 
proceed  to  (e) . Otherwise, 


with  the  case  chosen 
repeat  (c)  with  the  proper 


case  equations. 


e)  Calculate  the  exact  source  electric  field. 


If  Esource  < Ec  Proceed  with  device  properties 
calculations.  Otherwise,  repeat  the  calculations  with 
case  3 equations. 
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4 . THE  COMPLETE  FET  COMPUTER  MODEL  STRUCTURE 

4 . 1 INTRODUCTION 

The  modeling  of  the  active  part  of  the  FET  (under  the 
gate  metalezation) , which  is  the  most  important  part,  was 
described  in  some  detail  in  chapter  3 . In  the  present 
chapter  the  complete  model  of  the  FET  is  presented.  This 
model  includes  in  addition  to  the  active  region,  various 
parasitic  elements,  some  of  which  are  internal  to  the 
device,  and  some  of  which  are  external.  The  evaluation  of 
these  elements  is  discussed,  and  equations  are  derived  for 
the  internal  elements.  Also  considered  in  this  chapter  is 
the  structure  of  the  computer  model,  which  is  constructed  in 
modular  form.  The  various  subroutines  are  described  and 
their  function  discussed.  The  input  and  output  parameters 
are  stated  and  their  meaning  explained.  The  use  of  the 
model  for  various  applications  is  described. 

4.2  THE  MODEL 

The  complete  model  of  an  assembled  FET  chip  is  given 
in  Figure  4.1.  The  box  labeled  "basic  FET  chip"  represents 
the  active  region  of  the  device,  and  is  characterized  as 
described  in  the  previous  chapter «.  As  far  as  the  external 
terminals  are  concerned  the  equations  characterizing  the 
"chip"  are: 
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= GVSG 


+ GVDS 


dV  dV 

:con  + DVSG  -ai2  + DVDS  -&r 


dV„_  dV 

= Icon+SVSG  + SVDS  -d^ 


(4-1) 


(4-2) 


(4-3) 


Idis 

The  internal  source  current  is  designated  Is  in 

int 

equation  (4-3)  to  distinguish  it  from  the  total  source 

current,  which  includes  also  the  diode  current  (to  be 

discussed  below) . Of  course,  since  I = I_  - I , the 

s D g 

following  hold:  *nt 


SVSG  = DVSG  - GVSG 


SVDS  = DVDS  - GVDS 


(4-4) 

(4-5) 


The  conduction  current,  I , and  the  displacement 

con 

current  coefficients  (GVSG,  GVDS,  etc.)  are  calculated  by  the 
model  as  functions  of  Vg^,  VDg  using  the  expressions  derived 
in  the  previous  chapters.  For  example,  GVSG,  which  is  the 

dVSG 

coefficient  of  ^ - in  the  gate  current  equation,  is  found 
from  appendix  8.2  for  mode  A (Equation  8.2-11): 


GVSG  = CQ  (KGS  + KG*K1) 
(mode  A) 


(4-6) 


I 
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The  other  network  elements  in  Figure  4.1  are  basically 
parasitics,  which  are  unwanted  but  also  unavoidable.  The 
external  parasitics  are  as  follows: 

input  stray  capacitance 
output  stray  capacitance 
gate  wire  bond  inductance 
drain  wire  bond  inductance 
source  wire  bond  inductance 

The  wire  bond  inductances  can  usually  be  estimated  by  the 
physical  size  of  the  wire  bond.  The  stray  capacitances  are 
more  difficult  to  estimate.  For  a packaged  device  they  can 
be  characterized  by  measuring  an  empty  package. 

Rg  is  the  gate  metallization  resistance.  Usually  it 
can  be  estimated  quite  accurately  from  the  geometry  of  the 
gate  electrode,  as  well  as  the  metallization  thickness  and 
metal  properties.  The  resistance  Rg  as  well  as  the 
external  parasitics  are  inputs  to  the  model,  namely,  they 
have  to  be  estimated  and  specified  by  the  user. 

The  diode  in  Figure  4.1  represented  the  gate-source 
Schottky  barrier  junction.  Usually  the  gate-source 
junction  is  reverse-biased,  and  the  diode  is  not  necessary. 
However,  in  some  cases  (such  as  overdriven  power  amplifiers 
or  frequency  multipliers,  oscillators,  etc.)  the  gate-source 
voltage  becomes  positive  for  a portion  of  the  period,  and 
there  is  gate  conduction.  To  simulate  this  gate  conduction 


sm 


SO 

Lg 

L . 
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the  diode  was  added  to  Figure  4.1.  This  approach  (an 

external  diode)  is,  of  course,  approximate  and  does  not  take 

into  account  any  effect  internal  to  the  device  due  to  gate 

conduction.  The  characterization  of  the  diode  is  up  to  the 

user.  Throughout  the  simulations  presented  in  chapter  5 

an  ideal  diode  was  assumed.  The  diode  is  assumed  to  be  a 

short-circuit  for  V. . . > V and  an  open  circuit  for 

diode  on 

Vdiode  - Von*  The  turn-on  voltage,  VQn,  is  specified  by  the 
user  (a  typical  number  is  0.5  V).  In  all  practical  cases 
the  drain  to  gate  voltage  is  positive,  so  the  drain-gate 
diode  is  excluded  from  Figure  4.1.  In  any  special  case,  in 
which  VDG  may  be  negative  it  is  necessary  to  add  this  extra 
diode  to  the  model  (see  oscillator  example  in  chapter  5) . 

The  two  resistors  Rg<j,  Rgs  represent  the  bulk 
resistance  of  the  epitaxial  layer  from  the  edge  of  the  gate 
electrode  to  the  drain  and  source  electrodes,  respectively 
(see  Figure  1-4) . These  resistors  are  calculated  by  the 
simple  formula: 


R 


(4-7) 


where,  the  conductivity  is: 


a = qNDyo 


(4-8) 
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The  cross-sectional  area  A is: 

A = W a (4-9) 

The  effective  length  of  the  resistor  is  estimated  as 
the  distance  between  the  electrodes  plus  one-quarter  of  a 
circle  of  diameter  a.  Thus: 


lgs+  JS. 


ss  qNdyQWa 


sd 


lgd+  J2. 
qN^Wi" 


(4-10) 


(4-11) 


The  resistances  Rss,  Rsd  are  not  computed  by  the  FET 
computer  model  package.  They  must  be  calculated  by  the 
user  in  the  main  program  using  Equations  4-10,  4-11. 

The  resistance  Rsubst  is  the  substrate  resistance. 
Since  the  semi-insulating  substrate  on  top  of  which  the 
epitaxial  layer  is  grown  is  not  ideal  and  has  some 
non-zero  conductivity,  some  of  the  device  current  is 
carried  via  the  substrate.  In  the  present  model,  this 
phenomenon  is  not  included  physically.  Instead,  Rgubst  is 
added  externally  as  in  Figure  4.1.  The  value  of  Rsubst 
must  be  specified  by  the  user.  For  good  substrates  this 
effect  is  negligible,  and  a large  value  should  be  specified 

for  R . . . 

subst 
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R is  an  effective  charging  resistor  through  which  the 
c 

depletion  region  is  charged  and  discharged.  R is  in  series 
with  the  gate  to  source  capacitance,  thus  only  the  displace- 
ment source  current  is  carried  by  it,  as  shown  in 
Figure  4.1.  The  calculation  of  Rc  is  described  in  the  next 
section. 

4.3  THE  CHARGING  RESISTOR 
4.3.1  Introduction 

As  discussed  in  chapter  2 (section  2.2) , the  boundary 
condition  2-9,  which  is  used  to  evaluate  the  parameter 

(source  and  drain  conduction  currents  equal) , is  strictly 
valid  for  DC  drive  only.  For  AC  drive,  the  size  of  the 
depletion  region  changes  with  time,  and  it  is  charged  and 
discharged.  The  charge  necessary  to  accomplish  the  above  is 
extracted  from  the  current  between  source  and  drain.  Thus, 
instantaneously,  the  source  and  drain  conduction  currents 
are  not  equal,  but  there  is  a time  delay  due  to  the  finite 
velocity  of  the  charge  carriers.  Therefore,  to  compensate 
at  least  partially  for  the  error  introduced  by  2-9,  two 
parameters  are  introduced  into  the  model,  which  are 
calculated  for  each  pair  of  VSG*  VDS*  Parameters  are: 

a)  the  time  delay  (transit  time  via  the  channel)  and  b)  the 
effective  charging  resistor  through  which  the  depletion 
region  is  charged.  The  calculation  of  the  time  delay  and 
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its  use  in  the  model  are  presented  in  section  4.4.  The 
calculation  of  the  charging  resistor  is  described  here. 

4.3.2  Charging  Resistance  Of  a Distributed  Network 

The  FET  channel  is  basically  a distributed  combination 
of  resistors  and  capacitors  as  shown  in  Figure  4.2(a). 

Part  of  the  charge  carriers  flowing  from  source  to  drain  are 
used  to  charge  the  distributed  capacitors.  To  investigate 
this  "leakage"  current  consider  the  discrete  distributed 
network  in  Figure  4.2(b).  For  simplicity,  at  this  stage 
the  channel  is  assumed  to  be  subdivided  such  that  all  the 
capacitors  are  equal.  This  assumption  is  removed  later. 

For  the  network  in  Figure  4.2(b)  assume  that  the 
frequency  o)Q  is  not  too  high,  so  that  the  voltage  across 
each  capacitor  is  approximately  equal  to  VM  (small  voltage 
drop  across  the  resistors) . Then  the  network  in 
Figure  4.2(c)  can  be  considered  equivalent  to  the  network 
in  Figure  4.2(b).  The  effective  charging  resistor,  Rc,  is 
computed  as  follows: 

Subject  to  the  assumption  above  (u)Q  low  enough)  , the 
current  in  each  capacitor  is  I =wC  V„.  Then,  the  current 
in  Rjj  is  IQ,  in  RN_1  is  2IQ...in  is  NIQ.  So,  the  total 
power  dissipated  in  the  network  is: 

2 2 
P “ !0  + 4rN_1+9Rn_2+“*+n  Ri>  (4-12) 


•9 


(0 


Figure  4.2  (a)  FET  channel  as  a distributed 

network  (b)  discrete  distributed 
network  considered  in  text  (c)  equi 
valent  network 
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However,  the  power  dissipated  in  the  equivalent  network  is 

PE  = (IoN,2rc  - "W  (4-13) 

Equating  P=PE  one  gets  an  expression  for  Rc: 


(4-14) 


A similar  result  can  be  achieved  now  by  returning  to 
the  continuous  case.  In  the  discrete  case  it  was  shown 
that  the  current  in  each  resistor  is  larger  the  further  to 
the  left  it  is.  Assuming  now  a homogeneously  distributed 
capacitor  in  the  continuous  case,  the  current  is  simply  a 
linearly  decreasing  function  of  distance: 

I(x)  = Ig  (1-  ^)  (4-15) 

where  Ig  is  the  generator  current. 

For  the  case  of  a non-homogeneously  distributed 
capacitor  the  current  is  given  by  the  capacitance  ratio: 

I(x)  = Ig  (1-  ^1)  (4-16) 

where  Cfc  is  the  total  channel  capacitance  and  C(x)  is  the 
portion  of  channel  capacitance  between  0 to  X (x=0  at 
the  source) . 
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Using  the  same  approach  as  in  the  discrete  case,  the 
power  dissipated  in  the  network  is  calculated.  The  power 
dissipation  in  the  small  portion  of  the  channel  between 
X and  X+dX  isr 

2 2 

dP  = Ig  (1-  R(x)  dx  (4-17) 

ct 


R(x)  is  the  channel  resistance  per  unit  length. 
Thus, 

ig 


p - Ig 


/ 


(1-  £i£L)  R (x)  dx 


(4-18) 


In  the  equivalent  network  (Figure  4.2(c)): 


p*  = ig  R 


(4-19) 


Equating  P=P£  get: 


R 


c 


2R  (x)  dx 
ut 


(4-20) 


Equation  (4-20)  gives  the  general  expression  for  the 
charging  resistor  of  the  channel.  To  actually  evaluate  the 
resistance  the  proper  expressions  for  C (x) , Cfc,  R(x)  must  be 
used  in  4-20.  The  above  expressions  are  different  for 
each  operational  mode,  therefore,  the  charging  resistor 
expression  is  different  for  each  mode. 
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For  the  purpose  of  calculating  C(x),  and  R(x)  for 
each  mode,  the  FET  is  divided  into  channel  and  depletion 
regions.  The  transition  region  is  not  v.  nsidered  here. 
Actually,  half  of  it  is  "attached"  to  the  channel  and  the 
other  half  to  the  depletion  region,  as  demonstrated  in 
Figure  4.3.  Thus  the  effective  channel  height  is 
di(x)+  and  the  rest  of  the  device  is  assumed  to  be 
depletion  region. 

4.3.3  Mode  A 

The  capacitance  of  an  element  of  the  depletion  region 
between  X and  X+dX  is 


thus, 


dC 


eW 


dx 

a-d1(x)-d/2 


C (x) 


dx 

a-d1(x)-d/2 


(4-21) 


(4-22) 


and  so  Cfc  = C(lg) 


(4-23) 


Using  Equation  2-30  for  dx (x) , the  integral  in  (4-22) 
is  solved,  the  function  C (x)  determined  and  the  ratio  is 
found: 


1 


C(x)  _ J VVC°10  - Jl5*  Wsg+cmo 

°t  Jv^V~KOlO  - J VSG  + CD10 


(4-24) 
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Figure  4.3  FET  subdivision  for  charging 

resistor  calculation  (a)  mode  A 
(b)  mode  B (c)  mode  C 


The  channel  resistance  per  unit  length  R(x)  is 


calculated  by  definition: 


R(x) 


1_ 

aA 


1 

qNDWy [|  +d1 (x) ] 


(4-25) 


a - conductivity,  A - cross-sectional  area. 

Due  to  the  non-linear  relationship  between  electron 

velocity  and  electric  field,  the  dynamic  mobility  y to  be 

used  in  4-25  is  electric  field  dependent.  Using  the 

piecewise  linear  V-E  curve  (Fig.  2.2(b) -curve  1)  and 

defining  X . the  point  at  which  E (x  ) =E  : 

o o c 


P = 


X < x_ 

- o 


X > x_ 

o 


(4-26) 


Inserting  4-24,  4-25  in  4-20  yields  the  value  of 

Rc.  In  the  computer  model  the  integral  is  evaluated 

numerically  (10  divisions) . The  calculation  of  X is 

o 

described  in  section  4.4. 

4.3.4  Modes  B And  C 

The  evaluation  of  R„  for  modes  B,  C follows  the  same 

c 

steps  as  described  above  for  mode  A.  Listed  below  are  the 
expressions  derived  and  comments  regarding  specific 
problems . 
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The  problem  in  mode  C is  that  the  channel  is  partially 
depleted.  It  gets  more  and  more  depleted  as  V__  increases 
up  to  cutoff,  at  which  point  the  channel  disappears.  Thus, 
a criterion  has  to  be  established  for  a boundary  between 
the  channel  and  the  depletion  region.  Since  the  drain  is 
always  more  depleted  than  the  source,  it  is  possible  to  find 
a point  X < lg  at  which  the  charge  carrier  density  is  so 
small  that  it  can  be  practically  considered  as  fully 
depleted.  Thus,  for  the  purpose  of  the  charging  resistor 
calculation  the  channel  length  is  Xc  rather  than  lg. 

Xc  is  defined  conveniently  as  the  point  at  which 
n(x)  = 0.1  Nd.  Using  Equation  4-30  this  yields: 

xc  = ^ (0.93  + Vp)  (4-32) 

Thus,  for  X > lg  the  integral  for  R calculation  is 
c c 

evaluated  from  0 to  lg.  However,  for  X < lg  the  upper 

c 

limit  for  the  integral  is  Xc>  This  shortening  effect  is 
incorporated  in  mode  B as  well  since  in  mode  B there  is  also 
partial  depletion  in  the  range  y < X < lg. 

4 . 4 TIME  DELAY 

Due  to  the  finite  velocity  of  the  charge  carriers  in 
the  channel,  an  electron  entering  the  channel  at  the  source 
side  travels  for  a time  period  t before  it  reaches  the 
drain.  This  time  period  is  the  time  delay  or  transit 
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time,  and  is  a basic  limitation  of  the  device.  Of  course, 

the  shorter  the  gate,  the  smaller  is  t.  For  a typical 

microwave  device,  t is  on  the  order  of  10  psec. 

To  calculate  t it  is  first  necessary  to  calculate  the 

point  X at  which  E (X  ) =E„  along  the  X axis.  From  the 
o o c 

analytic  solution  expressions  of  chapter  2,  the  electric 
field  along  the  x-axis  is  given  by 


V irv  cosh(g-) 

V*'°>  ’ I?  + 1ST 


sinh(^L) 


Set  E (X  , 0)  = E and  get  an  expression  for  X 

AW  C U 

v 

2a  (V  la-)  2a*sinh  (£^0 

Xq  = |«  arc  cosh  [ -c  ^ ^— ] 


(4-33) 


(4-34) 


Recall  that: 

, 1/2 

arc  cosh  (X)  = ln[X  + (X  - 1)  ] (4-35) 

So,  the  channel  is  subdivided  into  two  regions: 

1)  XQ  > X > 0 : vx  = P0Ex  constant  mobility 

2)  X > Xn  : v = v constant  velocity 

The  transit  time  for  region  2 is  simply 


t2 


lg-XQ 


V 

s 


(4-36) 
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In  region  Is 


dx 

x 5t 


V--  = tt  = u e 


O X 


dx 

y E 
o x 


= dt 


(4-37) 


(4-38) 


(4-39) 


Equation  4-39  is  used  to  calculate  by  inserting  in 
it  E^  from  equation  4-33.  The  resulting  integral  can  be 
evaluated  analytically  by  using  integral  tables.  The 
solution  is: 


"X_ 

2a  r:r  . B+A  cosh(^) 

— - arC  310  irv  1 

o 


irx0  J rz 

A+B  cosh(T-H.)  \/b  -A 


2 2 
for  B >A^ 


2a 

1 \*y. 


In  [ 


r~3 — t ™ 

A+B+^/A  -B^tanh  (^i 


2a 

) 


A+B 


tt. 

tanh (^a  ) 


Ja2-b2 


1 for  A2>B2  (4-40) 


2alg 

^oVl 


nx 

tanh  <*1T>  for  A=B 


where: 


(4-41) 


7” 
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(7TVQ/2a) 

sinh(^£) 

The  total  transit  time  is: 


(4-42) 


T 


T1  + T2 


(4-43) 


In  modes  B and  C the  shortening  effect  of  the  channel 

(section  4.3.4)  is  taken  into  account  and  the  transit  time 

is  calculated  from  x=0  to  X rather  than  to  lg. 

c 

As  obvious  from  the  above,  the  time  delay  x is  a 
function  of  the  driving  voltages  V0_,  V_ _ . However,  from 
the  actual  calculations  it  was  found  that  x is  nearly 
constant  over  most  of  the  dynamic  range.  To  incorporate 
the  time  delay  in  the  device  equations  4-1  to  4-3,  note 
that  the  main  effect  of  this  phenomenon  is  to  delay  the 
electrons  entering  at  the  source  end.  Thus  in  equation  4-2 
the  conduction  current  term  has  to  be  delayed  by  x . 

So,  the  instantaneous  device  equations  are:  ^ j 


ig(t) 

vt} 


dvcp  dv 

GVSG(t)  + GVDS(t) 


dV 


rCon  (t_T)  + DVS<3(t) 


+ DVDS (t) 


dV 


DS 


dt 


(4-44) 

(4-45) 


Is  (t)  = In(t)  - Ig (t) 
int  D 


I 


1 

4 


(4-46) 
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When  t is  incorporated  in  the  equations  as  above, 
equations  4-4,  4-5  are  invalid,  and  the  source  current  must 
be  calculated  by  equation  4-46. 

4.5  COMPUTER  PACKAGE  STRUCTURE 

The  new  FET  model  presented  in  this  report  is 
practically  implemented  in  the  form  of  a computer  subroutine 
package.  The  user  has  to  write  his  own  main  program,  which 
is  to  reflect  the  particular  network  to  be  analyzed  or 
snythesized.  In  this  section,  the  structure  of  the  package 
is  presented  as  well  as  the  "communication"  between  the 
package  and  the  user's  main  program.  The  package  is 
written  in  FORTRAN  IV. 

The  following  statements  must  appear  in  the  calling 
program: 

COMMON/ INDAT/ND , MUO , A , LG , PHB , W, EC , ER , VS 

COMMON/FETCUR/ ICON , GVSG , GVDS , DVSG , DVDS , SVSG , SVDS , CINCS , 
*COCS , CINCG , COCG , CINCD , COCD , GM , GDD , TAU , RCHARG , RCHEFF 
COMMON/NOIT/KKVSA , KKVSB , KKVSC , INVSAT , KNNVSA , KNNVSB , 

* KNNVSC , KKNLDA , KKNLDB , KKNLDC 
COMMON/CHARAC/CHAADR , CHADEX 
REAL  ND , MUO , LG , ICON 

The  common  INDAT  contains  the  parameters  describing 
the  device  geometry  and  metallurgical  properties.  These 


I 
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parameters  must  be  specified  by  the  user  before  any  of  the 
package  subroutines  are  called.  All  physical  quantities 
should  be  specified  in  MKS  units,  and  so  are  the  output 
variables.  The  meanings  of  the  symbols  are: 

Nd  - doping  level,  MUO=yQ  - low  field  mobility  (does 
not  have  to  be  specified  by  the  user  - it  is  calculated  as 

vs 

y = 5—)#  A - device  thickness,  LG  - gate  length,  PHB  = <J> 

0 ° 

- built-in  potential,  W - device  width,  EC  - critical 
electric  field  for  velocity  saturation,  ER=er  - relative 
dielectric  constant  («12.5  for  GaAS)  , VS  - saturated 
velocity. 

The  common  FETCUR  contains  the  analysis  results  as 
follows:  ICON  - conduction  current,  GVSG  to  SVDS  - 
displacement  current  coefficients  (equations  4-1  to  4-3) , 
CINCS=GVSG  - input  capacitance  (common  source) , 

COCS— COCG=DVDS  - output  capacitance  (common  source  and 
gate) , CINCG=SVDS-SVSG  - input  capacitance  (common  gate) , 
CINCD=CINCS  - input  capacitance  (common  drain) , 

COCD=CINCG  - output  capacitance  (common  drain) , 

GM  - transconductance,  GDD  - drain  conductance,  TAU  - time 
delay,  RCHARG  - charging  resistor, 

RCHEFF=RCHARG+TAU/ (2tt  *CINCS)  - effective  resistance,  which 
simulate  both  the  charging  and  the  time  delay  effects. 
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Upon  exit  from  the  package  the  above  parameters  are 
defined,  and  may  be  employed  by  the  user  in  the  network 
equations.  The  common  NOIT  contains  variables  defining  the 
number  of  numerical  iterations  to  be  performed  for  each  mode 
and  case  (section  3.2)  as  follows: 

KKVSA  - mode  A case  3,  KKVSB  - mode  B case  3,  KKVSC  - 
mode  C case  3,  KNNVSA  = mode  A case  1,  KNNVSB  - mode  B case 
1,  KNNVSC  - mode  C case  1,  KKNLDA  - mode  A case  2,  KKNLDB  - 
mode  B case  2,  KKNLDC  - mode  C case  2.  The  above  9 
parameters  must  be  specified  by  the  user  before  calling  the 
package  subroutines.  This  can  be  done  conveniently  by 
subroutine  SETITR  to  be  described  below.  The  variable 
INVSAT  contains  a code  number,  which  specifies  the  mode  and 
case  of  the  last  analysis  performed.  Thus,  INVSAT  is  an 
output  parameter.  The  value  of  INVSAT  returned  for  each 
mode  and  case  is  listed  in  table  4.1.  For  mode  D (cutoff) 
INVSAT=4 . 

The  two  parameters  in  common  CHARAC  are  related  to  the 
charge  accumulation  correction  (section  3.7.2)  as  specified 
in  equation  3-95.  These  parameters  must  be  specified 
before  calling  the  package  subroutines. 

In  addition  to  the  above  commons,  the  package  employs 
several  more  commons  which  are  used  internally.  The  user 
is  prohibited  from  using  commons  bearing  those  names.  The 
names  are:  ACFCHC , SOURC , DRCOM, DRAIN, TRSTIM, CALC02 , UPCOM, 
CALCON , SMWDEC . 
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Table  4.1  Value  of  INVSAT  for  each  mode  and  case 


mode 

case^^_^ 

A 

B 

C 

1 

10 

20 

30 

2 

100 

200 

300 

3 

1 

2 

3 

The  FET  package  is  composed  of  14  subroutines,  out  of 
which  only  three  are  called  by  the  user.  The  subroutines 
called  by  the  user  are: 

SETITR  - this  subroutine  sets  the  number  of  iterations 
to  be  employed  in  the  computations  for  each  mode  and  case, 
namely,  it  sets  the  parameters  in  the  common  NOIT.  The 
calling  sequence  is: 

CALL  SETITR (JJ) 

JJ  is  a vector  of  length  9,  which  has  to  be  specified 
before  calling  SETITR.  The  first  3 locations  of  JJ 
correspond  to  case  1:  mode  A,B,C,  the  next  three  correspond 
to  case  2:  mode  A,B,C,  and  the  last  three  correspond  to 
case  3:  mode  A,B,C. 

FIRST  - this  subroutine  calculates  all  the  necessary 
constants  used  in  the  analysis.  It  has  no  parameters  and 
the  calling  sequence  is  simply: 

CALL  FIRST 

The  subroutines  SETITR  and  FIRST  must  be  called  in  the 


A 


main  program  (in  this  order)  before  the  actual  analysis 
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starts  (However,  before  calling  these  two  subroutines  the 
parameters  in  the  commons  INDAT,NOIT  and  CHARAC  must  be 
defined) . These  two  subroutines  are  used  only  once  before 
the  analysis  starts,  and  are  not  called  again,  unless  any  of 
the  device  parameters  are  changed  or  the  user  wishes  to 
change  the  number  of  iterations. 

The  actual  analysis  is  performed  by  calling  the 
subroutine  MAYN: 

CALL  MAYN (VSG,VDS,V1) 
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Figure  4.4  Flow  chart  diagram  for  the 
FET  computer  model 
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MOD A SC  - The 

analysis 

for 

operational 

mode 

A 

with 

smooth  V-E  curve, 

MODAPL  - The 

analysis 

for 

operational 

mode 

A 

with 

piecewise  linear 

V-E  curve. 

MODEB  - The 

analysis 

for 

operational 

mode 

B 

• 

MODEC  - The 

analysis 

for 

operational 

mode 

C 

• 

MODED  - The 

analysis 

for 

operational 

mode 

D 

(cutoff) . 

The  subroutine  MODED  is  used  just  to  set  the  output 
parameters  to  the  constant  values  corresponding  to  cutoff, 
which  are  the  limiting  values  of  mode  C parameters. 

In  addition  to  the  above,  there  are  four  auxiliary 
subroutines,  which  are  used  to  calculate,  numerically, 
integrals  needed  for  the  calculation  of  the  conduction 
current  and  its  various  derivatives,  as  discussed  in 
chapter  3.  Also,  there  are  two  subroutines  used  to  solve 
numerically  (Newton-Raphson)  two  simple  algebraic  equations. 
The  subroutines  are: 

ROOT  - solves  the  equation  x-sinx=A. 

ROOT  1 - solves  the  equation  x-sinx+B (x+sinx) =A 
UPADA  - computes  the  source  current  integral  and  its 
derivatives  of  the  transition  region  for 


modes  A and  B. 
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UP  AD  ER  - same  as  UPADA  except  for  mode  C. 

DRCURA  - computes  the  drain  current  integral  and  its 

derivatives  of  the  transition  region  and  , 

channel  for  mode  A. 

DRACUR  - same  as  DRCURA  for  modes  B and  C. 

The  new  FET  computer  model  is  valid  only  for  VDg  >0.  / 

If  a negative  value  for  VDg  is  used  when  calling  MAYN  the 
program  stops  with  a message  "VDS  NEGATIVE",  and  the  value 
of  VDS  is  printed.  This  is  not  a basic  limitation  of  the 
model,  since  VDg  < 0 simply  means  exchange  of  source  and 
drain.  Thus,  whenever  VDg  < 0,  the  user  can  exchange  the 
roles  of  drain  and  source,  such  that  the  new  VDg  is  positive 
and  then  use  the  model  (see  oscillator  example  in 
chapter  5) . 

There  is  no  upper  limit  in  the  range  of  V-._  (it  can 

bvj 

extend  to  °°)  . However,  for  VgG  < 0,  the  gate-source 
junction  becomes  forward  biased  and  gate  conduction  occurs. 

VgG  cannot  become  more  negative  than  the  turn-on  voltage 
of  the  (clamping)  diode.  This  effect  is  taken  care  of 
automatically  by  the  ideal  diode  used  in  the  model 
(Figure  4-1) . Since  the  diode  is  taken  external  to  the 
active  region,  the  computer  model  is  set  such  that  for 
VSG  < “ ^b'  VSG  set  e<3ual  to  $3  and  the  device  properties 
calculated. 
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5.  SIMULATION  RESULTS 

5 . 1 INTRODUCTION 

In  this  chapter  the  usefulness  of  the  model  is 
demonstrated  with  several  simulations.  The  results  are 
presented  in  the  form  of  graphs  generated  by  the  graphics 
package  of  the  PDP-11/45  mini  computer.  The  graphs  were 
plotted  on  a Tektronix  terminal  and  permanent  records  made 
with  a hardcopy  unit.  The  results  presented  here  contain 
the  following: 

a)  DC  analysis  - I-V  curves; 

b)  small  signal  analysis  - gm,  g d,  capacitances, 
S-parameters; 

and  c)  large  signal  analysis  - oscillator  build  up  and 
properties. 

The  effect  of  the  number  of  iterations  used  is  shown 
by  performing  analysis  of  the  same  device  with  varying 
number  of  iterations.  The  procedure  for  setting  the  charge 
correction  coefficients  is  demonstrated  on  a practical 
device. 

All  the  computations  presented  in  this  chapter  apply 
for  a typical  microwave  device,  the  parameters  of  which  are 
presented  in  Table  5.1.  The  first  12  parameters  correspond 
to  the  geometrical  and  material  properties.  The  last 
6 parameters  correspond  to  the  external  parasitics,  as  well 
as  the  characteristic  impedance  used  in  the  S-parameter 
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Table  5.1 


Typical  microwave  FET  parameters 
used  in  the  computations 


CASE  2 DEVICE  PARAMETERS: 

1 ) GATE-SOURCE  SEPERATION  ( LGS  ):  1 70  MICRONS 

) GATE  LENGTH  (LG):  1.  70  MICRONS 

) GATE-DRAIN  SEPERATION  ( LGD  ):  1.70  MICRONS 

) GATE  WIDTH  ( W ):  600.00  MICRONS 
) DOPING  LEVEL  ( ND  ):  7 50E+16  CM-3 

) CRITICAL  ELECTRIC  FIELD  (EC):  3.20  KV/CM 

) SATURATED  ELECTRON  VELOCITY  (VS):  1.  36E+07  CM/SE 

) RELATIVE  DIELECTRIC  CONSTANT  ( ER  ):  12.50 

) BUILT-IN  POTENTIAL  ( PHB  ):  0.700  VOLTS 
) GATE  METALLIZATION  RESISTANCE  ( RG  ):  2 000  OHMS 
) SUBSTRATE  LEAKAGE  RESISTANCE  ( RSUBST  ):  2.  00E+05 

) EPITAXIAL  LAYER  THICKNESS  (A):  0.3000  MICRONS 


CASE  2 PACKAGE /CIRCUIT  PARAMETERS: 
CHARACTERISTIC  IMPEDANCE  ( ZO  ):  50.0  OHMS 

SOURCE  INDUCTANCE  (SIND):  1.  00E-01  NANOHENRIES 

GATE  INDUCTANCE  ( GIND  ):  5.  00E-01  NANOHENRIES 

DRA I N I NDUCT ANCE  ( D I ND  ) : 2.  70E-0 1 NANOHENR I ES 

PARASITIC  INPUT  CAPACITANCE  ( SCAR  IN  ):  1.50E-01 

PARASITIC  OUTPUT  CAPACITANCE  ( SCAPO  ):  1.  50E-01 


P 


P 


OHMS 


I COFARAD 
I COFARAE 


I 
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calculations.  These  parameters  correspond  to  the  network 
components  in  Figure  4.1. 

5.2  NUMBER  OF  ITERATIONS 

As  described  in  chapter  3,  the  computer  model  first 
solves  for  the  parameter  V1  analytically,  and  then  employs 
numerical  iterations  to  improve  the  accuracy.  Of  course, 
since  the  numerical  iterations  are  time  consuming  it  is 
desirable  to  minimize  the  number  of  iterations  which  yield 
reasonable  results.  To  determine  this  minimum  it  is 
necessary  to  generate  I-V  curves  for  the  device  with  several 
number  of  iterations.  The  I-V  curves  change  as  the  number 
of  iterations  is  changed.  When  a reasonable  convergence  has 
been  achieved  the  change  in  the  curves  becomes  negligible 
as  the  number  of  iterations  increases.  The  smallest  number 
for  which  convergence  is  obtained  is  chosen  as  the  useful 
number  of  iterations.  The  number  of  iterations  for  each 
operational  mode  is  requested  by  the  program. 

The  above  procedure  is  demonstrated  for  the  device 
considered  in  this  chapter  for  operational  mode  A. 

Figure  5.1(a)  presents  the  I-V  curve  in  mode  A (VSG=0) 
without  any  corrective  iterations.  The  inaccuracy  in  the 
analytic  solution  causes  the  "jump"  in  the  curve  at  the 
point  where  the  set  of  equations  used  is  switched  from 
case  1 to  case  2.  Figure  5.1(b)  presents  the  I-V  curve 
for  the  same  case  with  one  corrective  iteration.  This 
curve  is  smooth  and  odes  not  exhibit  any  discontinuities. 
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It  was  further  found  that  adding  more  than  one  iteration  did 
not  cause  any  appreciable  change  in  the  curves 
(see  Figure  5.1(c)  - 4 iterations).  Thus  it  was  decided 
to  use  one  corrective  iteration  for  mode  A.  Using  the  same 
procedure  it  was  found  that  the  necessary  number  of 
iterations  for  mode  B is  one,  and  for  mode  C is  two. 

5.3  CHARGE  CORRECTION  COEFFICIENTS 

As  described  in  chapter  3,  the  computer  model  employs 
a charge  correction  in  operational  mode  A.  This  is  a first 
order  correction,  which  allows  a small  charge  accumulation 
at  the  drain  side  of  the  gate  active  region  (Equation  3-95)  . 
The  procedure  to  set  the  coefficients  CHAADR,  CHADEX  in 
Equation  3-95  is  described  in  section  3.7.2.  This 
procedure  is  demonstrated  for  the  device  considered  here. 

The  steps  taken  are  as  follows  (all  calculations  were 
performed  with  the  number  of  iterations  as  stated  in  the 
previous  section) : 

1)  choose  Vds=3V.  This  value  is  large  enough  compared 
to  the  "knee"  (=0.5V). 

2)  plot  gm  versus  VgG  for  the  above  value  of  VDg 
without  charge  correction  (the  two  coefficients  set  to 
zero)  - Figure  5.2(a). 

3)  set  CHADEX* 2 and  plot  gm  versus  V_„  (with  V_=3V) 

in  ovj  Do 

with  CHAADR=0 . 2 , 0.15,  0.1,  0.08.  The  "best"  curve 
(according  to  section  3.7.2  guidelines)  is  for  CHAADR=0.1 
(Figure  5.2(b) ) . 
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A 


5.2  gm  plots  (c)  CHAADR  = 0.1,  CHADEX  =0.5 


4)  Thus  choose  CHAADR=0 . 1 . This  value  yields  a smooth 


curve  for  g with  relatively  small  variation  of  g in 
m m 

modes  A,  B and  a sharp  decrease  for  mode  C (VgG  > 2.8V). 

5)  set  CHAADR=0 . 1 and  plot  g versus  V__  (with  Ve  =0) 

m DS  SG 

with  CHADEX=2,  1,  0.5,  0.3.  The  "best"  curve  (according 
to  section  3.7.2  guidelines)  is  for  CHADEX=0.5 
(Figure  5.2(c)). 

6)  Thus  choose  CHADEX=0 . 5.  This  value  yields  a smooth 
curve  with  a minimum  peak  at  the  knee. 

5.4  GENERAL  PROPERTIES 

Once  the  number  of  iterations  and  the  charge  correction 
coefficients  are  set,  the  device  is  uniquely  defined  as  far 
as  the  computer  model  is  concerned,  and  its  general 
properties  (DC  and  small  signal  AC)  can  be  calculated  and 
plotted.  An  example  of  the  interactive  "conversation" 
between  the  operator  and  the  computer  (PDP-11/45) , in  which 
the  necessary  data  is  entered,  is  given  in  Figure  5.3. 

The  first  data  item  entered  is  the  case  number.  The 
computer  disk  contains  a pre-prepared  disk  file,  which 
contains  all  the  necessary  information  for  a given  device 
(Table  5.1).  Each  device  is  assigned  a case  number.  Thus, 
when  the  case  number  is  entered,  the  device  parameters  are 
read  from  the  disk  file.  Next,  the  values  of  VDg  are 
entered  (number  of  values,  initial  value,  step) . Likewise 
the  values  of  VgG.  The  frequency  is  entered  only  for  the 
scattering  parameters  calculation.  When  no  scattering 
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Figure  5.3  "Conversation"  between  operator 
and  computer  PDP-11/45 
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parameters  are  calculated  this  data  is  item  not  needed. 

Next,  the  charge  correction  coefficients  and  the  number  of 

iterations  are  entered.  The  number  of  iterations  vector,  ( 

JJ,  contains  9 components  in  the  order  as  explained  in 

section  4.5.  Finally  enter  choices  for  plotting  (yes  or  no), 

type  of  plotting  (characteristics  or  scattering  parameters) , 

i 

and  line  printer  output  (yes  or  no) . 

The  results  of  the  analysis  corresopnding  to  the 
"conversation"  in  Figure  5.3  is  given  in  Figure  5.4(a). 

This  Figure  displays  the  I-V  curves  of  the  device  for  most 
of  its  dynamic  range.  More  properties  of  the  device  are 
presented  in  Figure  5.4(b)  and  5.5.  Figure  5.4(b)  displays 
the  drain  conduction  current  versus  V_,_  for  a constant 

SCj 

Vds=3V.  This  Figure  exhibits  the  linearity  of  the  FET 
across  most  of  its  dynamic  range.  There  is  a great 
deviation  from  linearity  as  the  device  approaches  cutoff 
(mainly  in  mode  C - VgG  > 3V) . The  variation  of  the  input 
capacitance  with  VgG  for  VDS=3V  is  given  in  Figure  5.5(a). 

The  capacitance  is  seen  to  be  a monotonically  decreasing 
function  of  VgG>  The  drain  conductance,  g^,  is  plotted 
in  Figure  5.5(b)  versus  VDg  with  VgG=lV.  It  is  seen  to 
be  quite  large  for  low  VDg  (high  low-field  mobility) , but 
it  decreases  sharply  above  the  knee. 

For  comparison,  the  I-V  curves  of  the  FET  were 
computed  using  the  Shockley  model  (see  chapter  1) . The 
results  are  presented  in  Figure  5.6.  They  were  computed 
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Figure  5.4  I-V  curves  of  the  FET  (a)  In  vs.  VDg 
(b)  ID  vs.  VSG 
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Figure  5.6  I-V  curves  for  the  FET 
using  Shockley's  model 
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| 

2 

using  the  low  field  mobility  of  GaAS  (4250  cm  /V-sec)  as 
the  constant  mobility  for  the  Shockley  model.  The  Shockley 
pinchoff  voltage  of  this  device  is  4.88V,  thus  the  currents 
predicted  are  very  high.  Of  course,  comparison  can  be  done 
only  for  VDg  below  the  knee  in  Figure  5.4  («  0.4  5V). 

Comparing  Figure  5.4  and  5.6  in  the  range  0 < VDg  < 0.45,  the 
conclusion  is  that  the  Shockley  model  predicts  somewhat 
higher  current  than  the  model  presented  here.  This  is  due 
to  the  smooth  V-E  curve  used  in  the  new  model  compared  to 
the  constant  mobility  in  Shockley  model.  Also,  the  existence 
of  the  transition  region  in  the  new  model  at  the  expense  of 
the  neutral  channel  tends  to  predict  lower  currents  compared 
to  the  two-region  approach  taken  in  Shockley's  model. 

The  measured  I-V  curves  of  a device  similar  to  the  one 
used  in  the  computations  are  given  in  Figure  5.7.  This 
device  has  the  parameters  as  listed  in  table  5.1  except  for 
an  epitaxial  layer  thickness  of  0.32  microns.  The  slightly 
larger  thickness  causes  the  saturated  current  to  be  somewhat 
higher  compared  to  the  computed  results  in  Figure  5.4. 
Comparison  of  Figure  5.4  and  5.7  yields: 

1)  The  predicted  current  is  slightly  lower  than  the 
measured  current,  as  explained  above. 

2)  The  predicted  knee  is  lower  than  the  measured  knee. 
This  fact  is  due  to  the  fact  that  the  external  voltage  in 
Figure  5.7  is  compared  to  the  internal  voltage  in 
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Figure  5.4.  To  get  the  external  voltages  from  the  known 
internal  voltages  the  voltage  drop  across  the  source  and 
drain  regions  should  be  added: 


’SG  ’SG 
ext  int 


VSG  " RSS  XS 


(5-1) 


V =V  + R I+V 
VDS  VDS  SS  AS  vdomain 

ext  int 


(5-2) 


Vdomain  the  volta9e  drop  across  the  high  field 
accumulation  region  between  gate  and  drain  (see  discussion 
in  section  1.3.6  and  also  references  [14] , [15] , [16] , [20] , 
[21]).  The  characterization  of  this  region  is  not  included 
in  the  new  model.  Suggestion  for  future  work  on  this  topic 
is  outlined  in  chapter  6. 

3)  The  measured  characteristics  include  a negative 
conductivity  region  for  VgG  around  zero  and  VDg  near  the 
knee.  This  effect  is  also  due  to  the  high  field  domain. 

5 . 5 SCATTERING  PARAMETERS 

The  most  common  method  used  to  characterize  a microwave 
linear  two-part  network  is  by  its  scattering  parameters. 

When  the  FET  is  used  as  a small  signal  device,  such  as  a 
low-noise  amplifier,  it  is  considered  a linear  device  and 
can  be  characterized  by  its  scattering  parameters  which 
relate  the  reflected  (b^)  and  incident  (a^  waves: 
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(5-0) 

I 


The  measured  scattering  parameters  of  the  device 
discussed  in  the  previous  section  (with  0.32  micron 
epitaxial  layer  thickness)  at  bias  VgG=2V,  VDg=6V  are 
presented  in  Figure  5.8.  S12  is  quite  small  and  ranges 

from  0.01  at  1GHZ  to  0.07  at  10  GHz  (not  shown). 

The  scattering  parameters  of  the  device  given  in 
Table  5.1  were  calculated  using  the  computer  model.  The 
general  approach  employed  in  the  derivation  of  the 
small-signal  parameters  is  a small  perturbation  around  the 
DC  bias.  First,  the  general  equations  for  the  network  in 
Figure  4.1  are  considered.  Then,  each  current  and  voltage 
is  expressed  as  the  sum  of  a DC  value  and  a small  AC  term. 
This  is  inserted  in  the  general  equations  above,  and  thus 
each  equation  is  separated  into  DC  and  AC  equations.  For 
the  small-signal  parameter  calculation  only  the  AC 
equations  are  considered.  In  these  equations  express  each 
variable  in  the  form  of  amplitude  and 
(i.e.  vSGe-Ja,t,  v(jse'3aJt)  • Rearrange  the  resulting 
expressions  into  the  standard  form  used  for  small-signal 
parameters  definition  (i.e.  the  external  voltages 
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expressed  in  terms  of  the  external  current  for  the 
Z-matrix) , and  get  the  expressions  for  the  parameters. 

For  the  network  considered  here  (Figure  4.1)  it  was  found 
useful  to  consider  first  a simplified  network,  and  then 
add  the  additional  elements.  The  procedure  is  as  follows: 

1)  Consider  the  network  as  in  Figure  4.1  excluding 
the  elements  Lg,  Lg,  Lfl/  Csin,  CSQ,  Rgubst.  The  diode  is 
excluded  also  since  in  small  signal  applications  it  is 
always  off.  The  Z parameters  of  the  simplified  network  are 
derived  as  discussed  above.  The  results  are: 


(1) 

Z 

11 


Rg  + 


E-D-C »F 
A *D-B  *C 


(1) 

Z 

12 


A *F-B »E 
A *D-B  *C 


(1) 

Z 

21 


D *G  C*H 
A*D-B*C 


(1) 

Z 

22 


A »H-B  *G 
sd  A »D-B  *C 


where : 

A = - jwGVSG 

B = - jwGVDS 


(5-4) 

(5-5) 

(5-6) 

(5-7) 

(5-8) 

(5-9) 
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C = - gmcos  (ux)+j  (o)*DVSG+gmsin  (ojt)  ) (5-10) 

D = gd  + jwDVDS  (5-11) 

E = - (1+g  Roccos  (ujt)  )+j  [w  • SVSG  (R  +R  ) + 

m So  c ss 

+gmRss  sin(oiT)  ] (5-12) 

F = 9dRss+j«‘SVDS.(Rc+Rss)  (5-13) 

G = g R cos(oit)  + 
m ss 

+ j[“Rss*SVSG+gmRss  sin(cox)]  (5-14) 

H - 1 + gdRss  + jwRss-SVDS  (5-15) 


t is  the  time  delay  (transit  time) . The  Z parameters 
equations  were  first  derived  without  considering  t,  and 
then  instead  of  g the  expression  g was  used,  thus 

yielding  the  above  expressions. 

2)  Add  to  the  network  considered  in  step  1 the 

resistor  Rsubst-  Since  this  resistor  is  parallel  to  the 
output  it  directly  adds  to  * Thus»  calculate  the 

inverse  of  the  Z ^-matrix  to  find  the  matrix,  and 

add  i/Rgu^gt  to  *22^ • yielding  the  matrix. 

3)  Add  to  the  network  considered  in  step  2 the  3 

inductances:  L , L . L . . Since  the  inductors  are  in 
g s d 

series  with  the  terminals  of  network  2 it  is  easy  to 
show  (after  converting  Y^  into  Z^): 
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(3) 

(2) 

Z 

Z 

+ 

jw  (L  +L  ) 

(5-16) 
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(3) 
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jwL 

(5-17) 
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12 
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(3) 

(2) 
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j«Ls 

(5-18) 

21 

21 

(3) 

(2) 

Z 

22 

Z 

22 

+ 

(LD+Ils) 

(5-19) 

4)  Add  to  the  network  considered  in  step 

3 the  2 

capacitors:  C . , 
sm 

Cso' 

thus  yielding  the  final 

network. 

As  in  step  2,  here  too 

the 

! capacitors  admittances  add 

directly  to  , 

Y (3) 
22 

, respectively: 
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Thus,  by  successive  transformations  from  Z to  Y and 
Y to  Z the  parameters  of  the  assembled  FET  chip  are  derived. 
The  Y parameters  computed  by  Equations  (5-20)  to  (5-23)  are 
used  to  calculate  the  scattering  parameters  by  the  standard 
transformation  equations  (see  [24]): 


’ll 


(1~yll)  (1+y22)'>~y12y21 

d+y11)  d+y 


22'  Jr12y21 


(5-24) 


’12 


-2y 
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(1+Y11) (1+y22)_y12y21 


(5-25) 
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21 


(1+y11) (1+y22)-y12y21 


(5-26) 


22 


(l+yll) (1~y22)+y12y21 
(1+yll) (1+y22)-y12y21 


(5-27) 


where  y y12,  y21'  y22  are  normal;‘-zed  Y 
parameters,  namely,  the  Y parameters  divided  by  Yo  (=l/Zo) . 

The  calculated  scattering  parameters  of  the  device 
considered  here  (Table  5.1)  for  applied  DC  bias  of 
Vds=5.8V,  VgG=*2.1V  (internal  voltages  - as  discussed  in 
section  5.4)  are  presented  in  Figure  5.8  for  the  frequency 
range  1-10  GHz.  Comparing  the  calculated  and  measured 
parameters  one  concludes  that  there  is  a reasonable 
agreement  for  (the  increased  difference  in  the  high 

frequency  may  be  due  to  inaccuracy  in  the  assumed 


A 
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VOS *5.800  VOLTS  (b)  CALCULATED 
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Figure  5.8  Scattering  parameters  of  the  FET 

(a)  (b)  S12  (frequency  in  GHz) 
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Figure  5.8  Scattering  parameters  of  the  FET 

(c)  S21  (d)  S22  (frequency  in  GHz) 
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parasitic  elements) , however,  there  are  noticeable 
differences  for  S21,  S22.  The  calculated  |S21|,  I s 2 2 1 are 
somewhat  larger  than  the  measured  results.  These 
differences  can  be  accounted  for  as  related  to  the 
high-field  static  Gunn  domain  between  gate  and  drain,  as 
discussed  in  chapter  1.  The  present  model  does  not  take 
into  account  the  Gunn  domain.  The  main  effects  of  the 
domain  are:  1)  a large  portion  of  the  applied  VDS  bias  is 
dropped  across  it,  thus  decreasing  appreciably  the 
internal  VDg  bias  across  the  active  region  under  the  gate. 
2)  it  presents  an  additional  AC  impedance  (capacitance  and 
resistance  in  nature)  instead  of  the  simple  resistance 
Rgd  in  Figure  4.1. 

The  decrease  in  the  internal  VDg  bias  causes  an 
increase  of  gd  (increased  loading  of  the  drain) , which 
decreases  both  I s22 1 and  | s2i I * This  is  demonstrated  in 
Figure  5.9  in  which  the  calculated  S21,  S22  are  presented 
for  a reduced  bias  VDg=0.8V.  From  Figure  5.9  it  is  seen 
that  S21  for  the  reduced  bias  matches  reasonably  with 
the  measured  results.  S22  is  also  changed  in  the  correct 
direction.  However,  since  the  domain  impedance  is  not 
included  in  the  calculation,  it  still  differs  from  the 
measured  results. 

5.6  LARGE  SIGNAL  ANALYSIS 

The  preceeding  section  compares  the  calculated  DC  and 
small  signal  performance  of  the  FET  with  experimental 
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Figure  5.9  Scattering  parameters  f^r  reduced  bias 
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results.  In  this  section  the  large  signal  performance  of 
the  device  is  examined.  For  this  purpose,  a simple  power 
amplifier  was  analyzed,  which  consists  of  the  FET 
(Figure  4.1)  with  a 50Q  resistive  load  and  an  ideal  voltage 
generator  with  50ft  source  resistance.  The  operating 
frequency  is  2 GHz.  No  input  or  output  matching  networks 
are  employed. 

The  state  equations  of  this  network  were  derived  (for 
details  see  next  section)  and  their  solution  obtained  with 
the  use  of  a standard  library  subroutine  RKGS  (SSP  library 
of  PDP-11/45) . The  waveforms  thus  generated,  were  Fourier 
analyzed  with  the  library  subroutine  FORIT.  The  analysis 
was  done  for  the  bias  conditions:  VgG=2V,  VDS=6V. 

The  results  of  the  analysis  are  presented  in 
Figure  5.10.  Shown  is  the  fundamental  output  power  versus 
the  fundamental  input  power,  compared  with  measured  data. 
From  the  figure  it  is  seen  that  the  model  predicts  quite 
accurately  the  saturated  power  (to  within  0.4db).  In  the 
linear  region  the  predicted  output  power  is  about  1.5db 
higher  than  the  measured  data.  This  agrees  with  the  small 
signal  analysis  which  predicts  higher  J S2l I fc^an 
measurement,  as  discussed  in  the  previous  section.  This 
is  a consequence  of  neglecting  the  Gunn  domain  in  the 
model . 

Since  the  results  obtained  for  the  DC,  small  signal, 
and  large  signal  analysis  are  reasonable  and  fairly  close 
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to  the  measurement  data,  we  may  turn  now,  with  some 
confidence,  to  the  application  of  the  model  to  a practical 
large  signal  network. 

5.7  OSCILLATOR  ANALYSIS 

To  demonstrate  the  use  of  the  model  in  a large  signal 
case,  an  oscillator  structure  was  analyzed  and  its 

i 

properties  evaluated.  The  network  considered  is  suggested 
by  Clarke  and  Hess  (page  241  in  [23]).  The  circuit  diagram 
is  presented  in  Figure  5.11. 

The  network  shown  is  a tuned-gate  oscillator.  The 
feedback  is  achieved  via  the  coupled  coils  which  are 
designed  such  that  M/L2<<1.  This  prevents  loading  of  the 
tuned  circuit  and  maintains  a high  Q,  thus  increasing  the 
frequency  stability,  and  decreasing  the  dependance  on 
device  parameters.  The  network  employs  a self-bias 
arrangement  for  the  gate.  The  combination  of  CQ , R_,  and 
the  gate-source  diode  clamps  the  positive  peak  of  the  gate 
to  source  voltage  to  VQN,  the  turn-on  voltage  of  the 
gate-source  diode  (as  0.5V).  Thus,  the  gate  to  source  bias 
is  automatically  set  to  approximately  the  negative  of  the 
oscillation  amplitude,  which  is  a self-limiting  mechanism. 

An  attempt  was  made  to  design  a 2GHz  oscillator,  using 
the  guidelines  of  Clarke  and  Hess  [23] , and  then  analyze 
it  in  the  time  domain  using  the  FET  model,  starting  at 
t=0  (turn-on  of  VDD)  up  to  the  point  at  which  steady 
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state  is  reached.  The  component  values  chosen  are  (design 
number  1 ) : 

Vdd=6V,  L1=L2=1  nHy,  M=0.2  nHy,  CR=5.5pF, 

RL=lkft,  CB=10pF,  Rb=3K« 

To  derive  the  network  equations  of  the  oscillator,  the 
FET  equivalent  network  of  Figure  4.1  is  used  (with 
Rsubst=  °°  “ ideal  substrate).  Figure  5.12(a)  presents 
the  network  for  VDg  > 0,  and  Figure  5.12(b)  for  VDg  < 0 
source  and  drain  exchanged) . The  network  is  analyzed 
using  the  state-space  approach.  There  are  nine  state 
variables,  which  uniquely  define  the  network  status.  The 
state  variables  chosen  for  the  analysis  are: 


Note  that  the  variables  VgG,  VDg  are  different 
in  the  two  networks  in  Figure  5.12.  In  (a)  they  are 
designated  by  superscript  S,  and  in  (b)  by  super- 
script D.  Since  the  difference  in  VDS  is  just  in 
(VBg^  = - V^)  it  does  not  present  a problem,  and  the 
variable  V^s'  was  chosen  as  the  state  variable  for  both 


p 

I 
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networks.  However,  Vg^  and  Vg^  are  entirely  different, 
and  must  be  treated  as  different  state  variables  in  each 
network.  This  fact  does  not  present  a problem,  since  when 
switching  occurs  between  the  two  networks  at  VDS=0, 

Vg^  = VSG^  * continuity  condition  enables  switching 

from  Vg^  to  Vg^  and  vice  versa.  This  is  necessary  if 
VDg  becomes  negative. 

The  state  equations  for  the  two  networks  were  derived. 
The  first  four  state  equations  are  different  for  each 
network.  Listed  below  are  the  last  five  equations,  which 
are  common  to  the  two  networks: 


© 


dV, 


C +C  C 

(CB+CR>  I4-CBIL2'V1  ^ -V0  R^ 
CBCR+CSi„(CB+CR> 


(5-28) 


r\  dv2 

© - ‘VV^so 


(5-29) 


© 


dv„  CBI4“(CB+Csin) TL2~  Rg  Vl"  Rg  V0 


0 
3t" 


CBCR+Csin(CB+dR) 


(5-30) 


© 


dl 


L1  L2(VDD“V2)-I4V0 


dt 


(5-31) 
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/T  dIL2  L1V0-M(VDD-V2> 
^ ^ ’ Vl-2  ' ' 


(5-32) 


The  first  four  state  equations  for  the  case  Vnc  > 0 

Do 

(Figure  5.12(a))  are  as  follows: 


/T'  dVSG  _ GVDS(I3-Icon)-DVDS  IG 
•O'  dt  GVDS  *DVSG  -DVDS*GVSG 


(5-33) 


^ dVDS  _ GVSG(I3-Icon)-DVSG  IQ 
Vi  dt  DVDS -GVSG-GVDS* DVSG 


(5-34) 


where:  I_  = I.  for  diode  OFF 

G 4 


V ^ ® ) +v 

Ir  = 1,-1 — — for  diode  ON 

o 3 con  R_ 


(5-35) 


The  diode  state  is  determined  by  comparing  the  diode 
voltage  against  Vn  : 


DIODE 


“ VSG  + VWW 


(5-36) 
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(SJ 

0 *_*<  - S,v2-vDs|-(VLd"vi-v0N»*»i'‘-<i”,,-1-,R,d)tl4‘1-,R„-'l-,*‘-D|l"atB,.»'  (5-40) 

d » d q+^»^9 


The  first  four  state  equations  for  the  case  of 
VQS<  0 (Figure  5.12(b))  are  as  follows: 


© 


(D) 

dVSG  _ GVDS(ID-Icon)-DVDS.IG 
dt  GVDS • DVSG-DVDS • GVSG 


dV 


(S) 

DS 


6VS6(ID“ICon,-DVSG,IG 


\ l.  d t GVDS* DVSG  - DVSD • GVSG" 


where: 


JD  “ X4  - *3 


1^  = 1^  for  diode  OFF 


(D) 


V-  r+vr»N 

ZG  * J4  * X3  “ Xcon  " "I—  for  diode  0N 


DIODE 


(D) 

" " VSG  - Rc(I3  + Icon) 


(5-41) 


(5-42) 


(5-43) 


(5-44) 


(5-45) 


(5-46) 
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The  third  and  fourth  equations  are: 

Diode  OFF 


<S) 


as  d y s 9 


(5-47) 


(S)  (D)  (S> 

dI4  . L.  1 VVDS>  ~ 1 VLd>  ‘ VVSC*Rc1con-VDS)  - « I 'L.  '".d*11,.1  * ‘V-d1  ‘V",.1  1 414  'L,R..-  'W  ( V",.1  1 

^ ^ • LdL.**-d'-,*L.L9 


(5-48) 


Diode  ON 


© 


dl 


(S)  (S) 

1 . 'W  (V2-VDS|-L»(VrVUS-V0N|-,]|l>.rt<",»>  ,Lg~>L.)-L,"..l*14|l,,.La-|,gL,) 

^ L.1*-.“-dL,*L.L9 


(5-49) 


© JF1  * 


(S)  (S) 

L,  «''2-VDs'  * <L.*1-dl(Vl-VDS-VOH|-I3|L.(R»d*R..|-|V,-d>"..|414|L.R..-tVtd)l\.<l,a)' 


(5-50) 
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In  the  above  equations  (5-33  to  5-50)  the  various 

circuit  element  values  are  given,  and  the  other  parameters 

(GVDS , DVDS,  R , I , etc)  are  evaluated  by  the  FET  model 
c con 

for  the  given  pair  of  (VgQ,  VDg)  in  the  case  of  VDg>  0, 
and  for  the  pair  (VgG,-VDg)  the  case  VDS<  °* 

The  above  set  of  state  equations  is  solved 
numerically  on  a digital  computer.  The  computer  used  is 
a mini  computer  PDP-11/45,  and  the  library  subroutine  used 
to  solve  the  equations  is  called  RKGS . The  solution  is 
based  on  the  Runge-Kutta  algorithm.  For  the  specific 
problem  outlined  above  the  initial  conditions  for  the 
state  variables  are  all  zero  at  t=0.  At  t=0  the  bias 
voltage  VQD  is  applied  and  the  oscillation  build-up  starts. 
The  computations  were  performed  up  to  t=50nsec  (100  cycles 
at  the  oscillation  frequency  of  2 GHz) , at  which  point  the 
oscillator  has  practically  reached  steady  state. 

Figure  5.13  presents  the  build  up  process  of  the 
oscillator.  Shown  is  the  output  voltage,  Vq,  as  a function 
of  time.  It  grows  from  zero  at  t=0  to  an  amplitude  of 
approximately  7 volts  after  about  80  cycles  (40nsec) . A 
two  cycles  period  of  the  output  voltage  in  steady  state 
(91st  and  92nd  cycles)  is  given  in  Figure  5.14.  It  is 
seen  that  the  output  waveform  is  practically  distortion 
free.  This  waveform  was  Fourier  analyzed  using  a standard 
library  subroutine,  FORIT,  in  the  mini  computer  PDP-11/45. 
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T ME" 5 O00E-10  SECONDS) 

(b) 


Figure  5.13  Oscillator  build-up  (a)  3rd  to  10th  cycles 
(b)  11th  to  20th  cycles 
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Figure  5.13  Oscillator  build-up  (c)  21st  to  30th  cycles 
(d)  81st  to  90th  cycles 
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Figure  5.14  Oscillator  output  voltage  in  steady  state 
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The  results  of  the  Fourier  analysis  confirm  the 
visual  impression  of  Figure  5.14.  The  fundamental  output 

power  is  24.3  mW  (13.9  dbm) . The  second  harmonic  power  is  ( 

0.035  mW  (-14.6dbm),  namely,  28.5  db  below  the  fundamental. 

The  other  harmonics  are  even  lower  than  the  second 
(i.e.  the  third  harmonic  is  -21.6dbm).  The  total  DC  power 

; 

supplied  by  the  DC  source  VDD  is  245  mW.  Thus,  the 
oscillator  efficiency  is  about  10%. 

There  are  two  problems  with  the  oscillator  circuit 
considered  here:  1)  There  exists  a low  frequency  parasitic 
oscillation  in  addition  to  the  desired  oscillation  at 
2 GHz.  This  low  frequency  oscillation  modulates  the 
amplitude  of  the  output  waveform.  This  can  be  seen  in 
Figure  5.13(d),  in  which  the  "peaks"  of  the  output  sine 
wave  are  seen  to  vary  slightly  from  one  cycle  to  another. 

The  modulation  frequency  can  be  roughly  estimated  as 
250-390  MHz.  The  amount  of  amplitude  modulation  is  quite 
small.  The  peak-to-peak  change  of  the  carrier  amplitude 
is  about  0.17  V (compared  to  7V  amplitude).  2)  There 
exists  high  frequency  parasitic  oscillations  in  the  drain 
circuit.  This  can  be  seen  in  the  waveforms  of  VDg  and 
1^,  presented  in  Figure  5.15  and  5.16.  The  high  frequency 
oscillations  do  not  appear  in  the  output  waveform  due  to 
the  filtering  effect  of  the  resonant  circuit.  The 
frequency  of  these  parasitic  oscillations  can  be 
estimated  at  about  10  GHz. 
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The  above  mentioned  parasitic  oscillations  can  be 
reduced  by  decreasing  the  value  of  L^.  This  causes  the 
AC  component  of  VDS  to  decrease  (the  FET  works  into  a low 
inpedance  load) . It  is  also  desirable  to  decrease  the  load 
resistance,  R_  . This  makes  the  structure  more  practical 
for  high-frequency  applications,  in  which  50 ft  systems  are 
utilized.  Thus,  the  new  design  values  chosen  are  (design 
number  2) : 

L1  = O.lnHy  , - 200ft 

This  design  was  checked  by  analyzing  the  network  in  the 
time  domain,  as  described  above.  The  steady-state  waveforms 
are  presented  in  Figure  5.17.  From  these  Figures  it  is 
seen  that  the  low  frequency  oscillation  does  not  exist  in 
the  new  design  (Figure  5.17(b)).  The  high  frequency 
oscillation  still  exists  but  it  is  greatly  reduced 
compared  to  the  previous  design  (compare  Figure  5.15,  5.16 
to  5 . 17 (c) , (d) ) . The  output  waveform  distortion  is  very 
small.  The  Fourier  analysis  results  are: 

fundamental  power  : 18.5  mW  (12.7dbm) 

second  harmonic  power  : 0.0033  mW  (-24.8dbm) 

third  harmonic  power  : 0.0024  mW  (-26.2dbm) 

total  DC  power  dissapation 


428  mW 


tnnroc  z-«  -coc  ohtoc  z*-«  hccc 
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Figure  5.17  Oscillator  waveforms  (second  design) 

(a)  output  voltage  (91st  and  92nd  cycles) 

(b)  output  voltage  (91st  to  98th  cycles) 
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Thus,  less  distortion  was  achieved  with  the  new  design. 

However,  the  output  power  for  the  new  design  is  smaller, 
and  the  efficiency  much  worse  (aobut  4.3%)  . 

Several  other  designs  were  evaluated  in  an  attempt  to 
achieve  better  efficiency  without  low  frequency  oscillations. 

In  addition  to  designs  number  1 and  2,  above,  the  following 
designs  were  analyzed:  j 

design  3 - same  as  2 except:  M=0.3nHy,  1^=500^ 

design  4 - same  as  3 except:  VDD=  5V 

design  5 - same  as  3 except:  VDD=  9V 

In  Table  5.2  the  performance  of  the  various  designs  is 

summarized.  It  is  seen  that  design  number  4 has  the  best 
efficiency.  Design  number  3 has  a slightly  lower 
efficiency  and  slightly  higher  output  power.  All  designs 
have  good  spectral  purity.  In  the  above  simulations  the 
average  computation  time  per  one  cycle  of  the  output 
voltage  is  about  2 minutes  on  the  PDP-11/45.  This 
corresponds  to  about  6 seconds  on  the  IBM  360. 
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6.  SUMMARY  AND  CONCLUSIONS 

In  this  report  a new  AC  large  signal  model  for  the 
GaAS  FET  is  presented.  The  model  is  based  on  basic 
principles,  namely,  electric  field  and  charge  transport  in 
the  device.  The  differential  equations  characterizing  the 
device  are  formulated  and  their  solution  obtained.  In 
contrast  to  previous  simulations  reported  in  the  literature, 
the  solution  obtained  here  is  analytic.  This  is  achieved  by 
introducing  several  approximations,  the  most  important  of 
which  is  assuming  an  analytical  expression  for  the  charge 
carrier  concentration.  The  approximate  analytic  solution 
which  results  is  quite  close  to  the  exact  solution.  However, 
to  improve  accuracy,  corrective  numerical  iterations  are 
employed  as  necessary. 

The  important  effect  of  velocity  saturation,  which  has 
a great  impact  on  the  device  properties,  is  included  in  the 
model.  Device  operation  is  discussed  and  analyzed  in 
terms  of  operational  modes,  which  clarifies  the  physics  of 
the  device  behavior.  The  model  is  "true"  AC  in  the  sense 
that  the  total  terminal  currents  are  computed  as  functions 
of  the  total  terminal  voltages  and  their  time  derivatives. 

The  fact  that  the  model  is  nearly  analytic  makes  it  very 
efficient  and  fast  when  implemented  on  a digital  computer. 
Thus,  the  model  is  practical  for  the  analysis  and  design 
of  networks  containing  FETs. 
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The  usefulness  of  the  model  is  demonstrated  in 
chapter  5 , in  which  the  computed  FET  performance  is 
presented.  The  model  is  shown  to  be  useful  for  small 
signal,  DC  and  large  signal  analysis.  The  most  important 
feature  of  the  model  is,  of  course,  the  large  signal 
capability. 

Summarized  below  are  the  main  limitations  of  the  new 
model,  along  with  suggestions  for  future  work,  which  may 
remove  some  of  these  limitations: 

1)  In  the  present  form  the  me  Ml  can  be  used  for  the 
doping  level  range:  5.1015  cm-3  to  2.1017  cm-3 

(see  appendix  8.3).  This  is  not  a basic  limitation  and 
depends  on  the  approximation  range  in  appendix  8.3.  If 
necessary  the  range  can  be  expanded  by  extending  the  above 
approximation.  Similarly,  the  model  uses  the  numerical 
values  of  the  diffusion  coefficients  of  GaAS . Thus  it  can 
be  strictly  used  for  GaAs  only.  However,  the  model  can 
be  easily  modified  for  any  material,  by  using  the  proper 
values  for  the  diffusion  coefficients. 

2)  The  model  can  handle  devices  with  uniform  doping 
concentration  only.  This  may  be  an  important  limitation 
in  the  future,  since  there  is  a present  tendency  to 
develop  non-uniformly  doped  devices.  (Investigators  expect 
to  improve  the  linearity  and  noise  performance  of  the  FET 
by  the  non  uniform  doping) . 


3)  The  model  assumes  an  ideal  substrate  (namely,  a 
perfect  insulator) . In  practice,  the  substrate  is  not  ideal, 
and  there  is  some  current  conduction  through  it,  resulting 
in  power  loss.  Of  course,  the  substrate  loss  can  be  handled 
as  an  external  resistor,  as  done  in  Figure  4.1.  A better 
representation  of  the  substrate  can  be  probably  achieved, 
when  the  problem  of  non-uniform  doping  level  is  solved. 

4)  The  model  in  its  present  form  does  not  take  into 
account  the  high  field  static  Gunn  domain  present  in  the 
region  between  gate  and  drain  as  discussed  in  chapter  1. 

As  shown  in  section  5.5  in  the  calculation  of  the  small 
signal  scattering  parameters,  the  main  effect  of  the  static 
Gunn  domain  is  to  drop  across  it  a large  portion  of  the 
externally  applied  VDg.  Thus,  the  internal  VDg  across  the 
active  gate  region  is  generally  small  (the  operating  point 
near  the  knee) . This  decreases  the  gain  (S21)  and  the 
output  impedance  (S22) . To  characterize  the  static  domain, 
it  is  necessary  to  assume  a domain  shape  and  match  the 
boundary  conditions  within  the  active  gate  region.  As  a 
first  approximation  it  may  be  sufficient  to  introduce  the 
voltage  across  the  domain  in  the  model.  A better  simulation 
of  the  domain  could  be  achieved  by  considering  also  the 
domain  impedance,  which  is  generally  a combination  of 
resistance*,  and  capacitance  in  parallel. 
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1 The  model  presented  in  this  report  is  a powerful  tool 

for  the  analysis  and  design  of  components  built  with  FETs. 

It  allows,  for  the  first  time,  the  prediction  of  output  power 
levels  for  large  signal  microwave  components  in  practical 
circuits  starting  from  a physical  description  of  the  FET 
device.  Its  use  for  analysis  purposes  is  straightforward  by 
following  the  instructions  in  chapter  4.  The  analysis 
results  combined  with  an  optimization  program  can  be  used  to 
design  networks  with  desired  performance.  Thus,  the 
existence  of  this  model  should  stimulate  much  network  design 
research  in  the  future. 
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8 . APPENDICES 
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ndix  8.1 


Electric  Field  Expressions 


This  appendix  lists  the  expressions  for  the  electric 
field  inside  the  FET  for  each  of  the  three  operational  modes. 
The  expressions  include  the  electric  field  on  both  the 
source  and  drain  side,  and  for  each  one  of  the  regions 
(channel,  transition,  depletion)  as  well  as  at  the  gate. 


Mode  A 


source  side  (x=0) 


channel 


,vi  „ "V2* 
15 


(8.1-1) 


transition 


Ex“  “[Ig(1" 


2ir(a-|  -d^O)  ) 


"(y-d.  (0))  ir  (y-d.  (0) ) 

{ -3 sin  3 }) 


irV  /2a 

+ ^|sroos 


(8.1-2) 


qN_d  it  (y-d.  (0))  it  (y-d.  (0)) 

Ey  “ 2Tre”  < 5T sin  3T ] 


(8.1-3) 
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depletion 


V.  y-d/2-d,  (0)  irV  /2a 

Ex  = ~d > + S-TTZ  cos<»] 

a-y  - (0 ) sinh(j^)  2a 


(8.1-4) 


Y e 


- (y-  | - dx  (0) ) 


(8.1-5) 


drain  side  (x=lq) 


channel 


Vl 

r 1 O 


x " "CTg  + 21“  coth(2a  } cos®] 


(8.1-6) 


Ey  ‘ 2iT  sin(St> 


(8.1-7) 


transition 


Ex  “ -IlfU- 


irty-dj^dg)) 


2ir  (a-j  -d^  (lg) ) 


- sin 


(y-d1  (lg) ) 


>)  + 2iT  coth  cos  (JJ)  (8.1-8) 


qNDd  7T  (y-d.  (lg)  ) 
Ey  * 2ire  ^ 3 


* (y-d. (lg) ) irV^ 

3in  a 1 ♦ .in  (JJ) 


(8.1-9) 
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depletion 


E 


x 


V.  y-d  (lg)-d/2 
“[Tg(1"  a-dinTg)-d/2)  + 


cos(^) 


qN_ 


(y- 


I -d^lg)) 


2a~ 


sin(S) 


(8.1-10) 

(8.1-11) 


gate  (y=a) 


qND  d 

E=Ey*  -^(a-  J -d1(x))  + 


Sinh(ji2) 


sinh  (^J) 


(8.1-12) 


Mode  B 

source  side  (x=0) 
same  expressions  as  for  mode  A. 


drain  side  (x=lg) 


(8.1-13) 


transition 


(8.1-14) 


-185- 


qNDd  V-V  irV 

Ey  = 2iT"t0-sine+  — 3— (9+sin0)]+  ^ sin(g)  (8.1-15) 

depletion 

V1  <3Nnd 

Ex  = ?ef-  tzy-aa- 

(8.1-16) 


V 


1 


+ 


wo 

7T- 


(8.1-17) 


gate  (y«a) 

region  0 < X < y 

same  as  for  mode  A (8.1-12) 

region  y < X < lq 


qNDd  a 

V -zr-'2 f -1 


V 


F~  + Uq] 1 + 


(irVo/2a) 


s inh  ( ^^) 


(8.1-18) 
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Mode  C 


source  side  (x»0) 


transition 


V,  qN_  .2  , ffV  /2a 

■'lg(1+  *eT(F>  {2cos9~9  >>  + cos(j5>> 


(8.1-19) 


Ey  = 2~~  t0-sin0"  (9+sin0)] 


(8.1-20) 


depletion 


vl  , ttV  /2a 

-[Ig(1-  4FiT  (2y-d<1-  ;2>>>+  slnh(Iia)  =°=<2f>) 

2 cl 


(8.1-21) 


qNOd  r2  Z , VP, 

"2T"  [2  d 1 B-1 


(8.1-22) 


drain  side  (x=lg) 


same  expressions  as  for  mode  B 


jate  (y=a) 


same  expression  as  for  mode  B (8.1-18) 
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Appendix  8 . 2 

Displacement;  current  expressions 

This  appendix  includes  a of  the  displacement 

current  expressions  at  the  FET  terming s in  the  three 
operational  modes.  The  expressions  have  be-.n  derived  as 
described  in  chapter  2.  The  following  notations  are  used: 

d,  8 as  defined  in  Equations  2-13,  2-33 


CQ  i ew 


K1 


CD10  = <PB 


0.0237 


^d2 


(8.2-1) 

(8.2-2) 

(8.2-3) 

(8.2-4) 


K1  and  K2  are  functions  of  VSQ,  VDg  and  are  different 
for  each  operational  mode.  They  are  calculated  as 
described  in  chapter  3.  The  sign  notations  are  gate  and 
source  - outward,  drain  - inward,  thus: 


(8.2-5) 


Therefore,  only  the  gate  and  drain  expression  are  given 
here.  The  source  current  is  calculated  from  8.2-5. 


i 

\ 
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Mode  A 


Define:  KGS  (^V1+VSG+CD1°  " ^vSG+CDl°) 


KDS  = [coshf^S.)-  l]/sinh(ji2.) 


ulc 


KG 


, iia  {^- 

N 2e  V,  l^vl 


x .+VSG+CD1°  ‘ Wl  [(Vl+VSG+CDl0) 


3/2 

- (Vsg+CD10)  ]}  - KDS 


KD  * coth(ji^)- 


. 0 . 1487d  +0 . 5a  (a-d) -d.  (lg) [a-0.5d.  (lg)-0. 

15  «• 1 1 


a(a-d1(lg)-d/2) 


+ CSG1A 


eV,  ta-H  (lg)]  +d (d, (lg) -a+0 . 2d) 
2qNDlg 


[a-  I " d1(lg)]3 


(8.2-6) 


(8.2-7) 


3/2 


(8.2-8) 


5d] 


(8.2-9) 


CSG1A 


(8.2-10) 
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then, 

gate: 

drain: 

Define: 


KGSC 

KGS  = 

+ 


IDIS=Co[(KGS+KG‘K1)~3Fi  +(KDS+KG.K2)  (8.2-11) 

dV__  <jv 

IDIS-Cot"(KD'K1+CSGlA)_at-  +(coth(5i2)-KD*K2)  -5^] 


Mode  B 


(8.2-12) 


4 


KDS  = same  as  equation  (8.2-7) 


qN 


mB  “ ^ Ig[1~  Of  d2l  - 

a-d  qNDd 

_ Ig~fl"  4£3“te+0.2026d)  J 


(8.2-13) 


(8.2-14) 


[2qN  

V ~e  vf  1 ^ ^p+VSG+CDl°  ‘ 


+CD10  ] + 


KGSC  (1- 

V1 


(8.2-15) 


then, 

gate: 

drain: 
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■^IS^cJ  (kgsC+KGC»K1)  d^G  + (KDS-f KGC  • K2 ) ■ ^°S] 
IDIS!=Co,:"KDC’K1*"3t“  + (coth  (^2.)  -KDC  *K2) 


(8 


.2-21) 


(8.2-22) 
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Appendix  8 . 3 

Ratio  Function  For  Mode  C 

This  appendix  contains  the  numerical  calculation  of  the 
ratio  function  needed  for  the  source  conduction  current 
approximation  for  mode  C (see  section  3.4.1).  Also  presented 
is  a polynomial  approximation  to  the  above  function. 

The  function  considered  here  is  given  by: 


II 

I"  f (l+cos8)d8  ~l 
L“'  Jl+  (Ev/Ex)  2 -I 


/ [8+sin8] 


1+ (Ey/Ex) 


where , 


2 


8-sin8+K  (8+sin8)  2 

£ 2 1 

i+Ryl  (2cos0-(j  ) 


qNDd 

2ireE 

c 


(8.3-1) 


(8.3-2) 


(8.3-3) 


K^1  is  given  in  Equation  3-4 
8 is  given  in  Equation  3-46 

Kyl  is  a small  number  (typically  0.02-0.03)  and  was 
found  to  have  negligible  effect  on  the  function  value  in 
the  range  of  the  other  parameters  considered: 
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0 < Kc  1 1 

5 . 1015cm-3  < ND  < 2*1017cm~3 

The  values  of  the  constants  used  are: 

q=1.6*10“19coulomb,  E =3 . 2 • 105V/m,  e=12.5e  . 

C - o 

The  computation  of  the  integral  in  8.3-1  was 

performed  numerically  using  100  partitions.  The  results  are 

presented  in  Table  8.3-1  for  11  values  of  K and  6 values 

c 

of  Np.  The  same  data  is  also  presented  in  graphical  form 
in  Figure  8.3-1.  It  is  evident  that  the  function  considered 
is  a monotonically  increasing  function  of  both  Kc  and  . 
Therefore,  it  is  reasonable  to  assume  that  a polynomial 
approximation  can  be  found.  An  attempt  was  made  to  find  a 
fifth  order  polynomial  approximation  of  the  form: 

R (K  ) =RR0+RR1 *K  +RR2  *K2+RR3  *K3+RR4  *K4+RR5  *K5  (8 . 3-4 ) 

c c c c c c 


L9 


Table  8.3-1  The  ratio  function 


.2  0.3 


.5  0.6 


0.8  0.9 


5.1015 

1.095 

1.132 

1 1.207 

1.32S 

1.473 

1.632 

1.1016 

1.129 

1.210 

1.356 

1.572 

1.814 

2.063 

3.1016 

1.170 

1.325 

1.617 

1.954 

2.279 

2.571 

6.1016 

1.219 

1.459 

1.883 

2.308 

2.672 

2.986 

1.1017 

1.230 

1.559 

2.065 

2.525 

2.909 

3.222 

2.1017 

1.263 

BBS 

2.378 

2.895 

3.295 

3.616 

3.262  3.443  3.606 


4.500  4.664 


0.1  0.2  0J  0.4  0.5  0.8  0.7  0.8  0.9  1.0 

Kc 


Figure  8.3-1  The  ratio  function  (N^  as  a parameter) 


I 
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The  coefficients  RRO,  RRl,...were  determined  by  curve 
fitting  the  polynomial  to  the  data  in  Table  8.3-1.  This  was 
done  separately  for  each  value  of  NQ.  The  curve  fitting  was 
achieved  by  using  a standard  IBM  optimization  routine  named 
FMCG.  This  routine  is  included  in  the  Scientific  Subroutine 
Package  (SSP) . The  resulting  values  of  the  6 coefficients 
for  each  ND  are  presented  in  Table  8.3-2.  The  values  of  the 
coefficients  for  any  intermediate  value  of  ND  are  calculated 
by  linear  interpolation. 

The  polynomial  approximation  was  found  to  be  very 
accurate  (within  2-3%  of  the  numerical  solution) . In  fact, 
the  graphs  of  the  polynomials  are  practically  coincident 
with  the  graphs  of  the  actual  function  in  Figure  8.3-1. 


Table  8.3-2  Coefficients  for  the  polynomial  approximation 


~ t 

ND(cm"J)j  RRO 

\ 

RR1 

RR2 

RR3 

RR4 

RR5 

5.1015 

1.095  ! 

] 

0.066 

2.824 

-1.123  | 

-1.510  1 

s 1 

I 

| 0.982 

1.1016 

1.129 

0.364 

4.621 

| 2.065 

3.1016 

1.170 

! 1.464 

1 

4.992 

i-3.833  I -3.162 

L _ I i 

! 

’ 2.982 

6.1016 

1.219 

2.652 

j 

4.481-4.894  | 

1 ! 

i 

| -2.806  j 

! j 

3.362 

1.1017 

’ 

1.230 

1 

3.740 

i 

3 . 138 | -5 . 429  ! 

J 

1 1 

J -0.957  j 

2.531 

1 

2.1017  I 1.263  j 5.702 


-0.443-3.179 


-0.559  1.895 
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Appendix  8.4 

Integral  Function  For  Mode  C 


l This  appendix  contains  the  numerical  calculation  of  an 

integral  function  which  is  part  of  the  expression  of  the 
drain  conduction  current  for  operational  mode  C - case  2 
(section  3.4.2).  Also  derived  is  an  analytic  approximation 
for  the  above  function.  The  function  considered  is 
ir 


I 


/ 


(1+cosQ ) d6 
Jl+ (Ey/Ex) 2 


(8.4-1) 


where: 


0-sin8+£3 (0+sin0) 


1 ^coa^') + tan«2911 


(8.4-2) 


ir2eV 


L1  qN„da  coth(?S2') 


(8.4-3) 


'2  * la 


(8.4-5) 


*1  vp 

(8.4-6) 

It  is  obvious  that  I+O  for  Also  for 

Ey 

■*  tan(-f20)  and  the  integral  becomes: 


I 
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ir 

/ 


(l+cos0) cos (£29)d0 


(8.4-7) 


‘1  o 


This  integral  can  be  evaluated  analytically: 


I — ^1^2^  * X-  s^n(^2ir^  + 2T1V& — J”  sinI  (1+^2^  ‘n’]  + 


J^-*® 


+ 2 (1 -l2)  sint (1"£2) ^ 


(8.4-8) 


The  function  considered  here  was  evaluated  numerically 
(100  partitions)  for  the  following  range  of  parameters: 


0 < *1  < 5,  l2  = 0.15,  0.4  , 0 < l3  < 1 

The  results  are  given  in  Table  8.4-1,  8.4-2  and 
Figure  8.4-1.  Also  included  is  the  value  of  (Equation 
8.4-8),  which  is  the  limit  for  This  limit  is  shown 

also  in  Figure  8.4-1. 


1 
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Table  8.4-1 


Table  8.4-2 


o 1 2 3 4 5 6 7 

A 


Figure  8.4-1  Numerical  evaluation  of  the 
integral  for  l2  = 0.15 


-201- 


From  Figure  8.4-1  it  is  evident  that  the  function 
considered  here  is  similar  to  an  exponential  with  a limit 
of  y^  for  The  slope  depends  on  Jl-j.  Thus  an  attempt  was 

made  to  approximate  the  function  by  the  following: 


I 


=s  y1  (*2)  t1  “ V*2  (*3)  e 


_Jcl*l 


(8.4-9) 


y^  is  given  in  Equation  8.4-8  and  for  simplicity  y2 
is  assumed  a linear  function  of  l^z 

y2(£3)  - 1 + k3  • Z3  (8.4-10) 


The  three  coefficients  k^,  k2,  k3  are  chosen  by  trial 
and  error  to  yield  a good  fit  to  the  tabulated  data.  The 
goal  was  to  achieve  good  curve  fitting  especially  for 
> 1/  since  as  discussed  in  section  3.4.2,  in  practice 
^ is  rarely  below  1.  The  coefficients  chosen  are: 

kx  * 2.414,  k2  = 0.3835,  k3  ■ 4.798 

The  function  in  8.4-9  was  evaluated  for  the  same 
range  of  parameters  as  the  actual  function  considered. 

The  results  are  tabulated  in  Tables  8.4-3,  8.4-4. 
Comparison  between  Tables  8.4-1,  8.4-3  and  8.4-2,  8.4-4 
shows : 


A 


I 
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Table  8.4-3 

Evaluation  of  the  function  in  8.4-9  for  £2=0.15 


Si 

0 

0.1 

1 

2 

5 

0.0 

0.0 

1.955 

2.819 

2.964 

3.061 

0.3 

-4.456 

0.313 

2.421 

2.772 

3.010 

0.6 

-8.913  | -1.330 

2.022 

2.581 

2.960 

1.0 

-14.855 

-3.520 

1.490 

2.326 

2.892 

Table  8.4-4 

Evaluation  of  the  function  in  8.4-9  for  &2=0.4 


as 

0 

0.1 

1 

1 

2 

5 

0.0 

0.0 

1.787  ] 2.577  j 2.709 

2.798 

m 

-4.075 

! 

0.286  ] 

2.213  | 2.534  j 

2.750 

wm 

-8.150 

-1.215 

j 

1.848  j 

1 2.360 

2.705 

1.0 

-13.583 

-3.217 

1 

1.362  | 

i 

2.127 

2.644 
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For  tj  * 0 the  approximation  is  reasonably  accurate 
for  all  vaues  of  A^,  in  the  range  considered. 

For  Aj  71  0 the  approximation  fails  for  A^  < 1 . 

The  accuracy  is  improved  the  lower  A3,  the  lower 
l 2 and  the  higher  A^  are. 
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